关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
" a+ W2 H5 @! \  H3 U& V( K' `h=b-a;
$ T  V/ D/ n, L, f3 f' ]n=1;
eps=1e-3;
1 c2 h& \( H* Q/ a7 ^" LT1=0;
6 ?+ Q  @1 D8 s: _T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
8 M% I1 Y. U; r1 ?/ u3 |( {/ wwhile abs(T2-T1)>eps
    h=h/2;
    T1=T2;
    x=a+h;
    s=0;
    while x<b
9 F/ C2 }1 Y; v4 ~8 c, A        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
    T2=1/2*(T1+h*s);
4 U/ S0 m5 C) I$ H- v# d9 l. d$ W    n=n+1;
end
9 h2 Q; `0 K( }3 }* O1 bI=T2+1/3*(T2-T1);
* N" [7 i" o# {( Z4 ^6 J+ v 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
1 p8 W7 s3 i; p7 c% Wfunction [I,n]=FHTX(f,a,b,ep)
' D) c4 `- F- K' w5 f% r$ \4 Na=a+eps;
$ E5 f& F, d  Oh=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));
, u8 R' S: f; Q: R% \tol=1;
% S* d9 a- J  V! Q, w! D# ~8 `, Wwhile tol>ep
& Z* {% C$ K# q' q    u=h/2;
. Z- k5 r0 e4 V* R$ f% {1 {" N    x=a+u;
- ^" Z' c  i( }, ^    H=0;
9 _. z) ~! ?% K& b9 Q; Y6 F    while x<b
& U% L4 \( ^4 }, {5 G5 t2 a- U2 z        H=H+subs(sym(f),findsym(sym(f)),x);
7 h" A4 a! k0 T/ U4 d; l        x=x+h;   
, }; Z+ F# l2 M1 e! k5 V' k+ e" p1 L    end
    T2=1/2*(T1+h*H);
, c9 ^! T' o. W1 g. S2 Y0 }    tol=abs(T2-T1);
/ L* H3 I% N& E% w9 c/ _1 t4 P+ }    h=u;
    T1=T2;
    n=n+1;
( u3 y' n" b! mend
" t9 A: E0 e0 t; ]/ Q7 `I=T2+1/3*(T2-T1);