关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
- W% ?$ Z/ K9 Q0 F6 {6 ~; ih=b-a;
% A# w/ N8 Q( L* D2 zn=1;
eps=1e-3;
4 t1 t. `8 o3 S9 K, v' sT1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
% i& ]9 U" V* d* z* Kwhile abs(T2-T1)>eps
0 M8 x2 @* {* P$ B& r" E* R4 D    h=h/2;
' e) R0 D4 `! G- K    T1=T2;
    x=a+h;
    s=0;
; N! Y) h- p3 p% W    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
) V2 _( r4 G9 l        x=x+h;   
7 t- h1 I# b( a8 D1 q  M9 w$ Q3 D    end
( f; W! O6 v. I! o    T2=1/2*(T1+h*s);
    n=n+1;
  s3 m9 g, C; j# |4 }/ Fend
9 F# @, v$ Y4 b: P( P1 ^: nI=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
$ S* Q, t' i! [5 i% ]a=a+eps;
8 G; w0 ^# d% S. U6 T6 J9 }h=b-a;
n=1;
ep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
# G( \& A' r$ c+ a% v2 M' c! ~  ltol=1;
( V! J' E7 J7 D0 K% Y/ Dwhile tol>ep
    u=h/2;
% X9 a# C  c: m* k. b% ^: c  r    x=a+u;
    H=0;
2 `! h  X; g; g/ J    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
5 G8 H% N: s9 W1 r        x=x+h;   
    end
; w$ Z- o! s; ~    T2=1/2*(T1+h*H);
# A6 b" E9 J& l/ c0 x# J    tol=abs(T2-T1);
    h=u;
    T1=T2;
    n=n+1;
end
I=T2+1/3*(T2-T1);
* p$ {; o8 B+ Q, K# o1 @% a