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.
> -------------------------------------------------------------------------=
----------------------