关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
2 o+ E8 n' A3 m. @8 uh=b-a;
2 c5 x& H' d/ n" Rn=1;
4 b2 E% {; Q' A9 i/ s$ ueps=1e-3;
2 S/ v- e8 h6 D2 c5 kT1=0;
( l) k: w' o. R0 IT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
5 T9 d: T- S6 t' n- \0 N" B    h=h/2;
    T1=T2;
    x=a+h;
# Z7 O4 D6 E8 O2 {0 D4 h3 [+ `    s=0;
    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
( v' T# s' W" P        x=x+h;   
    end
' K1 P3 D+ ^; D! \* g/ \4 g    T2=1/2*(T1+h*s);
    n=n+1;
( |; g$ E' Q  s! o7 a3 a0 M5 uend
I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
1 R+ p- Q) m" }/ J$ U& `a=a+eps;
, x: v  O; x  R/ ?1 th=b-a;
n=1;
ep=1e-6;
6 Y( N0 S- c8 s( }- R! x0 V5 b* Y9 w9 nT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
while tol>ep
. f9 j. T  q4 i2 N/ f9 Q. X    u=h/2;
, P+ E. N- l: A1 `; \* a9 x    x=a+u;
& \: H) e0 o: Y2 ]3 K% b    H=0;
    while x<b
8 G" J. A- s4 t9 O1 A$ l. E        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
3 z  Y2 B' O- y, J6 o    end
. ]5 `- T# K2 a    T2=1/2*(T1+h*H);
4 P4 C# _9 s  B, S0 N    tol=abs(T2-T1);
$ N$ z1 J+ t1 [. p1 h9 [    h=u;
; Q( B" q! \( k+ U    T1=T2;
    n=n+1;
end
. V6 r: G2 U  y& \, V* P: B0 K7 s( iI=T2+1/3*(T2-T1);