MathGroup Archive 2009

[Date Index] [Thread Index] [Author Index]

Search the Archive

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: