Re: Problems in plotting Lyapunov Exponents:

*To*: mathgroup at smc.vnet.net*Subject*: [mg98700] Re: Problems in plotting Lyapunov Exponents:*From*: "Sjoerd C. de Vries" <sjoerd.c.devries at gmail.com>*Date*: Thu, 16 Apr 2009 04:17:36 -0400 (EDT)*References*: <gs47mm$7q5$1@smc.vnet.net>

I don't think this is going to work this way. n and ndrop shouldn't be of type Complex, they're Integers. Compile needs to know of xlist as well (3rd argument of Compile), but there is a problem, since it seems that Lists needs to have fixed size (I might be wrong here). Cheers -- Sjoerd On Apr 15, 11:00 am, CreativeSolution... at gmail.com wrote: > 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)&, xi= nit, > 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. > -------------------------------------------------------------------------= ----------------------