关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
; r9 e  t4 U: Dh=b-a;
n=1;
. b4 b, g0 P! a% X$ Heps=1e-3;
+ G6 l9 U( A: b' \( ]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
8 l3 T- k5 E! Q2 D    h=h/2;
# g" q9 z8 t( q8 L/ P! V% L' \    T1=T2;
    x=a+h;
    s=0;
8 C/ l( z; v* ^, ?3 \9 Q3 t/ a    while x<b
' u7 k; x1 |8 g* E6 R! h$ @        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
2 N( r. [  F2 Z/ Q( U; T    end
    T2=1/2*(T1+h*s);
) s: t$ O# |4 X8 D6 `6 h    n=n+1;
end
3 O( \. f6 M1 Q( g4 }7 U$ t$ LI=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
  h6 [; |+ m# W/ T" `  R' x* \function [I,n]=FHTX(f,a,b,ep)
a=a+eps;
# ~; J2 ^/ I, i* |/ H1 zh=b-a;
: g" T, B1 D7 g& dn=1;
ep=1e-6;
- ^2 [5 p3 Y+ I% _. |T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
0 p' U- [1 ]1 e# k; Ttol=1;
while tol>ep
! c/ x% S6 ~7 g& l. R1 ^6 ~2 V    u=h/2;
    x=a+u;
    H=0;
    while x<b
/ ~$ j3 v# b9 h2 }/ L9 q; D7 r        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
1 _" F  R$ O" u3 j    end
( E  S: ~& ?! H! i3 `    T2=1/2*(T1+h*H);
1 q' l$ K( t. U& R0 K    tol=abs(T2-T1);
    h=u;
    T1=T2;
    n=n+1;
1 M( Q" T% I, T7 K9 _6 d" Kend
I=T2+1/3*(T2-T1);