Tuesday, April 12, 2016

Lorenz Attractor in Python with Matplotlib.

 Using three mapping.



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()



.