File:Douady rabbit, plotted with Matplotlb.svg

from pylab import *
from numpy import NaN
from matplotlib.colors import LogNorm

s = -.123, .745

res = .001    # grid resolution
max_it = 1000 # maximum number of iterations

figure(figsize = (10, 10))

def jul(a):
	z = a
	for n in range(1, max_it + 1):
		z = z**2 + s[0] + 1j * s[1]
		if abs(z) > 2:
			return n
	return NaN

X = arange(-2, 2 + res, res)
Y = arange(-1.5,  1.5 + res, res)
Z = zeros((len(Y), len(X)))

for iy, y in enumerate(Y):
	print (iy + 1, "of", len(Y))
	for ix, x in enumerate(X):
		Z[-iy - 1, ix] = jul(x + 1j * y)

save("douady", Z)	# save array to file

imshow(Z, cmap = plt.cm.viridis, interpolation = 'none', norm = LogNorm(),
  extent = (X.min(), X.max(), Y.min(), Y.max()))
xlabel("Re(z)")
ylabel("Im(z)")
axis((-1.5, 1.5, -1.5, 1.5))
savefig("douady_rabbit_python.svg")
show()