import numpy as np #untuk operasi array
import matplotlib.pyplot as plt #untuk gambar grafik
import matplotlib.animation as animation #untuk menggerakkan grafik
fig, (ax, ay, az) = plt.subplots(3,sharex=True)
x = 1.
y = 1.
z = 1.
#plt.ylim(-43,43)
#plt.xlim(-43,43)
#membuat garis/kurva dengan sumbu-x adalah x, sumbu-y adalah y
line, = ax.plot(x, y, 'o' )
line, = ay.plot(x, z, 'o' )
line, = az.plot(y, z, 'o' )
def animate(i):
global x,y,z
dt = 1./64.
s = 10.
b = 8./3.
r = 28
xdot = s * (y-x)
ydot = x*r -x*z -y
zdot = x*y -b*z
x = x+xdot*dt
y = y+ydot*dt
z = z+zdot*dt
plt.figure(1)
line, = ax.plot(x, y, 'o' )
line, = ay.plot(x, z, 'o' )
line, = az.plot(y, z, 'o' )
return line,
ani = animation.FuncAnimation(fig, animate, frames=2000, interval=100, blit=False)
#ani.save('Lorenz.mp4',bitrate=1024)
plt.show()
Pages
▼
Tuesday, April 12, 2016
Lorenz Attractor in Python with Matplotlib.
Using three mapping.
.

No comments:
Post a Comment