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

```

• Prev by Date: Re: Fitting a Function with no closed form
• Next by Date: Re: axis label position
• Previous by thread: Problems in plotting Lyapunov Exponents:
• Next by thread: Re: Problems in plotting Lyapunov Exponents: