关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
( T! `( @2 }' S$ [  Uh=b-a;
2 x' u! Y1 m1 s7 U/ F8 \7 In=1;
eps=1e-3;
T1=0;
( T% v4 K! L+ L6 c0 aT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
    h=h/2;
& h7 f3 V! R4 v; N    T1=T2;
    x=a+h;
    s=0;
+ m' p. }" I2 ?8 j    while x<b
2 D# {3 B/ ]' u        s=s+subs(sym(f),findsym(sym(f)),x);
1 S$ C0 L: V* R        x=x+h;   
    end
; {/ W) q7 A0 ]: T8 @: D; i7 l    T2=1/2*(T1+h*s);
    n=n+1;
; L, @/ g. X( D$ qend
I=T2+1/3*(T2-T1);
; ^( a9 W5 z4 L) b) G  R6 J 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
, r  t2 s0 S& @% G3 X# Xa=a+eps;
h=b-a;
n=1;
. o2 y; u7 q6 ]9 f1 M7 Kep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
" V1 b# [* h' u& G/ Xwhile tol>ep
% P9 o5 I# i- `) E8 B    u=h/2;
    x=a+u;
) Z" F  L& h" h    H=0;
0 x/ ~; ^3 e2 `) o7 F& i/ v( B1 A    while x<b
3 K' w# G# ]6 X$ `2 e        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
0 r  F% I# F$ j9 Y, q* d6 t    end
    T2=1/2*(T1+h*H);
2 _4 E, o( h+ e' R, D    tol=abs(T2-T1);
    h=u;
* q. v4 k" J2 t3 y% }8 U    T1=T2;
9 U( ^4 |2 R: ]9 \! b! C2 f( ]    n=n+1;
% ?5 r! t2 A) h4 uend
I=T2+1/3*(T2-T1);
9 U$ D5 t  y$ G7 k5 M+ V' `0 q1 X8 _