数学建模社区-数学中国
标题:
关于数值积分的程序
[打印本页]
作者:
753085848
时间:
2010-11-22 21:08
标题:
关于数值积分的程序
function [I,n]=TX(f,a,b,eps)
) [6 a$ B5 i1 {7 o' j# k# O. U- P
h=b-a;
5 O4 x9 t2 S3 `" D. y- Z
n=1;
4 @% ?& B+ T* l3 a4 B3 F7 d7 s/ Q
eps=1e-3;
5 ^9 C+ W7 g, I! h4 h
T1=0;
5 q. R7 t8 m5 [! W8 ]
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
- D( n, N6 S. _% A( z
while abs(T2-T1)>eps
6 ~/ Z) K% }( Z) Q4 W( w! S
h=h/2;
5 C% D! X9 I5 D
T1=T2;
1 z+ n( m3 D( S4 W; X! w
x=a+h;
7 i5 b5 L- b ^: A) |# o' C+ s
s=0;
% L. p) a' s3 H
while x<b
6 e" t2 R; ~1 p9 L2 Q/ |
s=s+subs(sym(f),findsym(sym(f)),x);
8 `! ^5 `+ x o5 k
x=x+h;
+ m5 i+ _: O, b$ a7 F) g
end
5 Z: l4 S6 A7 M' Y9 j, ]
T2=1/2*(T1+h*s);
. g: D& S; p' I. d% C* W
n=n+1;
{8 W0 J( }/ B4 i( J. F& O
end
$ M! p8 c' U5 O
I=T2+1/3*(T2-T1);
, P5 X( L8 Z0 U8 z2 M& R5 J( y
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！
作者:
753085848
时间:
2010-11-23 17:43
自己修改过的，是正确的！
) h! ?/ \3 T6 i7 F9 K& k$ E% O5 j
function [I,n]=FHTX(f,a,b,ep)
" F' O- p; [' ]+ O
a=a+eps;
) y4 H. n+ C3 g% F H2 {4 C) p2 p
h=b-a;
: U* T1 @, {5 e- `! S8 ~1 R9 U
n=1;
0 h- p( R7 S/ N- w* ~
ep=1e-6;
' H7 P/ X2 U; F7 {8 b- Z
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
6 A7 F" \* i- N/ C( i `9 ~8 l
tol=1;
% I b' T/ V/ Z* Q; D
while tol>ep
. X' ~( D5 _4 m
u=h/2;
6 _( D! W1 s" M8 u' j
x=a+u;
8 u6 N" l, Y0 \- m: c7 ?# o
H=0;
. P( I, i% ]* [" v* ]
while x<b
4 j2 f7 U! K. B- `7 h6 m9 Q1 C+ o
H=H+subs(sym(f),findsym(sym(f)),x);
, |! @2 @8 T. _( ]/ a
x=x+h;
+ b2 \4 W9 u7 x! x# v y, W: a
end
+ N6 ?2 Z- ~3 h
T2=1/2*(T1+h*H);
* x0 j8 N/ y* |( k% t" W3 g
tol=abs(T2-T1);
$ S2 }; c8 X# X( Q6 O( H
h=u;
& m3 Y; }6 i! Z5 k) ?
T1=T2;
5 A- p1 z7 P8 i0 a1 B" z
n=n+1;
0 `) h! g9 N- O* p2 S
end
8 W& o' | x# U8 b9 m
I=T2+1/3*(T2-T1);
3 Z- m3 I# ?0 `
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！