- 在线时间
- 0 小时
- 最后登录
- 2012-12-5
- 注册时间
- 2012-12-1
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 5 点
- 威望
- 0 点
- 阅读权限
- 10
- 积分
- 2
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 1
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 2
升级   40% 该用户从未签到 - 自我介绍
- 大神
 |
r = 3.04042*10^(-6);* f. K8 _8 E- z9 n
p[x_] := x^5 - (2 + r) x^4 + (1 + 2 r) x^3 - (1 - r) x^2 +
9 v6 O$ J4 [* {, _1 ` R 2 (1 - r) x + r - 11 N' v4 }6 z. T: R* S1 p6 g
q[x_] := p[x] - 2 r*x^2;, V' t: e( A7 T6 ?1 n& I5 F) Y
a0 = 0.5; b0 = 1.5; delta = 10^(-6); k0 = 100;6 b, T4 D0 L( h6 b2 k/ C
a = a0; b = b0;0 x: u. E4 h) @# ?3 O5 _/ u5 r
Do[x = (a + b)/2;
$ [. g! p- L: F1 k+ Q If[p[x] == 0, Break[], If[N[p[x]*p[b]] < 0, a = x, b = x]];% j5 K* `' h7 n8 p, b; h% W
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]) w* h3 |3 R1 B
x1 = x
7 L9 I6 l( K# a3 }' J. m4 |4 IClear[x]) @' V" M1 C. Y/ N, ?& O
Do[x = (a + b)/2;
$ T" x: ]/ a: S( H If[q[x] == 0, Break[], If[N[q[x]*q[b]] < 0, a = x, b = x]];: {/ s" a7 W9 L9 v" |
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
: m0 X, x" t; y# |* S5 xx2 = x+ ^" B/ y" S9 y: n6 I' J
Clear[x]
* j+ Q. q' J6 V- o2 |3 d' L2 s, O# ]7 V- h8 j5 p* @* \ t4 I
我想让p q 两函数的解都求出来 每次都是两个一样的解(是第一个的) 明明用了clear的 怎么做" Z+ i3 k" O+ q( }
|
zan
|