新手跪求指导

heugaodie

r = 3.04042*10^(-6);
; E# d- L8 s" K1 w2 k0 gp[x_] := x^5 - (2 + r) x^4 + (1 + 2 r) x^3 - (1 - r) x^2 +
! l! ?+ V. a! j1 b+ X  2 (1 - r) x + r - 1
7 _: e6 ]! m7 bq[x_] := p[x] - 2 r*x^2;
a0 = 0.5; b0 = 1.5; delta = 10^(-6); k0 = 100;
1 i0 L5 E( K% T% ]2 n! h) [$ wa = a0; b = b0;
0 i" ^! L, H& \/ Y( mDo[x = (a + b)/2;
If[p[x] == 0, Break[], If[N[p[x]*p[b]] < 0, a = x, b = x]];
  O1 n* T8 \1 |5 b7 `3 _7 G If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
x1 = x
Clear[x]
! Q* }5 z; f- y3 \. t- zDo[x = (a + b)/2;
+ _$ E- w' G2 S$ d0 Q* _ If[q[x] == 0, Break[], If[N[q[x]*q[b]] < 0, a = x, b = x]];
0 t+ T+ z. b2 ^; {" o5 S8 A If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
x2 = x
Clear[x]
: y2 F; K$ V% q/ c  o6 W
/ h9 x2 C& y( C/ N, B7 q我想让p q 两函数的解都求出来 每次都是两个一样的解（是第一个的）  明明用了clear的 怎么做
3 p$ m# m* d7 b2 |4 `# k" ^

qudazi

你试一下将Clear[x]换成x=.看看

花如梦

别只用x,换个变量试试

小小灰

看不懂，看不懂