[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
lorenz.png
Description: PNG image