关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
8 w8 j9 X$ W0 m3 T* i0 T/ hn=1;
) \' @# p5 _' U/ t- Ueps=1e-3;
* J. C! |" F( sT1=0;
* G4 o' r, e" b. T2 [$ U, W/ dT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
. z( @  \0 [! M* C; I# f    h=h/2;
- f: I  t/ @0 @( F' @9 h    T1=T2;
    x=a+h;
6 u1 J# H2 x& I# r7 w- l' ?    s=0;
, H; n$ P8 P, x# i9 g3 A    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
9 o/ A2 E% u( f2 Z( z        x=x+h;   
, \3 x0 b; ?. C+ z. Q4 u9 F    end
. g1 U& P3 r2 c# a) h4 S    T2=1/2*(T1+h*s);
, b/ {4 F* i0 T) G: v# `    n=n+1;
end
8 h/ ~. v' C. u( G! e0 }I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
4 L" ~( U* Q1 N. |* Vfunction [I,n]=FHTX(f,a,b,ep)
4 i& F" O6 H  h- U9 Z- s# i  ^a=a+eps;
+ m6 m& n' d- J0 P# c: j3 B  Kh=b-a;
n=1;
ep=1e-6;
: F- C8 P# S/ I# K" b& Q" `) V# g- mT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
$ f: l" J: D/ q( M7 R. _6 atol=1;
while tol>ep
1 w( w5 s% z5 n5 E/ G- I; |    u=h/2;
9 l. q7 _9 z2 F" ]/ ?    x=a+u;
4 x- `) r. N' D% ]' T) a- L    H=0;
1 @: K1 \; N$ T: {6 L    while x<b
% m) e9 Q( D. u, X, P. F) N        H=H+subs(sym(f),findsym(sym(f)),x);
& C5 g. ~$ n( P  }. f7 s        x=x+h;   
    end
! _4 w: j( D: I  d; p( D* a% ^    T2=1/2*(T1+h*H);
2 \2 G/ H. L! a/ M% L& o    tol=abs(T2-T1);
  m" s4 l% F2 e1 P    h=u;
' F; L. k4 _/ ]/ C1 t+ t! E4 J4 f    T1=T2;
    n=n+1;
+ q0 F! C" Y* e: K) v6 y8 ~; hend
I=T2+1/3*(T2-T1);
4 Y  K3 A, b0 R/ j$ W5 L, i% I# p