关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
% w) [; |8 P+ T. Th=b-a;
n=1;
eps=1e-3;
, U- j% k# s$ @! n( qT1=0;
3 [1 \! U5 J$ j1 L6 k0 J' BT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
7 ?$ |" d5 A& H9 @9 W; |* |+ `' jwhile abs(T2-T1)>eps
    h=h/2;
2 a  z5 n: N, A# a: b; ^1 F    T1=T2;
    x=a+h;
2 x6 X4 B8 o3 l1 O$ N    s=0;
( Z7 u% x/ N; s  c    while x<b
! Q/ d* c3 u8 u        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
    T2=1/2*(T1+h*s);
) n) f1 _+ }: h, Z2 }6 }, V    n=n+1;
. K- A0 [5 p/ A& Send
: G; u; o& F4 g/ j6 `* XI=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
a=a+eps;
) n  o* p9 d* a' ~' A$ @( Yh=b-a;
n=1;
- r/ [% y1 H; nep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
2 Y; `8 C0 N9 ^tol=1;
while tol>ep
    u=h/2;
3 y" m( ]. t0 R1 \    x=a+u;
    H=0;
7 j/ y$ b/ H7 k; S/ o    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
$ S1 d# M0 D9 }/ ]; \0 g        x=x+h;   
    end
    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
    h=u;
) o- s4 e7 ]- S; m; ^  f, M# r5 G    T1=T2;
# j4 r6 z7 P* S* b# ^; b    n=n+1;
end
, V2 H) a  d, n6 I7 N( |. a0 cI=T2+1/3*(T2-T1);
9 d6 d, O6 V  ~  a' w