关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
# n  w( o& Q( [7 {2 q0 f; sn=1;
# }! L' H  T4 Seps=1e-3;
% q# B( K- L: ?: ~- b5 V5 t$ C+ jT1=0;
" g9 a8 k" \2 VT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
) S/ K* f& {: ]4 mwhile abs(T2-T1)>eps
    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;   
+ r6 o4 D) y: C8 ~    end
    T2=1/2*(T1+h*s);
    n=n+1;
end
+ \; J% f6 J9 o+ }) cI=T2+1/3*(T2-T1);
7 `  g* g: g5 |& ~, o3 A 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
( Y( z* a) ^* U! cfunction [I,n]=FHTX(f,a,b,ep)
a=a+eps;
/ b! R- }+ q5 @7 ]0 ?4 l4 J. ~. Wh=b-a;
n=1;
ep=1e-6;
3 R  q5 Z+ Z; \& rT1=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;
    H=0;
. S5 k1 ]/ l3 S4 _- n% [    while x<b
9 Q% D* b: K% H9 j2 ?& q$ U        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
4 t* N7 n; g: V% s! K3 g    end
: ^0 M1 P/ i* e* z    T2=1/2*(T1+h*H);
3 P7 M' I3 s. e2 r* n- i: }/ Q    tol=abs(T2-T1);
; o  Y# u( h- [- L: v# M    h=u;
) N, G6 Y& p$ k8 c% \9 h! d    T1=T2;
    n=n+1;
end
& {9 V( u& J! `) f. q9 C# LI=T2+1/3*(T2-T1);
+ W/ \6 _; Q: |7 }- g