Lyapunov Exponent question about LCE package

Asked Mar 12 '22 at 10:21

Active Mar 12 '22 at 10:45

Viewed 149 times

This code gives result as given below on my computer

Needs["LCE`"]; (*or*) Needs["lce`"];
F[{x_, y_, z_}] := {16 (y - x), x (45.92 - z) - y, x y - 4 z};
LCEsC[F, {19, 20, 50}, 0.1, 800, 20, 0.02]
LCEsC[F, {19, 20, 50}, 0.1, 800, 20, 0.02, LCEsPlot -> True]

{{1.50085, 0.00327127, -22.4872}, 2.06689}

However, it should have been as shown below as mentioned in the tutorial.

How can I edit it? PLEASE GUIDE ME

edited Mar 12 '22 at 10:45

Roman

47,322
2
55
121

asked Mar 12 '22 at 10:21

user444

2,414
1
7
28

4

Try this Show[ LCEsC[F, {19, 20, 50}, 0.1, 800, 20, 0.02, LCEsPlot -> True], PlotRange -> All] – LouisB Mar 12 '22 at 11:19
Thank you so much – user444 Mar 12 '22 at 11:52
@LouisB or the author of the OP: out of curiosity, is this LCE package, the one written for the Sandri methods or another one? – Mar 12 '22 at 19:27
1

@kcr I used the Sandri package for Mathematica > 7 downloaded 2022-03-11 from www.msandri.it/soft.html to duplicate the results in the OP. – LouisB Mar 12 '22 at 21:31
@LouisB thanks a lot! – Mar 12 '22 at 21:32
1

FYI, I updated this Sandri's code here to use NDSolve instead of Euler method. It should be faster / more accurate. – Chris K Mar 12 '22 at 23:17

Lyapunov Exponent question about LCE package

0 Answers0