| Author |
Comment/Response |
Rimas
|
 05/19/12 1:25pm
Hello!
I need to calculate Lyapunov exponents and plot their dependence on control parameter r for a map:
x_n+1 = r*x_n-x_n^3.
Here is a source code
strike1 =
Compile[{{r, _Real}, {n, _Integer}},
1/n*Log[Abs[r - 3*#^2]] & /@ NestList[r*# - #^3 &, .2, n]];
ListPlot[{#, strike1[#, 600]} & /@ Range[0.1, 2],
PlotRange -> {{0, 3}, {-5, 2}}, AxesOrigin -> {0, 0}]
but something is wrong. No error messages appear, but no plot also - only axes.
Please, help. What's wrong with it?
URL: , |
|
