关于数值积分的程序

753085848

function [I,n]=TX(f,a,b,eps)
$ y" Y5 J+ b0 t/ z8 ^: bh=b-a;
n=1;
eps=1e-3;
T1=0;
/ i' a3 G6 p+ lT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
/ Z+ h% {- s: i" L; d  o- _while abs(T2-T1)>eps
- A' \6 O0 J( q, ^3 \    h=h/2;
    T1=T2;
    x=a+h;
6 Z) h* r% Z6 h$ l    s=0;
+ _$ Y4 U& X: k/ G; ^- N: K    while x<b
! y& q. Z8 i5 @9 B5 u        s=s+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
    end
; D' I, n0 r2 Z% y6 Z7 t    T2=1/2*(T1+h*s);
7 L8 w& Z7 c, {# S    n=n+1;
end
; E! C) i$ C5 X/ iI=T2+1/3*(T2-T1);
: V: F: X4 I. a 这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

753085848

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
a=a+eps;
( j5 I2 T! |8 U$ c" _h=b-a;
. V4 E6 g$ V7 H# `1 s! d) An=1;
; _7 s" w" c$ r( o- c# |' i2 qep=1e-6;
$ R9 G8 a$ Q5 q* @- dT1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
while tol>ep
    u=h/2;
    x=a+u;
    H=0;
4 M& F1 x: ?3 q: U* K4 X3 v    while x<b
* K% [0 S) S# @! O+ h        H=H+subs(sym(f),findsym(sym(f)),x);
        x=x+h;   
* `' L5 |, |, u' t# U7 v+ Q    end
. y/ J! A$ z$ F$ K2 [    T2=1/2*(T1+h*H);
$ ]" P9 ]# d7 G    tol=abs(T2-T1);
    h=u;
    T1=T2;
* m( P5 @& ~2 f    n=n+1;
end
9 X/ N( v- d8 g; d8 ]I=T2+1/3*(T2-T1);