Mathematica 9 is now available
Student Support Forum
-----
Student Support Forum: 'Lyapunov exponents' topicStudent Support Forum > General > Archives > "Lyapunov exponents"

Next Comment >Help | Reply To Topic
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: ,

Subject (listing for 'Lyapunov exponents')
Author Date Posted
Lyapunov exponents Rimas 05/19/12 1:25pm
Re: Lyapunov exponents jf 05/20/12 6:52pm
Re: Re: Lyapunov exponents Rimas 05/21/12 00:32am
Next Comment >Help | Reply To Topic