关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
n=1;
" s; P! c2 v$ j# O, v7 v* \: o5 geps=1e-3;
T1=0;
T2=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;
. ]( b0 ]# G1 q& n0 _    T1=T2;
    x=a+h;
% [3 o3 S9 P" w0 c    s=0;
    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
$ b7 C  u+ {' p4 J    T2=1/2*(T1+h*s);
; G$ g6 Q" r, i# X3 A    n=n+1;
" P( a/ j* A! yend
' I9 q) @- r( s, N' l$ r2 Y, TI=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
7 [' k. ~9 Z& @: q( {! Hfunction [I,n]=FHTX(f,a,b,ep)
a=a+eps;
h=b-a;
n=1;
ep=1e-6;
# t2 j* l* b7 n) c1 NT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
1 w& y$ }/ h  d, y, Htol=1;
9 o  o# P) t) v" |while tol>ep
8 J# {, ?2 k, Q2 t" i/ [+ v7 B3 N    u=h/2;
    x=a+u;
! S  G- k, q+ K    H=0;
1 \. {8 W' _8 s' O. X    while x<b
( {; }1 v) i1 n1 \/ @$ M# t2 I        H=H+subs(sym(f),findsym(sym(f)),x);
/ X8 n4 h. H$ o  G        x=x+h;   
* o4 B( }$ u8 R1 V* T1 S0 I+ C: d8 s    end
1 G: ]# o* e. T, A/ X- |    T2=1/2*(T1+h*H);
2 U1 P5 k/ |3 G: O    tol=abs(T2-T1);
. ~$ ~$ f5 l; B. y    h=u;
    T1=T2;
    n=n+1;
" H& ^$ h% P& Vend
I=T2+1/3*(T2-T1);
! T6 Q; w; `7 U