数学建模社区-数学中国
标题:
新手跪求指导
[打印本页]
作者:
heugaodie
时间:
2012-12-4 22:51
标题:
新手跪求指导
r = 3.04042*10^(-6);
. A( T, I3 D* L% ]4 s* e1 r
p[x_] := x^5 - (2 + r) x^4 + (1 + 2 r) x^3 - (1 - r) x^2 +
' y, }& N, [9 [5 p3 P* z4 B; m
2 (1 - r) x + r - 1
/ ~2 Q; E5 W' p
q[x_] := p[x] - 2 r*x^2;
# B' E/ x' B S m
a0 = 0.5; b0 = 1.5; delta = 10^(-6); k0 = 100;
( m9 }" s B9 r8 _
a = a0; b = b0;
5 s5 A l- C8 `- V
Do[x = (a + b)/2;
9 m j$ _# M' y# A0 G. Q
If[p[x] == 0, Break[], If[N[p[x]*p[b]] < 0, a = x, b = x]];
9 P# g% D) k! k0 E, @9 _4 |
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
7 K: t. z" s: f8 f, `# i# I& ]
x1 = x
0 l$ W2 K+ F) h& Y9 C2 {
Clear[x]
; V5 f: W4 Y. J3 {& B+ e7 o! t8 K
Do[x = (a + b)/2;
$ X9 Z) V# D y& s
If[q[x] == 0, Break[], If[N[q[x]*q[b]] < 0, a = x, b = x]];
: |! n' ?6 V6 T9 {
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
) X t8 D9 Y3 v
x2 = x
& T- q) V& ~* C0 |3 y
Clear[x]
6 f8 N+ K s8 g' [0 g- |
3 V5 o* K" I# D ^ ^# P
我想让p q 两函数的解都求出来 每次都是两个一样的解（是第一个的） 明明用了clear的 怎么做
" ]0 P4 b. `+ @. Y! 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