关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
( V1 [9 y6 M  G  Y0 d  dh=b-a;
n=1;
eps=1e-3;
T1=0;
% ]/ J+ r: `& ^! _8 r" xT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
2 _) M2 _7 Y$ C  K5 P% W; M    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;   
6 s+ @5 s  u0 P) R' C    end
    T2=1/2*(T1+h*s);
6 {  |- n4 e: X    n=n+1;
end
I=T2+1/3*(T2-T1);
' j  `# X! g8 I! V' ^2 h 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
' O$ }' s: k" a3 o0 R) Va=a+eps;
! `7 q2 C# L! X/ r5 Wh=b-a;
8 R8 ^2 ~8 C, o$ M3 w# cn=1;
ep=1e-6;
7 e2 B! Z  A- `3 l7 W) A. ST1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
' S% Y+ d% o6 c8 \3 etol=1;
# |% f4 I! Q* }# G; `6 J) xwhile tol>ep
    u=h/2;
    x=a+u;
    H=0;
: Q; x% W4 J! p+ |2 k    while x<b
9 R6 H0 a8 u/ J+ R        H=H+subs(sym(f),findsym(sym(f)),x);
7 ^% k# y# {" S0 y  Z/ c9 S' K) }        x=x+h;   
    end
8 a/ I/ G" x: Y/ S; u: X7 Y    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
' Q5 W% ?# z% V+ w, j    h=u;
    T1=T2;
    n=n+1;
end
' {+ G6 V8 D: ]1 E0 i  w' I% KI=T2+1/3*(T2-T1);
# L# V3 \  u. E* E. s4 I