关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
; f6 F# y6 v! `2 L5 u; B2 v3 {6 D  ^h=b-a;
& d" ^1 M; D* }( fn=1;
eps=1e-3;
! ^" U% `, ?+ R2 |8 RT1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
- g& Y/ H6 S: t- o: X( @+ Jwhile abs(T2-T1)>eps
; v, A5 ]' k: r" K! m    h=h/2;
& q5 Y2 ^" F; b; X4 \% Z    T1=T2;
    x=a+h;
    s=0;
    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
3 w' i9 X/ A' f% Y5 _2 p    end
    T2=1/2*(T1+h*s);
    n=n+1;
+ K, G9 L: r. L' J7 m; Q. hend
$ f" E  E# @8 C1 gI=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
a=a+eps;
& O: W1 G" \5 b+ R0 M2 o# Qh=b-a;
n=1;
- y) a" A9 q0 A& [! @) g, A0 z& n0 }ep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
# j# d9 u* }& C4 ~5 i+ e2 w* ^) J1 lwhile tol>ep
: y" w: O( X( ?8 v$ U    u=h/2;
    x=a+u;
    H=0;
! z& R+ C: _% }" @. k+ E    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
6 O8 @0 y" s. w- p% G/ q' l# [7 p$ h    end
) n3 x% g) k) e/ u    T2=1/2*(T1+h*H);
4 A+ B2 U5 |5 c    tol=abs(T2-T1);
% X' ]/ X) G$ C' n% J3 ^& l    h=u;
1 U- U6 y6 q) [! \    T1=T2;
6 A# F$ _' T3 ^3 M    n=n+1;
end
I=T2+1/3*(T2-T1);
) H1 p4 w$ j1 ~