关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
! W7 K* H! e" g( S5 R9 _: lh=b-a;
9 \1 _: H1 V# wn=1;
eps=1e-3;
T1=0;
" e/ \# y0 z2 d$ G2 r8 cT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
* E6 y2 g+ }- Y4 k5 w    h=h/2;
    T1=T2;
$ O4 U/ l( }+ G" f* T6 Y    x=a+h;
    s=0;
    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
( A) c& f, m- w2 N    T2=1/2*(T1+h*s);
; P- g1 s( ~6 c0 }* B    n=n+1;
. l; ]$ q, F. `8 rend
I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
  L; V  V, Z+ f0 N/ x* K. Kfunction [I,n]=FHTX(f,a,b,ep)
- a! p& {- W- D0 Ka=a+eps;
$ s/ `" y2 H* Jh=b-a;
) s9 y* _% Y' x. Yn=1;
5 t# [/ r+ e$ y. Y4 Tep=1e-6;
8 w% `  k. i' _4 ~; ], FT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
, _) v2 h6 i: ^, z$ r8 E- O" stol=1;
4 c. y- H: I2 t$ S- ^. vwhile tol>ep
    u=h/2;
/ A& e; k, ^! v" l- Q5 h" J: V/ r    x=a+u;
    H=0;
' ^' C8 o! t7 z7 a7 l    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
    h=u;
. f1 q1 k9 X2 _  ]    T1=T2;
    n=n+1;
end
; A4 B$ a5 G9 A! D" p: Z4 qI=T2+1/3*(T2-T1);
3 d% E6 k1 D: p