关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
7 D# U7 c& x' `$ Qh=b-a;
n=1;
eps=1e-3;
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
$ ^, P8 e7 }# ]& c2 h  C    h=h/2;
1 g) l  \+ j4 Z    T1=T2;
( Y% I9 A4 Q; w+ e3 F# @7 z    x=a+h;
    s=0;
1 Z" w, t* `4 s& `( A4 M: I    while x<b
% y% F0 m. m: Q/ |        s=s+subs(sym(f),findsym(sym(f)),x);
0 ~9 b' v$ I/ `        x=x+h;   
    end
( _+ ^8 t: U" c: |  R    T2=1/2*(T1+h*s);
    n=n+1;
end
I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
; M% Y. q4 y0 i6 O1 `" s& U2 `function [I,n]=FHTX(f,a,b,ep)
* [- Z: _5 \* @0 h& Ea=a+eps;
! g/ l! g2 B: R6 X" g% r, Uh=b-a;
6 S8 e8 V6 d5 W7 v& T- z# T) j  x9 y* f/ Wn=1;
" i: @& e: \2 d: Kep=1e-6;
/ U8 `2 Z' @9 R' ?' l/ Z" l- }* KT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
0 l: N$ D/ r, u- P7 Ctol=1;
. y$ m# J  N* j: N# Owhile tol>ep
, Q9 \( T' Z# q# |+ o2 e    u=h/2;
    x=a+u;
    H=0;
3 S+ b/ o7 I- @/ L) l1 W    while x<b
( {$ @" l# U5 p+ g/ |        H=H+subs(sym(f),findsym(sym(f)),x);
7 T( A) N: H8 R& o4 U& O        x=x+h;   
    end
( k& J2 F7 c; P( T: H9 A    T2=1/2*(T1+h*H);
% {( B, G" U! S7 n$ P    tol=abs(T2-T1);
    h=u;
    T1=T2;
3 k! [: x  t6 v9 t- u; }5 b    n=n+1;
2 G7 t- p# Y1 N" }, O/ r) E+ ]$ Z$ vend
1 x1 Y/ a1 f* sI=T2+1/3*(T2-T1);