Butterfly effect : Lorens equations
Page 1 sur 1
Butterfly effect : Lorens equations
- Code:
rem ============================================================================
rem Butterfly effect : Lorens equations
rem Adaptation d’un code en Yabasic
rem Papydall Juin 2017
rem ============================================================================
Init()
Boucle()
end
rem ============================================================================
SUB Init()
dim xa(2), ya(2), za(2)
dim c1, dt, a, b, c, i, x, y, z, x1, y1, z1, xd, yd
picture 10 : width 10,width(0)-50 : height 10,height(0)-60 : top 10,10: left 10,20
2d_target_is 10 : color 10,0,0,0 : caption 0,"Butterfly effect .... <ESC> pour arrêter ..."
xa(0) = 3.051522 : xa(1) = 3.051522
ya(0) = 1.592542 : ya(1) = 1.582542
za(0) = 15.62388 : za(1) = 15.62388
c1 = 0.292893 : dt = 0.02 : a = 5 : b = 15 : c = 1
END_SUB
rem ============================================================================
SUB Boucle()
repeat
for i = 0 to 1
x = xa(i) : y = ya(i) : z = za(i)
x1 = x-a*x*dt+a*y*dt
y1 = y+b*x*dt-y*dt-z*x*dt
z1 = z-c*z*dt+x*y*dt
x = x1 : y = y1 : z = z1
xd = y-x*c1 : yd = z+x*c1
if i = 1
draw(xd,yd,128)
else
draw(xd,yd,255)
end_if
xa(i) = x : ya(i) = y : za(i) = z
next i
until scancode = 27
END_SUB
rem ============================================================================
' Draw a shiny ball
SUB Draw(xd,yd,cl)
dim_local i1,j1,c,c1,c2,d1,d2,i,k
c1 = 19.3 : c2 = width(10)/2
d1 = 11.0 : d2 = 392.0
k = 10
i1 = c1*xd+c2
j1 = 0-d1*yd+d2
for i = 127 to 255 step 16
c = 0.09*(10-k)
2d_pen_color cl,i,i : 2d_fill_color cl,i,i
2d_circle i1+c,j1+c,k
k = k - 1
next i
END_SUB
rem ============================================================================
Page 1 sur 1
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum
Lun 1 Jan - 0:25 par Papydall-Admin
» A ceux qui célèbre Noël, bonnes fêtes
Dim 24 Déc - 10:49 par Papydall-Admin
» Joyeux Noël et Bonne Année
Ven 8 Déc - 1:34 par Papydall-Admin
» Planets of the Solar System : Tilts and Spins
Lun 20 Mar - 15:43 par Papydall-Admin
» Bonne Année 2023
Sam 31 Déc - 1:39 par Papydall-Admin
» Fractals - Mandelbrot
Ven 21 Aoû - 22:51 par Papydall-Admin
» Convertisseur Décimal ---> Binaire, Octal, Hexadécimal, ...
Mer 21 Nov - 1:08 par Papydall-Admin
» Balises {USER...}
Lun 19 Nov - 22:12 par Papydall-Admin
» Useful Dog
Ven 6 Avr - 14:25 par Papydall-Admin