关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
5 Q' ?0 K% i- M8 p) C6 `5 G. Nn=1;
- p* \, q# n# {eps=1e-3;
1 k1 M4 A: Q/ K- A+ w. D+ i/ ~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;
    T1=T2;
0 _, x3 u$ \+ X% Y: S7 B6 P; ~8 E    x=a+h;
    s=0;
! h& p! w2 Y) w# b+ n0 N( b( x    while x<b
/ |8 m7 t5 R. {; S& e  B# @' L5 P) `; p        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
    T2=1/2*(T1+h*s);
/ q4 n' ~+ [% p& h/ i+ u" g" Z    n=n+1;
6 B6 G7 s. m  _; Zend
I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
& P0 X' g+ h& o6 c. _# hfunction [I,n]=FHTX(f,a,b,ep)
3 D) O: h0 D3 @0 ba=a+eps;
! ]" v' ^6 N" e2 Z8 ~6 P: Sh=b-a;
n=1;
ep=1e-6;
% C& B* u9 q& g5 X% c- UT1=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;
& i5 }- I/ D- u9 w: S& ?* b+ V    x=a+u;
    H=0;
4 A! W% |& M8 l5 X. {/ a. I    while x<b
% o- N. M  u- }6 U3 O        H=H+subs(sym(f),findsym(sym(f)),x);
/ S& u4 q( ?( j        x=x+h;   
/ I) |( U1 c: Z+ W' x' y    end
    T2=1/2*(T1+h*H);
6 b* o! U3 A& d3 P2 u    tol=abs(T2-T1);
' g! g% \3 f% r# Y/ O' L    h=u;
5 V9 t8 Q4 j1 R; Z# K    T1=T2;
    n=n+1;
2 O, f: b- E4 Kend
I=T2+1/3*(T2-T1);