数学建模社区-数学中国
标题:
新手跪求指导
[打印本页]
作者:
heugaodie
时间:
2012-12-4 22:51
标题:
新手跪求指导
r = 3.04042*10^(-6);
: l! Q4 b" l0 K; d
p[x_] := x^5 - (2 + r) x^4 + (1 + 2 r) x^3 - (1 - r) x^2 +
Q1 k0 q$ G5 p, f6 G3 F2 B
2 (1 - r) x + r - 1
( r! p9 H% J) Y& b: |% |
q[x_] := p[x] - 2 r*x^2;
8 N$ F2 ~2 Y) z& t
a0 = 0.5; b0 = 1.5; delta = 10^(-6); k0 = 100;
$ P! ~4 V/ }$ y+ U+ @1 H7 g
a = a0; b = b0;
" `; l+ g0 j, L8 D
Do[x = (a + b)/2;
4 Z% t# z/ H1 J; R
If[p[x] == 0, Break[], If[N[p[x]*p[b]] < 0, a = x, b = x]];
5 X; [$ [1 K+ L- H/ `/ W0 ?
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
( J: V2 _' a0 ?6 s8 `7 V
x1 = x
8 R8 S+ d* `6 z/ ]( D
Clear[x]
0 P6 B) A1 [4 p5 i& s1 _
Do[x = (a + b)/2;
* d2 l5 p! |8 F' |; P8 l; s
If[q[x] == 0, Break[], If[N[q[x]*q[b]] < 0, a = x, b = x]];
0 J: w) d& m9 D# o- h) W8 u' s; T
If[N[Abs[b - a]] < delta, Break[], If[k == k0, Print[失败]]], {k, k0}]
& L: L; B+ e8 ?9 n9 P8 J" J
x2 = x
0 E5 `" n' g/ _! A6 v
Clear[x]
- ~5 D) I0 P7 V/ @
5 U3 @; t3 C4 _/ b: C
我想让p q 两函数的解都求出来 每次都是两个一样的解(是第一个的) 明明用了clear的 怎么做
_. O# v0 X a3 w4 @
作者:
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