关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
& V! J1 `) |' {; E1 k: u  yh=b-a;
+ R7 g( N8 Y3 `4 q  c" x2 Tn=1;
5 h' u. ^6 B+ o4 H: h3 x# H0 jeps=1e-3;
T1=0;
8 r1 D) Q$ s2 P# b* }T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
0 N# Q/ N5 h6 ^/ O1 {: |6 [. G" {while abs(T2-T1)>eps
! M6 H! ]! B9 a: I3 W- B. n    h=h/2;
6 g# |; w! c% Q+ r9 E/ _6 w    T1=T2;
3 |" {4 _7 Z2 `: `4 W; j    x=a+h;
, X5 A6 f% @& W* ^' \( b3 k    s=0;
    while x<b
# ^) Z( T7 P1 U+ w        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
8 e2 \0 t4 p( z# ~# Z    end
/ @) K6 @( Y$ W    T2=1/2*(T1+h*s);
    n=n+1;
5 g% w# g; _8 c0 h# J# h1 ]end
: N0 t: X# }' F9 C* mI=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
% f5 E1 {, u6 T' xfunction [I,n]=FHTX(f,a,b,ep)
$ ^/ f& a2 h( V1 D4 Ca=a+eps;
h=b-a;
! N* N: U* S( jn=1;
ep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
while tol>ep
9 l# B4 }' u5 h. I; N% ?0 ~    u=h/2;
    x=a+u;
& z" l) x0 d! l9 I2 `* P    H=0;
) I! {, m4 A, n3 q' V    while x<b
4 x. E  D; U/ z) z& L        H=H+subs(sym(f),findsym(sym(f)),x);
7 k' |8 f; o% ]- E* J$ Y. ^        x=x+h;   
0 X) w0 C! C, L+ s2 Z! V0 b' D" p    end
    T2=1/2*(T1+h*H);
0 ]1 M0 r5 ~; E& U% h2 l, m  X    tol=abs(T2-T1);
: a" ?3 I& F; |1 ^# a+ [0 n    h=u;
    T1=T2;
7 ~$ H6 M/ F8 X+ W  t7 d* w    n=n+1;
end
. g* W- j* M7 \& n6 E! m: kI=T2+1/3*(T2-T1);