Lorenz Attractor

# -*- coding: utf-8 -*-
# увидеть результат можно на сайте atm-irbis.tk
import matplotlib.pyplot as plt

class LorenzSystem:
    def __init__(self):
        self.x = [3.051522]
        self.y = [1.582542]
        self.z = [15.62388]
    def generate(self):
        i = 0
        di = 0.0001
        while i < 100:
            x = self.x[-1]
            y = self.y[-1]
            z = self.z[-1]
            x1 = x + 5 * (-x + y) * di
            y1 = y + (15 * x - y - z * x) * di
            z1 = z + (-z + x * y) * di
            self.x.append(x1)
            self.y.append(y1)
            self.z.append(z1)
            i+=di
    def coords(self):
        return self.x,self.y,self.z


v = LorenzSystem()
v.generate()
# не удивляйтесь, что здесь три координаты, так как планировалась версия для 3D построения
x,y,z = v.coords()
plt.plot(x,y)
plt.show()