关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
$ _( t' V/ C# N# K! uh=b-a;
0 |& D6 o, P# m$ L7 C. jn=1;
# p7 w7 s! Y; ]6 u  E; {eps=1e-3;
) n, ^  s6 P0 j& b8 ^T1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
    h=h/2;
# h: {' S0 I, q' m2 M0 v    T1=T2;
6 L. l. h1 T# l! v& _  _    x=a+h;
    s=0;
    while x<b
# Z2 |# I4 a3 |5 Y6 O8 ?1 Y) s        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
  \3 `7 w6 @9 R, \+ \    end
    T2=1/2*(T1+h*s);
! P$ y4 V5 {( |7 w- F    n=n+1;
- j1 Q. K: W8 T  ?: U. Pend
/ o' N3 `( k( T, f, Y7 L1 _I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
a=a+eps;
9 O' o, I  [1 T) _) E1 X! _h=b-a;
0 K4 s* R2 X; b0 W1 ?n=1;
# L' m, m; _# n; Q& U* n0 X9 Iep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
; j- U2 f4 b- F9 H" q! S9 @while tol>ep
    u=h/2;
    x=a+u;
. g# H: A' i" x3 q    H=0;
    while x<b
; D8 S- x. d7 x$ N        H=H+subs(sym(f),findsym(sym(f)),x);
$ h  m) F$ E/ R3 [6 T/ }: H        x=x+h;   
    end
    T2=1/2*(T1+h*H);
' r% y3 Y6 |. ]: x- T) J" }' U4 N& l% A    tol=abs(T2-T1);
    h=u;
    T1=T2;
& M: r2 n- E1 R9 I% x+ G! I    n=n+1;
end
4 W. O  w, \9 K8 }2 TI=T2+1/3*(T2-T1);