关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
# p7 T' x; f0 L5 }' d+ E& vn=1;
: L6 ]: ^6 \1 Y$ o( c. `eps=1e-3;
8 G- s9 @. D; h. T- `T1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
    h=h/2;
    T1=T2;
    x=a+h;
) |+ ?, z+ u! g. h  B7 h' p" Q    s=0;
  N' l4 m2 Y0 ^, s0 l9 d    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
% M& z( S& R' G1 N/ ^1 }& ?2 t8 N. }        x=x+h;   
5 B, v! R0 }5 u8 ~; F2 Z% @/ @    end
    T2=1/2*(T1+h*s);
3 |( O, k2 X5 ~! D5 Q    n=n+1;
& _9 T; Z3 [# I' U4 Kend
I=T2+1/3*(T2-T1);
) z1 v! R& A. o+ T 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
5 h2 j  T5 Z- h- Dfunction [I,n]=FHTX(f,a,b,ep)
a=a+eps;
7 E' P: t0 L/ U# n5 v7 Kh=b-a;
" h6 s& s3 Y% i+ g# dn=1;
ep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
while tol>ep
    u=h/2;
    x=a+u;
3 o0 E- l0 R" a/ X    H=0;
: D# {4 c+ o7 R    while x<b
5 V. S' }/ C4 ]6 E* @$ Z3 f& c1 U        H=H+subs(sym(f),findsym(sym(f)),x);
4 f9 ^: k% `' o7 q! N        x=x+h;   
    end
" n4 Y9 x' K1 k5 b; G' f    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
    h=u;
    T1=T2;
    n=n+1;
+ l+ o9 D; q2 l( V; X( |6 }6 P5 }end
3 O2 C* b+ ~7 z8 N2 a' B- z. {I=T2+1/3*(T2-T1);
7 t4 f: R& G: r