贴一下代码吧,不知道会不会被度娘改。
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}
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}




