数学建模社区-数学中国
标题:
新手跪求指导
[打印本页]
作者:
heugaodie
时间:
2012-12-4 22:51
标题:
新手跪求指导
r = 3.04042*10^(-6);
}+ A9 Y' L* c! G
p[x_] := x^5 - (2 + r) x^4 + (1 + 2 r) x^3 - (1 - r) x^2 +
( D* O+ \ E8 L& c$ ?* x/ l
2 (1 - r) x + r - 1
: o$ ]2 ] i9 u
q[x_] := p[x] - 2 r*x^2;
) Z7 v8 R& F6 v
a0 = 0.5; b0 = 1.5; delta = 10^(-6); k0 = 100;
+ V7 y6 U5 ]3 _8 w, Q1 X% h
a = a0; b = b0;
: c- Y, c3 U* Y( N7 @' U6 _
Do[x = (a + b)/2;
1 y7 |" W) E: z; Q
If[p[x] == 0, Break[], If[N[p[x]*p[b]] < 0, a = x, b = x]];
- z) m. Z! K5 r( ~6 u1 {8 r
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
/ L$ H$ N0 Y9 M+ M2 Y
x1 = x
2 d: @* f8 p ], k _; j' ^9 |3 }
Clear[x]
u: H2 g: F. |! M. f
Do[x = (a + b)/2;
! h: h) _) @& ~+ M5 [$ A# E
If[q[x] == 0, Break[], If[N[q[x]*q[b]] < 0, a = x, b = x]];
, X8 F' t- t( F; C9 L- W
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
! D: q# I6 F, T$ q: X) q# {2 m9 B
x2 = x
. y/ a' h* \( I
Clear[x]
4 G6 C- y; l, Q* T5 i! G1 o! R
. C) `2 t# B- A$ W; u9 n
我想让p q 两函数的解都求出来 每次都是两个一样的解(是第一个的) 明明用了clear的 怎么做
& U! X0 T0 d8 J4 w
作者:
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