Algorithme

x0 = .01 y0 = .003 N = INT(RND * 10) gamma = .9 + .5 * RND omega = RND - .5 x = x0 y = y0 iterates = 1 lambda = -1.8 alpha = 2 beta = 0 FOR nn = 1 TO 30000 skip first 100 iterations PlotPoint x, y zzbar = x * x + y * y zreal = x zimag = y FOR i = 1 TO N - 2 za = zreal * x - zimag * y zb = zimag * x + zreal * y zreal = za zimag = zb NEXT zn = x * zreal - y * zimag P = lambda + alpha * zzbar + beta * zn xnew = P * x + gamma * zreal - omega * y ynew = P * y - gamma * zimag + omega * x x = xnew y = ynew NEXT