关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
. h: O) f; K5 Zn=1;
7 `! A/ U3 ]! Q( jeps=1e-3;
T1=0;
; G1 H! B) Y  w$ D7 A( yT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
; ~. v" A9 X# {. n5 J3 I& y( T) Zwhile abs(T2-T1)>eps
( b% j7 G. m" q0 Q2 t  w, ~1 C0 z    h=h/2;
; F: [, t+ [$ }. u9 g; x    T1=T2;
    x=a+h;
    s=0;
    while x<b
5 P- f6 c3 c$ l& t  W" @        s=s+subs(sym(f),findsym(sym(f)),x);
# Z, _/ k- Q! ^" y6 [2 S! I        x=x+h;   
7 v2 I9 l5 D: Q0 e  e1 v    end
8 T2 O, s& G1 n    T2=1/2*(T1+h*s);
  X7 p4 d1 H) D( P2 P* G    n=n+1;
end
I=T2+1/3*(T2-T1);
& n8 ]+ h# Q- X" b0 u 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
a=a+eps;
1 r+ Z. ?+ l8 J: l2 kh=b-a;
n=1;
) c8 Z8 z7 f1 E9 Y0 k' L/ [* T: eep=1e-6;
2 P" Q+ q0 @" gT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
& z8 q& Q0 e/ H' e9 _8 v' k0 ~tol=1;
( A7 |6 ~2 {1 Qwhile tol>ep
* y. Q; U/ T- }( p% T4 ?    u=h/2;
& |, u! h* n8 y. I* K4 T- e+ v' d    x=a+u;
    H=0;
    while x<b
$ w# T9 [* _. s$ ^7 r0 i, _        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
+ @; {5 h/ {( H' ~2 J& j" N    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
8 e# R+ X) {# r. {& W$ ?/ q    h=u;
. b% N+ K3 \5 j' Z0 B4 _( M' O1 X7 ~    T1=T2;
    n=n+1;
/ W5 K" |. {8 X: @/ d6 Xend
% ^' K: \3 C) G8 wI=T2+1/3*(T2-T1);
1 q* O! N* {( U% a! Y- N2 u. \