关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
' K% s5 E/ v0 T$ a7 a4 H, c3 }5 M5 zh=b-a;
n=1;
8 y9 q  Z6 B- Y% P. M. Q, b# jeps=1e-3;
+ r4 q  s" \- g7 u0 I  l- WT1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
$ f# c6 a# g" s5 N, G: a: _9 p+ Bwhile abs(T2-T1)>eps
    h=h/2;
    T1=T2;
    x=a+h;
    s=0;
    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
+ D9 r( b! ]9 b+ K  f        x=x+h;   
: |3 @7 D7 \6 I# f  k1 y' x    end
    T2=1/2*(T1+h*s);
    n=n+1;
end
% c3 l7 {3 y: z& p4 ?% H4 ~0 J! V9 @I=T2+1/3*(T2-T1);
3 c# r* e" {4 J1 o1 W7 J9 R 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
8 S1 [) [8 O# x7 \2 @' g. |a=a+eps;
h=b-a;
% T: u4 {4 ~7 l% Jn=1;
2 i$ k: q6 ]1 Y$ k1 U, gep=1e-6;
. e; }- L3 {6 WT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
9 }! F$ S# K% S2 k  D+ l, Ttol=1;
while tol>ep
    u=h/2;
    x=a+u;
" L8 |. l0 n/ a    H=0;
$ f) H2 t' K: U* R    while x<b
8 W7 |- |4 t$ \- o% N! f' ~8 i5 l        H=H+subs(sym(f),findsym(sym(f)),x);
, C* w; D/ w( x        x=x+h;   
    end
! r3 [' J3 t' h& D6 H$ ?* b( g    T2=1/2*(T1+h*H);
, y' s$ D$ v) H" r9 F    tol=abs(T2-T1);
    h=u;
+ J: q( }4 f* ?/ d% A; _5 F    T1=T2;
    n=n+1;
4 a* T* q1 @) V7 m* P2 J8 ~end
4 ~' x9 M6 b: KI=T2+1/3*(T2-T1);
; R& P: N( {2 t& \( T