关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
/ v: v5 a0 c, j* J1 P$ V$ zn=1;
eps=1e-3;
: ~. M% ~! [: [9 w9 ~9 j1 eT1=0;
7 b* r+ z" V% _- x0 L/ ?5 sT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
7 g- [: k* I$ r( _' U' awhile abs(T2-T1)>eps
8 }4 b2 S% N3 P! T! F7 N, ?) ?+ ]    h=h/2;
6 Z0 l4 O. u! N6 N+ Q    T1=T2;
/ y  \  g: w5 a+ H* `    x=a+h;
4 }- B  A4 T5 m; _. g& m6 C) `1 o    s=0;
    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
. _# R0 A7 B: Z' r9 o        x=x+h;   
8 M, q& o4 C8 n    end
) O# ^9 Y  K+ o# {    T2=1/2*(T1+h*s);
/ Z$ @; Z, U& ]( K0 [; F# M    n=n+1;
end
  @% H& o& _1 O, s  y2 HI=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
* p3 n) n+ @9 w) Q3 ua=a+eps;
h=b-a;
! q2 W  w; v" k5 L- G, o& ln=1;
8 W  v3 X& m  d$ @/ l' u. K. Xep=1e-6;
5 R9 I. P' ^9 f' t) sT1=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;
* T. P7 h6 P3 q# u- {1 C6 `+ M# D: |2 ?    x=a+u;
    H=0;
! m& m/ a6 l5 N; b# j& C, Q# i5 j! M6 b    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
0 z% K9 \" x" }" e- ~0 _    end
5 X, u$ `$ _3 j0 C    T2=1/2*(T1+h*H);
  L6 S( D6 a! Y4 w/ a; @3 }. m    tol=abs(T2-T1);
4 z0 J- n4 j: L& t! k/ i5 B    h=u;
3 A/ Z' x2 v$ B4 G    T1=T2;
    n=n+1;
- N' Z2 {! U1 H5 D( A5 Oend
6 ^/ F  W1 P" A1 j6 HI=T2+1/3*(T2-T1);
1 z8 B5 Y/ M2 S