关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
( H9 U) b" p. X; Q/ ~2 {7 uh=b-a;
n=1;
eps=1e-3;
  H3 q3 M4 A9 u2 h+ g0 w% yT1=0;
3 v6 u+ v+ T0 FT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
# `, ^  {9 P! Dwhile abs(T2-T1)>eps
    h=h/2;
6 A" c( R# K1 Z- H6 R( H3 \    T1=T2;
    x=a+h;
+ s9 n- x' M: _4 ?6 i0 Y    s=0;
6 Q+ E  J, ]. B. B$ e7 b4 f    while x<b
! A0 e, [5 a" d) ^8 P  I* K        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
: ]/ l4 ~9 w7 _9 L3 ?2 U- F4 T    T2=1/2*(T1+h*s);
    n=n+1;
end
  M7 X* G, u# P( Z! h" k  [I=T2+1/3*(T2-T1);
: _5 p. b/ f0 x9 e+ C3 R 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
4 m' Q2 {% c/ N9 U' w5 u+ Ha=a+eps;
h=b-a;
" M9 C3 `6 V( X. fn=1;
ep=1e-6;
" b: d- D0 Q7 ]. {- \: {T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
while tol>ep
/ s* y* D2 \$ N- W4 n8 L    u=h/2;
& L5 j  E" |, w: [& T' Q' f0 E0 h/ H    x=a+u;
9 P* P. n; m8 N  l3 R! e4 l! n# \7 Q2 R    H=0;
    while x<b
        H=H+subs(sym(f),findsym(sym(f)),x);
7 b  X' I8 F! U8 w% L3 ^) R: C) Z7 A        x=x+h;   
    end
    T2=1/2*(T1+h*H);
* ?6 E, ^: X! f3 D3 R& C    tol=abs(T2-T1);
    h=u;
    T1=T2;
- u: O# I+ {  b# N7 {1 ]' b    n=n+1;
4 p) w( b3 I) X9 [9 v# Bend
( i& |1 @& Y7 y5 X2 pI=T2+1/3*(T2-T1);