关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
h=b-a;
n=1;
9 y+ }- B; `" d) E( [# Deps=1e-3;
T1=0;
& M- B, m. |, Z) ~3 h3 T6 ~$ k! ET2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
4 ]* h/ S( ]0 O. `% s8 A* Swhile abs(T2-T1)>eps
& k1 e& O4 }( a$ b, i. l7 [    h=h/2;
/ H, M5 l" n- `3 J/ a    T1=T2;
    x=a+h;
! l5 t. K' i- ?/ U- t    s=0;
+ U) |+ J7 k# J  @; Y    while x<b
! R. \$ I4 [$ r5 f        s=s+subs(sym(f),findsym(sym(f)),x);
. @. ]1 |- Z. l1 X4 o' t7 i        x=x+h;   
    end
    T2=1/2*(T1+h*s);
    n=n+1;
0 w% f1 R! X/ s0 s* p" \end
. H; Y9 W1 k  x; ~5 E4 bI=T2+1/3*(T2-T1);
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
, t. D) d8 `" G( z( qa=a+eps;
5 w4 e2 E& d; @0 I6 S& ph=b-a;
1 P' p1 H# c4 z; Cn=1;
; F2 k4 c( @7 E; ^3 w7 nep=1e-6;
( n  K$ F: T4 [T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
5 V  W! }+ {# ?) |3 awhile tol>ep
& B* I8 b/ V6 p. Y' d, k    u=h/2;
) O9 w5 ?# b3 b6 [  n# W* I4 Y    x=a+u;
    H=0;
, _1 b' ?: a0 P2 h, p    while x<b
* W+ p4 N* E" @  D) T        H=H+subs(sym(f),findsym(sym(f)),x);
4 V0 E7 }. n2 n: i) ?. Y5 e        x=x+h;   
3 M5 Q3 Z3 a1 ~    end
$ {1 y4 e/ k$ J; ~  ]    T2=1/2*(T1+h*H);
    tol=abs(T2-T1);
: d# x. w! J; T- T. @) [    h=u;
    T1=T2;
    n=n+1;
% f% e' s7 U" j. ]2 ~end
  H- V+ K# ~7 A7 T, V& OI=T2+1/3*(T2-T1);