关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
n=1;
' G1 d- e8 x. r  Q- ?eps=1e-3;
8 q# p& ~' a( d# F! _& {T1=0;
8 b% Q2 V2 K" W. O; pT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
5 N3 \6 {/ L7 Y    h=h/2;
    T1=T2;
    x=a+h;
; c: e$ x3 _5 ?7 k9 a" L; j    s=0;
    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
1 [% b/ P' r8 [0 w    end
    T2=1/2*(T1+h*s);
    n=n+1;
end
I=T2+1/3*(T2-T1);
: w' G, h  G9 K: z+ [ 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
" {% m6 F0 f- C# V/ cfunction [I,n]=FHTX(f,a,b,ep)
a=a+eps;
; ]1 `% N' G- E3 Vh=b-a;
* c( m+ ]0 j) N' in=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
    u=h/2;
6 Z- Q' f/ B( @, R# e" a6 t' @    x=a+u;
    H=0;
    while x<b
* l- e0 c* j6 q' Y5 S' t        H=H+subs(sym(f),findsym(sym(f)),x);
; G9 T# Y4 n9 e1 H' \- O7 k: o# A        x=x+h;   
( i( G8 S$ }" e' b4 O    end
( }* ?. X8 Z4 ^4 B. M: k4 o    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
    h=u;
    T1=T2;
    n=n+1;
end
; W% u4 T3 C; F: Y5 B1 xI=T2+1/3*(T2-T1);
, y- [/ U# f  q- P7 U+ I