"""
Moore Spiegel Attraktor
(c) Juergen Meier (www.3d-meier.de)
25.12.2021
"""

import c4d
import math

def CreateSplineObject():
  # Variablen und Konstanten
    Name = "Moore Spiegel Attractor"
    a = 100            # Konstante Moore Spiegel Attraktor
    b = 26             # Konstante Moore Spiegel Attraktor
    x0 = 1.0           # Startwerte
    y0 = 1.0           # Startwerte
    z0 = 1.0           # Startwerte
    Delta = 0.001      # Delta Wert fuer die Berechnung
    N = 20000          # Anzahl Punkte
    Faktor = 10        # Skalierungsfaktor
    x = x0             # Startwerte zuweisen
    y = y0             # Startwerte zuweisen
    z = z0             # Startwerte zuweisen

    obj = c4d.BaseObject(c4d.Ospline)       # Splineobjekt erzeugen
    obj.ResizeObject(N)                     # Groesse des Spline festlegen
    obj[c4d.SPLINEOBJECT_CLOSED] = False
    #obj[c4d.SPLINEOBJECT_TYPE] = c4d.SPLINETYPE_LINEAR
    obj[c4d.SPLINEOBJECT_TYPE] = c4d.SPLINETYPE_BSPLINE
    #obj[c4d.SPLINEOBJECT_TYPE] = c4d.SPLINETYPE_BEZIER
    obj.SetName(Name)


    for i in range(0,N):
       dx=y
       dy=z
       dz=-z-(b-a+a*x*x)*y-b*x
       x=x + Delta * dx
       y=y + Delta * dy
       z=z + Delta * dz
       obj.SetPoint(i, c4d.Vector(x*Faktor, y*Faktor, z*Faktor))

    obj.Message(c4d.MSG_UPDATE)
    return obj

def main():
    slinobj = CreateSplineObject()

    doc.InsertObject(slinobj, None, None, True)

    c4d.EventAdd()

if __name__=='__main__':
    main()