关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
8 Q% V6 X6 `3 |) w7 V) o3 L1 Th=b-a;
n=1;
7 S7 }: H" _% G1 N* c& A" Reps=1e-3;
  y! _6 v# G2 G- y, ~9 J2 r0 E# FT1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
2 T% K$ H6 I+ W0 m4 o! a# m. L" H$ e" ywhile abs(T2-T1)>eps
    h=h/2;
! e- E$ j. k6 ~! f8 _    T1=T2;
    x=a+h;
    s=0;
7 J- z+ w/ {& H, d2 D. U1 ?    while x<b
        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
. f9 O3 R8 ~  S" N% |$ I* o' W- _3 z    T2=1/2*(T1+h*s);
' d' U9 u; P' [3 z    n=n+1;
( x' L2 K5 N& i  O+ ^6 j: }6 v! Zend
# j3 K0 e" E! k- ]: oI=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
$ B. H# @8 k6 b, Y9 |! Lfunction [I,n]=FHTX(f,a,b,ep)
a=a+eps;
h=b-a;
0 [2 K% m% u2 K9 V" y  gn=1;
9 R6 ~$ N2 H, mep=1e-6;
# x* a/ n1 D; vT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
while tol>ep
1 a9 j/ S' l2 K+ V    u=h/2;
) j8 Z0 U1 V$ Y, F, r    x=a+u;
( v+ |3 K& s  l: G" V# G' U$ \. a) C    H=0;
. g3 K9 f2 q# }" y    while x<b
7 @' O( m! Q0 m) `- P% N/ _9 Q; d4 G        H=H+subs(sym(f),findsym(sym(f)),x);
8 V# v& A; Y+ I7 `% {% ^, g        x=x+h;   
1 `2 E  Y0 Z' T# y5 [  g    end
% X) W( k+ {- J    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
: K6 X" S( ]& I    h=u;
    T1=T2;
2 U, M" Q& K! y( y# |1 q    n=n+1;
end
I=T2+1/3*(T2-T1);
( d9 B' }# l7 P. q4 K7 a' A; j