关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
' ~" u( p+ S2 W3 xh=b-a;
n=1;
: ^/ L1 `2 M% I% }2 y& heps=1e-3;
3 c7 Q! k! U, G% N& M" {T1=0;
2 V1 _: K7 }/ ^/ k5 C& IT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
/ y& N; }" \( m5 B    h=h/2;
    T1=T2;
    x=a+h;
% c" ?6 W, R( B& h    s=0;
! s0 Q6 s- z* Y3 F7 \    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
1 K( }( F( K9 m, R  {" P. J    end
; D( ~7 D6 F" N    T2=1/2*(T1+h*s);
    n=n+1;
4 j1 i# L% a% n1 F5 N# Xend
I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
6 D; y; k4 z9 r2 z( @function [I,n]=FHTX(f,a,b,ep)
1 M- A" l- \6 M$ Na=a+eps;
: |  S: ~6 M/ J7 x) c& gh=b-a;
n=1;
3 w- u! ]! e$ @* l% `6 rep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
4 B& G2 m& t2 F+ G2 ]tol=1;
) N& m9 _3 _: x6 {5 Iwhile tol>ep
    u=h/2;
5 O% |$ F+ L6 X( c+ ~    x=a+u;
9 e5 f; d2 W* e4 M    H=0;
    while x<b
4 V/ O. l! B, \6 u  V        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
  l6 r7 _# n1 e. D1 K2 C/ @3 {2 J    end
& O' J$ _% b. O& C2 t: Q) l    T2=1/2*(T1+h*H);
0 {7 p, d0 l1 e; K$ ]3 w    tol=abs(T2-T1);
- f7 W7 _+ j" \3 o" c! v: ^    h=u;
" w# B4 Q, a1 E3 M* T    T1=T2;
    n=n+1;
end
/ B, [4 V" T& R( c$ C5 ZI=T2+1/3*(T2-T1);
: o- i7 Q% F% E1 P- W: ~* E0 C4 F* E