关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
2 t1 k- Y* Z) C/ Y" dn=1;
eps=1e-3;
. k/ m& u# R/ N' @6 T6 dT1=0;
! d; T- j" N1 A6 o$ S6 ^T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
% K( T, E# n( ~# e* Pwhile abs(T2-T1)>eps
    h=h/2;
+ B2 X2 |6 K$ Z    T1=T2;
& K0 B9 g' i, w, S, i- ^    x=a+h;
& G2 M2 H" [- R. N7 j    s=0;
    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
! f5 m4 T- @* V. k- q) D- N    end
1 i' ~2 Q1 R' [+ H+ a# o$ R    T2=1/2*(T1+h*s);
    n=n+1;
' H8 y2 _8 w+ s, C- B; Cend
I=T2+1/3*(T2-T1);
. v* w( w* Q; S) j! Z  O2 A6 ^3 } 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
. i5 e. w' [2 \0 t; G) mfunction [I,n]=FHTX(f,a,b,ep)
a=a+eps;
h=b-a;
4 X% w2 [' _, \, n% B- T& ^( {n=1;
1 `2 ?4 X; U, F8 G9 c/ ]ep=1e-6;
3 X3 L2 a; v! Q2 oT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
: `6 X) L) Y/ B/ a. G  J; Gtol=1;
while tol>ep
# ~4 p. g& j3 s/ \' o    u=h/2;
: M" ~; |8 h) ?# p2 H    x=a+u;
! d1 d9 e0 H% ]; `% q; T& j; S2 M    H=0;
$ s' P4 u5 E* Z4 L6 G2 J    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
4 x$ [3 q$ ]7 ?# u0 U        x=x+h;   
    end
5 U9 l  I$ w: H% D1 l5 g* p    T2=1/2*(T1+h*H);
7 T+ K2 V6 t, G3 j. G    tol=abs(T2-T1);
    h=u;
    T1=T2;
6 x5 U) E: L* V* J2 ~    n=n+1;
5 M# j6 D' F3 F7 ~end
4 p1 M7 P7 a$ L! \I=T2+1/3*(T2-T1);
! s0 O! J) m( _& K