关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
! O1 q1 C1 F$ m* h, hh=b-a;
/ g8 Q. f0 v! A+ X, I" m1 Wn=1;
% {6 f9 M2 Z, T; J0 [/ w! Neps=1e-3;
T1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
+ |: L- N8 r8 `/ bwhile abs(T2-T1)>eps
8 m9 D1 D6 b. ~, p. A- }+ E    h=h/2;
    T1=T2;
    x=a+h;
    s=0;
    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
6 H8 d, q  @$ ^1 \    T2=1/2*(T1+h*s);
3 p7 ]- T6 `* ]! N4 d    n=n+1;
end
I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
. h7 e7 B- o3 c" Z- @9 Nfunction [I,n]=FHTX(f,a,b,ep)
" O% C& k1 z8 O# x# O4 j* va=a+eps;
h=b-a;
( h5 m6 ~4 a1 ?6 B; Y3 a6 J( Kn=1;
: n/ ]  W' E. Eep=1e-6;
7 K# ]; f$ G4 g' J# NT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
5 @, S# m0 G8 t. ^2 Ctol=1;
while tol>ep
' x- a# o: J: ~    u=h/2;
" o2 w# ~$ _. J, w3 k( i" M    x=a+u;
! h( Y, |: ^" @0 z) t0 @7 `1 M    H=0;
5 o% b5 ~" L/ T; B( v& _, X    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
" e' u# d, G4 `1 U0 R        x=x+h;   
    end
    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
0 X* s: ?3 }. k3 X. a$ h& O8 n/ G    h=u;
) R- p) U( L, G: S! P    T1=T2;
, w- w7 c& X# V! y8 d9 V* A* @    n=n+1;
end
7 {3 e6 E. m+ GI=T2+1/3*(T2-T1);
; ]8 _9 O* \# Q2 H$ f/ [6 G8 J