Lorenz Attractor in Blender

lorenz.py

# as explained by QuantitativeBytes
# https://youtu.be/iDHWyt1v4T8

import bpy

class Lorenz:
    def __init__(self, sceneRef, objName, color, initX, initY, initZ):
        self.X, self.Y, self.Z = initX, initY, initZ
        self.dt = 0.005
        self.a, self.b, self.c = 10, 28, 8/3
        self.color = color
        self.objName = objName
        self.sceneRef = sceneRef
        
    def Step(self):
        self.X = self.X + (self.dt * self.a * (self.Y - self.X))
        self.Y = self.Y + (self.dt * (self.X * (self.b - self.Z) - self.Y))
        self.Z = self.Z + (self.dt * (self.X * self.Y - self.c * self.Z))
        
    def Generate(self):
        
        numPoints = 10000
        
        self.curve = bpy.data.curves.new("LorenzCurve", type='CURVE')
        self.curve.dimensions = '3D'
        self.curve.bevel_depth = 0.1
        
        attractorPoly = self.curve.splines.new('POLY')
        attractorPoly.points.add(numPoints - 1)
        
        for i in range(0, numPoints):
                attractorPoly.points[i].co = (self.X, self.Y, self.Z, 1)
                self.Step()
                
        self.body = bpy.data.objects.new('curve', self.curve)
        
        self.body.name = self.objName
        self.sceneRef.collection.objects.link(self.body)
        
scene = bpy.context.scene

attractor1 = Lorenz(scene, "attractor1", (0.4, 0.1, 0.1, 1.0), 0.1,0.0,0.0)
attractor1.Generate()