Problems in plotting Lyapunov Exponents:
- To: mathgroup at smc.vnet.net
- Subject: [mg98649] Problems in plotting Lyapunov Exponents:
- From: CreativeSolutionsNZ at gmail.com
- Date: Wed, 15 Apr 2009 05:00:17 -0400 (EDT)
Hi
I have been trying to debug the following piece of code
to plot Lyapunov Exponents, but getting several errors.
-----------------------------------------------------------------------------------------------
In[1]:=LypHS=Compile[{l, xinit, {n, _Complex}, {ndrop, _Complex}},
xlist=Drop[
NestList[-30000l+(1-l)#+l*(2858.16/(#-500)^0.82)&, xinit,
n],
ndrop+1];
Apply[Plus, Log[Abs[(1-l)-l(2343.71/(xlist-500)^1.82)]]]/
Length[xlist]];
Plot[LypHS[l, 10, 5000, 100], {l, 0.0, 1.0}]
-----------------------------------------------------------------------------------------------
Can anyone try to run it and help me to fix these errors ?
Or is any better code is possible for this problem ?
Best regards.
-----------------------------------------------------------------------------------------------