#/***********************************************************
#    lorenz.rb -- Lorenz (ローレンツ) アトラクタ
#***********************************************************/
require "plotter.rb"

A = 10.0
B = 28.0
C = 8.0 / 3.0
D = 0.01

gr_on
x = y = z = 1
for k in 0...3000
    dx = A * (y - x)
    dy = x * (B - z) - y
    dz = x * y - C * z
    x += D * dx;  y += D * dy;  z += D * dz
    if (k > 100); draw(x, z);  else; move(x, z); end
end

exit 0