关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
1 e& |; \+ s. N. An=1;
. Y6 K4 W: K5 j' aeps=1e-3;
T1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
' A& x, q% A9 z3 @% Wwhile abs(T2-T1)>eps
    h=h/2;
    T1=T2;
    x=a+h;
" v6 M% h. Y3 g    s=0;
, u- T, N3 w3 L* T& i% A    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
, Q- ~2 x) B2 A        x=x+h;   
, E' e+ W+ s+ u1 P9 H5 @1 q9 i    end
    T2=1/2*(T1+h*s);
: F0 h* {+ H6 W, c/ ?    n=n+1;
3 q! ]1 G: \7 g7 {+ l1 Uend
# M# V9 n  M8 gI=T2+1/3*(T2-T1);
: l5 I5 J2 a* b+ x# _) X% `- G! J( T6 K 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
a=a+eps;
h=b-a;
n=1;
: W- Z, X* Q/ uep=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- h9 z2 ]* m. X! |0 T    u=h/2;
  j) h2 J" B% M# x' N    x=a+u;
- U; }8 O  Q! E9 {; K    H=0;
+ O: O! ~" K0 v" B    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
' _! d' }0 X' R4 {    h=u;
    T1=T2;
    n=n+1;
- B/ P) X3 C% Gend
5 R, ]* Y8 X" Y+ M- j; [  M6 dI=T2+1/3*(T2-T1);