Le Tore Fractal de  KAM


Algorithme

Pick a [1 - 100]

Pick F close to 1.0

Pick OrbitIteration [e.g. 150]

Pick OrbitIncrement [e.g. 0.005]

Pick Orbit [e.g. -1.2]

Pick MaxOrbit [e.g. 8]

orbit = - 1.2

do

orbit = orbit + OrbitIncrement

x = orbit / 3

y = orbit / 3

do

xnew = x * COS(a) + (x * x - y) * SIN(a)

ynew = x * SIN(a) - (x * x - F * y) * COS(a)

x = xnew

y = ynew

Skip first 100 iterations

Plot x, y

loop OrbitIteration Times

loop until orbit = MaxOrbit


[Image]