关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
n=1;
/ ^& `6 b  j$ v6 xeps=1e-3;
T1=0;
' j* t/ c, n. P( E% @  T! FT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
* A2 N3 W1 s, B    h=h/2;
    T1=T2;
  v/ S7 [- ~  _$ A6 l$ t" _+ p% X4 n2 q    x=a+h;
    s=0;
1 a+ q  X1 r, |% v- A  w3 T    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
2 _: d* r% l8 {) C        x=x+h;   
    end
$ b" ?# O* h$ E/ d3 D0 ?    T2=1/2*(T1+h*s);
    n=n+1;
$ W6 p" z$ B/ o, {end
: }( I/ D. ^) Q+ b7 JI=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
& r$ Y) R, u9 t9 M1 V5 |5 rfunction [I,n]=FHTX(f,a,b,ep)
/ l7 b1 W4 I' Z4 I9 Aa=a+eps;
h=b-a;
3 G; e( i2 p# V( l1 }/ Mn=1;
4 G! I8 H" m; `* Eep=1e-6;
) K0 |, c# N/ W; q4 rT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
2 b) Q/ I( T4 T  dtol=1;
1 Z9 d3 G9 }2 b2 K7 q. G# T/ gwhile tol>ep
    u=h/2;
    x=a+u;
    H=0;
# o" f; G. s: B! S    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
9 q% N- W; g1 ]  K% o* z# X    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
! K$ c; {- ?) T: R5 ^% H+ i5 a    h=u;
    T1=T2;
    n=n+1;
end
I=T2+1/3*(T2-T1);