1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
# This example displays a beautiful strange attractor: the Lorenz
# attractor.
# You may run this example by doing:
#
# ode < lorenz.ode | graph -T X -C
#
# or alternatively, to get a real-time plot,
#
# ode < lorenz.ode | graph -T X -C -x -10 10 -y -10 10
#
# You may also produce and print a Postscript version by doing
#
# ode < lorenz.ode | graph -T ps -C -x -10 10 -y -10 10 -W 0 | lpr
#
# The `-W 0' sets the line width for the Postscript plot to
# be zero. That means that the thinnest line possible will be used.
# The Lorenz model, a third order system.
# Interesting cases are r = 26, 2.5<t<30, x = z = 0, y = 1
# and r = 17, 1<t<50, x = z = 0, y = 1.
x' = -3*(x-y)
y' = -x*z+r*x-y
z' = x*y-z
r = 26
x = 0
y = 1
z = 0
print x, y
step 0, 200
|
