关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
3 L/ j9 }; n2 X5 Q, Sh=b-a;
5 m, u$ H) j# {; U) H' tn=1;
" R7 I- _* H6 h/ reps=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
; G: z. r1 Q1 e0 ]7 v    h=h/2;
& ]3 l* R& ?, M3 B1 n* \    T1=T2;
# |( u( k0 |% g, u3 C& H9 ]  J& ]    x=a+h;
7 L" U. N2 m! c% Z/ `    s=0;
( ?6 F9 T& S" Z& r/ J. u- |    while x<b
  t3 N8 V6 q; e* b        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
! M: j0 X+ g- Z    end
    T2=1/2*(T1+h*s);
. u8 C  s5 T1 H3 K8 \+ [& _    n=n+1;
end
I=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
% E, L! R/ n- t' |/ j1 `2 G0 v1 Lfunction [I,n]=FHTX(f,a,b,ep)
a=a+eps;
4 N- _- a, f, [# u2 r1 Rh=b-a;
n=1;
4 S' n' ], B& Y% Zep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
! ?' i3 p& V& H) g+ c6 vtol=1;
while tol>ep
4 r7 d4 }# S, i; D9 `% U    u=h/2;
2 M  J" L' m! j8 X    x=a+u;
    H=0;
    while x<b
6 Q' B- p! w+ x- p0 K: ?+ Y  ~; V        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
8 R: d: \' t# d; D: n    h=u;
    T1=T2;
    n=n+1;
2 O+ e! I. g0 g% j9 L  [end
0 S  _% B, Q2 i7 ]1 S8 M. oI=T2+1/3*(T2-T1);
6 T- ?- O. ^( a' J