请问我为何无法重现这个程序?
mathematica吧
全部回复
仅看楼主
level 3
在英文wiki上看到这张图,想照猫画虎自己写一个却发现一运行就耗尽我E5机子的16G内存然后卡死。
但是看图上好像很快就能跑出结果,请问是我哪里写错了么?附上原图和我写的版本。谢谢了。度盘1bnpNfFx
2014年08月27日 17点08分 1
level 3
贴一下代码吧,不知道会不会被度娘改。
bifurcate[f_, a0_, k0_, k_, while_: (True &)] := NestWhileList[f, a0, while, 2, k0+k-1] // Drop[#, k0] &
logistic = {r, y} \[RightTeeArrow] r y (1-y);
Row@{"Initial value: ", a0 = RandomReal[{0.1, 0.9}]}
Row@{"Points per r: ", density = 10^2}
Row@{"Initial k: ", k0 = 10^4}
Row@{"Time taken: ", Timing[
plotData=ParallelTable[
{ConstantArray[r, density], bifurcate[logistic[r, #] &, a0, k0, density]}\[Transpose],
{r, 3.5, 4, 0.0001}] // Flatten[#, 1] &;][[1]], " s"}
Row@{"Data length: ", plotData // Length}
2014年08月27日 17点08分 2
1