关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
; p1 m0 T: z8 kn=1;
  W) `  B* J9 u9 F: K  V- O0 _eps=1e-3;
T1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
- c2 p$ S9 X* Z4 Y9 z3 Z1 P( bwhile abs(T2-T1)>eps
    h=h/2;
    T1=T2;
* n7 O+ {9 ^, ]7 B1 L. h. H4 r    x=a+h;
    s=0;
    while x<b
& a4 x9 R! q& H& F  G" j/ i        s=s+subs(sym(f),findsym(sym(f)),x);
* |7 k$ o+ h0 U8 b. ^) [& c        x=x+h;   
2 P( K7 {; h; t( _2 B, u6 P    end
    T2=1/2*(T1+h*s);
    n=n+1;
end
( {  a  C- O+ m- c! d/ xI=T2+1/3*(T2-T1);
( t, ?; D) O1 ?/ K' e; E9 ?2 M 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
3 [3 [/ g1 L" w( Mfunction [I,n]=FHTX(f,a,b,ep)
' t- d# E& C6 J( }( a5 ga=a+eps;
, U( ~# F8 R( G  {6 {h=b-a;
) b$ \3 H& @/ c9 h+ u& Tn=1;
: f+ J* x/ M3 f2 n$ R% Q8 K- R& W0 {ep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
: F3 N. f; n; E+ M# F4 Ztol=1;
while tol>ep
8 r( d( u8 ?8 l8 M    u=h/2;
# w2 ?5 C+ n1 u1 H3 G; E0 H    x=a+u;
    H=0;
  j7 b! f$ q1 o# {% }1 u: }# ?    while x<b
. J: ]4 ^' P5 q; e. K        H=H+subs(sym(f),findsym(sym(f)),x);
2 J1 u7 O' w, F) Z        x=x+h;   
5 D/ ]/ {1 M; ]2 |; C  K& }  `8 c    end
    T2=1/2*(T1+h*H);
# j7 w2 |0 ]; g& J$ B9 r" j% Y    tol=abs(T2-T1);
1 f: x8 [; s$ v: I+ e    h=u;
5 \  N- t+ m8 V8 G3 W# |+ z    T1=T2;
( m1 s& g" f: H" M    n=n+1;
end
& a& U  g6 v$ [# j' TI=T2+1/3*(T2-T1);