关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
  r6 T6 o+ h/ Y3 mh=b-a;
/ N2 P/ p6 P3 E, x/ u! |+ Bn=1;
4 G# U- e) V/ x" X, `eps=1e-3;
: [/ r7 w  S! G& ?0 V- y0 pT1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
# @* p1 x  _0 M    h=h/2;
0 R- Z$ B  t2 [2 j: m0 ]& P    T1=T2;
    x=a+h;
    s=0;
    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
9 ?3 Z2 Y* K2 d8 x& n        x=x+h;   
3 X0 D; ~: f( T1 g! |    end
6 D9 I( z2 C9 ^3 A# T& ^5 H& v    T2=1/2*(T1+h*s);
    n=n+1;
2 B' Z4 L- I# W/ B$ @# o0 G5 fend
I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
- q2 t2 \2 S3 _function [I,n]=FHTX(f,a,b,ep)
a=a+eps;
h=b-a;
+ u  {+ E9 D# G; \8 h2 Cn=1;
ep=1e-6;
* h: p3 s: D( e7 A; F5 h3 ST1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
2 H0 W$ w' P6 V; ttol=1;
# T* f6 T5 l( p% o! Gwhile tol>ep
    u=h/2;
( t+ `& J. S8 K0 h    x=a+u;
    H=0;
! \9 R9 @- @1 J2 {" u: P0 i    while x<b
: F9 ?) ?8 P  [        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
0 K% p0 E  p; {( R" t2 o! C3 u    end
0 g1 ]& ~1 o) B  H7 z    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
    h=u;
* G8 a* l+ O8 ~9 J0 M7 {6 f: l) x    T1=T2;
    n=n+1;
3 u9 J5 M. j' S; n; eend
I=T2+1/3*(T2-T1);