关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
" H8 y! ?1 m" d% [- jh=b-a;
8 C& i$ c3 W+ f% f6 |; W7 v* X! c% fn=1;
eps=1e-3;
7 h& ~2 D* j+ j4 }T1=0;
8 x) J: G% }$ g/ RT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
' g2 Z5 }: B7 j" S0 v    h=h/2;
    T1=T2;
) r4 [6 T) C) [5 f1 P    x=a+h;
$ Y8 I2 v/ d, V* G) ]4 ]' S* r( _    s=0;
    while x<b
6 t2 V. o6 }1 N' k- h* w        s=s+subs(sym(f),findsym(sym(f)),x);
6 B3 s, |  ?8 j" A6 r5 C5 |        x=x+h;   
    end
    T2=1/2*(T1+h*s);
) O( O+ V5 ^  |7 T4 }9 O    n=n+1;
end
I=T2+1/3*(T2-T1);
7 d4 C; G4 V  T* p3 R- u9 a 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
: B& Z1 i3 F6 H2 Y- oa=a+eps;
: q! X7 u. u% Z6 {h=b-a;
n=1;
ep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
) ?7 C6 A+ q+ j' p$ _while tol>ep
/ G) \1 B0 J* k$ w* e# ^    u=h/2;
+ u0 E/ z( P. M1 F, X: H- B( C* z    x=a+u;
    H=0;
    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
9 E& j6 C* p: O5 n' i        x=x+h;   
" C: f. c  H" ?" K    end
1 W& B1 G  m# u$ ~# J4 A0 x    T2=1/2*(T1+h*H);
0 e  l5 N+ Z5 M# ^& O    tol=abs(T2-T1);
    h=u;
    T1=T2;
    n=n+1;
end
7 j0 y4 I7 r, v3 A; R2 p1 AI=T2+1/3*(T2-T1);