新手跪求指导

heugaodie

r = 3.04042*10^(-6);
$ l# I  [8 y% u% b9 }$ B" p' w  q# @p[x_] := x^5 - (2 + r) x^4 + (1 + 2 r) x^3 - (1 - r) x^2 +
  2 (1 - r) x + r - 1
q[x_] := p[x] - 2 r*x^2;
a0 = 0.5; b0 = 1.5; delta = 10^(-6); k0 = 100;
a = a0; b = b0;
Do[x = (a + b)/2;
: G7 j7 e  T. U# u' G. C9 s5 { If[p[x] == 0, Break[], If[N[p[x]*p[b]] < 0, a = x, b = x]];
) T7 f" L4 Z$ A1 G7 Q2 t# t If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
  K' D# C% I5 t+ B2 b. Vx1 = x
Clear[x]
. G# H7 s5 d6 w' kDo[x = (a + b)/2;
0 U) u* c# w7 x/ `; w( z If[q[x] == 0, Break[], If[N[q[x]*q[b]] < 0, a = x, b = x]];
+ ?8 r' G7 L2 M  k: D/ L& Y( A0 S If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
# w5 k9 Y) y0 v, ]: Gx2 = x
3 m( g, Y" K! i' ?4 ]Clear[x]
0 ]" Q2 W& c6 y/ ~5 W8 }
/ }: A7 C3 i2 H$ ^4 u我想让p q 两函数的解都求出来 每次都是两个一样的解（是第一个的）  明明用了clear的 怎么做

qudazi

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

花如梦

别只用x,换个变量试试

小小灰

看不懂，看不懂