关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
& v; H0 I+ e) _- b) w1 Mh=b-a;
3 C. d9 A; J0 W- V; m' v, \3 x0 Mn=1;
eps=1e-3;
1 b  K' U- l! q0 z7 k1 {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;
    x=a+h;
- u! \5 M( E0 J8 {$ h- y( k3 X    s=0;
$ ?9 P5 k$ B. _* E2 h$ H! r+ f    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
% `1 ~$ a( ]2 S- \5 a- v        x=x+h;   
+ S0 F/ \6 p6 s( v0 g4 f    end
    T2=1/2*(T1+h*s);
. W7 i0 B3 {4 P: `    n=n+1;
end
I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
9 F$ Z* F! Q3 ~3 [6 da=a+eps;
h=b-a;
n=1;
ep=1e-6;
1 y' `6 Z/ S8 ^2 ^- wT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
. C) U6 @( S; w+ y, w0 ~8 Ctol=1;
8 I0 m( p$ N3 c6 F& ?& F! ^: x9 }while tol>ep
    u=h/2;
4 `7 X7 i% n$ ~7 y2 \4 ^4 P    x=a+u;
    H=0;
7 u/ y5 m2 {3 {% T) ]- S    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
  \+ f- q9 N9 G        x=x+h;   
  [0 |0 X' a" h" h% E! p6 M    end
+ g0 k# c2 g. Y    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
" x. ^# c6 w4 u& v9 e  |  `, z    h=u;
    T1=T2;
    n=n+1;
end
I=T2+1/3*(T2-T1);
! ]7 U$ k6 E0 Y! q8 A6 M- N; K