数学建模社区-数学中国
标题:
新手跪求指导
[打印本页]
作者:
heugaodie
时间:
2012-12-4 22:51
标题:
新手跪求指导
r = 3.04042*10^(-6);
; {1 O% M; |( h$ ?5 y
p[x_] := x^5 - (2 + r) x^4 + (1 + 2 r) x^3 - (1 - r) x^2 +
- @) Q( V) T; y: l6 ~# k7 }0 X3 J
2 (1 - r) x + r - 1
- j, V+ D* R* B6 f4 Y R2 D* X" k3 M
q[x_] := p[x] - 2 r*x^2;
) c8 B0 C. m8 D' y$ t3 S8 H7 v
a0 = 0.5; b0 = 1.5; delta = 10^(-6); k0 = 100;
- v8 [1 N* E' V5 O- @! b
a = a0; b = b0;
6 Z) p) J+ L9 Z
Do[x = (a + b)/2;
% K2 U9 R- F- L% ^
If[p[x] == 0, Break[], If[N[p[x]*p[b]] < 0, a = x, b = x]];
8 a" |4 f' S7 n" G) l$ Z. B- p
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
- n5 s/ l2 u! C: r4 z7 Z }
x1 = x
* ?: g. ^: c7 }
Clear[x]
' I, u$ }$ q5 a8 J" J
Do[x = (a + b)/2;
* o9 Z2 M# x5 W2 P- R
If[q[x] == 0, Break[], If[N[q[x]*q[b]] < 0, a = x, b = x]];
# R. H1 a6 N& h; U/ V
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
: t6 ]5 Q' i& [3 D) X9 ~/ X# b8 t6 w
x2 = x
7 s0 W* Q) R' V! X, Y
Clear[x]
! M5 d1 ?6 N, X- C
+ X- `6 r% A5 i" t, r/ V
我想让p q 两函数的解都求出来 每次都是两个一样的解（是第一个的） 明明用了clear的 怎么做
( |- B6 k7 y. H9 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