help-octave
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Help with some kickass demo


From: Nicolas Borboën
Subject: Re: Help with some kickass demo
Date: Tue, 14 Jul 2015 18:59:41 +0200
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:31.0) Gecko/20100101 Thunderbird/31.7.0

On 14. 07. 15 17:47, Doug Stewart wrote:
On Tue, Jul 14, 2015 at 11:40 AM, Doug Stewart <address@hidden
<mailto:address@hidden>> wrote:
    On Tue, Jul 14, 2015 at 6:01 AM, Nicolas Borboën
    <address@hidden <mailto:address@hidden>> wrote:
[...]
        I've found "Visualizing the Lorenz Strange Attractor with Octave" on
        http://www.ibm.com/developerworks/library/l-datavistools/
        (Listing 4), but it is done with Octave 2.1.50 and I'm running
        the last one (3.8.1 in a ubuntu-like docker). >

try this

1
  function y = lorenz( x, t )
y = [10 * (x(2) - x(1));
x(1) * (28 - x(3));
x(1) * x(2) - 8/3 * x(3)];
endfunction

  x = lsode("lorenz", [3;15;1], (0:0.01:25)');
  plot3(x(:,1),x(:,2),x(:,3))

Dear Doug,
        looks awesome, thanks a LOT !

Attached as "plot3(x(:,001),x(:,3),x(:,2))" screenshot.

Have a nice day/night,
Nicolas


--
Nicolas Borboën        Responsable informatique
EPFL   École polytechnique fédérale de Lausanne
IGM                 Institut de génie mécanique
CH-1015 Lausanne 15  -  Station 9  -  ME C3 397
fix. +41 21 693 54 55     mob. +41 76 615 25 80
http://igm.epfl.ch           http://www.epfl.ch

Attachment: lorenz.png
Description: PNG image


reply via email to

[Prev in Thread] Current Thread [Next in Thread]