关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
- K; g  g% s' x/ m3 d( Q( w8 |h=b-a;
4 y% m, @+ x  hn=1;
) O9 X8 {( q- N. @- `eps=1e-3;
T1=0;
5 j0 G; ~/ D) k. ]' UT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
- M! j3 Z1 V4 K6 D" ^( R9 x! vwhile abs(T2-T1)>eps
    h=h/2;
    T1=T2;
    x=a+h;
' p* N/ Y: {7 J" ^8 [    s=0;
5 m; `/ P$ P$ h! _6 H! d, k    while x<b
& x: h. k8 G" y4 z: |        s=s+subs(sym(f),findsym(sym(f)),x);
# _: B4 B1 E1 A( x7 _3 R        x=x+h;   
    end
) ?9 ?' A9 x+ A1 I    T2=1/2*(T1+h*s);
    n=n+1;
2 U; w; I* m9 N3 Q* _) k6 vend
% O2 {2 o# B' @/ `I=T2+1/3*(T2-T1);
- M, ~# t2 R/ D 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
a=a+eps;
h=b-a;
7 m, n: z* \% e9 k) On=1;
ep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
+ j* b5 i. M4 s& V5 L+ W2 Owhile tol>ep
$ x% S% C; F/ v. X8 J/ J    u=h/2;
    x=a+u;
7 l4 }4 U/ t" ]3 T1 ^' C' M    H=0;
4 h3 g7 v# J4 e, e$ @8 b    while x<b
. q* ~/ W% V- C% v/ t6 X1 V9 q        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
4 S$ T) x9 L' f$ U3 E' w7 S& B    T2=1/2*(T1+h*H);
5 Q6 h& m9 h4 e$ E  j    tol=abs(T2-T1);
    h=u;
# m! [3 V% t: ^0 H' z( u9 O9 A( s    T1=T2;
    n=n+1;
; F1 f& }. O& g# w  bend
I=T2+1/3*(T2-T1);
: u( T) p( Z' M0 @