新手跪求指导

heugaodie

r = 3.04042*10^(-6);
p[x_] := x^5 - (2 + r) x^4 + (1 + 2 r) x^3 - (1 - r) x^2 +
1 v7 F' H# n" @+ N/ p  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;
  I& S3 q" S4 u( o: z7 _. s# }5 |Do[x = (a + b)/2;
If[p[x] == 0, Break[], If[N[p[x]*p[b]] < 0, a = x, b = x]];
( u& R+ m3 A% d+ v: J If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
x1 = x
. ?7 s0 D8 f+ \/ v& }' h8 ?Clear[x]
Do[x = (a + b)/2;
6 [  u3 V8 b9 X3 S) N- X5 z) Z: e 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}]
$ v# }5 z2 o" Z/ \, r% W% ], [% |# c7 Hx2 = x
) f/ A  u, P, r  v; w, HClear[x]
2 d$ v' w; L. }- Z' h7 S
我想让p q 两函数的解都求出来 每次都是两个一样的解（是第一个的）  明明用了clear的 怎么做
- F) f1 ]. q  u0 ?7 N; y

qudazi

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

花如梦

别只用x,换个变量试试

小小灰

看不懂，看不懂