数学建模社区-数学中国

标题: 新手跪求指导 [打印本页]

---

作者: heugaodie    时间: 2012-12-4 22:51
标题: 新手跪求指导
r = 3.04042*10^(-6);
p[x_] := x^5 - (2 + r) x^4 + (1 + 2 r) x^3 - (1 - r) x^2 +
$ E- y" [% O' \* I( {" b# P: H2 `0 d  2 (1 - r) x + r - 1
q[x_] := p[x] - 2 r*x^2;
$ W8 z- j: r) d1 z8 y! X7 xa0 = 0.5; b0 = 1.5; delta = 10^(-6); k0 = 100;
a = a0; b = b0;
! [0 S& z3 L0 I4 DDo[x = (a + b)/2;
If[p[x] == 0, Break[], If[N[p[x]*p[b]] < 0, a = x, b = x]];
) j/ q* i+ c1 r+ R$ h! U If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
3 e# h3 W/ W( t0 nx1 = x
Clear[x]
! d- l* M. i$ k- M' D# J# z  }6 XDo[x = (a + b)/2;
: `' N0 c" O3 G8 u( i; a9 w: F If[q[x] == 0, Break[], If[N[q[x]*q[b]] < 0, a = x, b = x]];
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
# c) E5 P4 ~. J, b4 Zx2 = x
Clear[x]
) Z( T( F+ _' j7 d5 E) X
我想让p q 两函数的解都求出来 每次都是两个一样的解（是第一个的）  明明用了clear的 怎么做
: \+ f* j4 K/ s' P

---

作者: 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