数学建模社区-数学中国
标题:
新手跪求指导
[打印本页]
作者:
heugaodie
时间:
2012-12-4 22:51
标题:
新手跪求指导
r = 3.04042*10^(-6);
7 s# @- b6 i) f3 k* T4 v% v6 t
p[x_] := x^5 - (2 + r) x^4 + (1 + 2 r) x^3 - (1 - r) x^2 +
1 i5 I" q* j& O) K" ?
2 (1 - r) x + r - 1
/ h% l) R8 l( Z$ |
q[x_] := p[x] - 2 r*x^2;
% y7 a" B9 A" P9 u0 m% X+ m* R
a0 = 0.5; b0 = 1.5; delta = 10^(-6); k0 = 100;
, ]' b0 L% n0 M: {2 C! f- Q d
a = a0; b = b0;
" O1 {& w! K+ c) v4 U( J
Do[x = (a + b)/2;
8 m/ J. a- v6 c% M5 R8 X( M
If[p[x] == 0, Break[], If[N[p[x]*p[b]] < 0, a = x, b = x]];
; a* C' ?; U* P4 G* N
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
' w+ ~* Q: }; s) |
x1 = x
. O5 }) M- T3 F1 l$ n
Clear[x]
% G" k, {8 ]! o& B. m
Do[x = (a + b)/2;
9 {5 U( [5 f- ]% Z8 X( Q8 D
If[q[x] == 0, Break[], If[N[q[x]*q[b]] < 0, a = x, b = x]];
. h# t3 Y& k' j$ K
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
3 L9 b) p2 I& L. |# l! ?/ }3 h
x2 = x
0 z6 M( _. w! W; O( j
Clear[x]
, d2 z9 i _( K* c$ g. m0 R9 B5 H
, E1 J7 y6 \; e# N
我想让p q 两函数的解都求出来 每次都是两个一样的解（是第一个的） 明明用了clear的 怎么做
% S, H8 `& i+ D
作者:
qudazi
时间:
2012-12-6 09:08
你试一下将Clear[x]换成x=.看看
作者:
花如梦
时间:
2014-4-12 20:35
别只用x,换个变量试试
作者:
小小灰
时间:
2014-4-28 16:57
看不懂，看不懂
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5