Mandelbrot set: Difference between revisions

m
(Solution in q)
Line 8,980:
 
for k in range(n): # you can also try range(2*n) due to rebasing
M, R = Z.real ** 2 + Z.imag ** 2 < r ** 2, Z.real ** 2 + Z.imag ** 2 < E.real ** 2 + E.imag ** 2
E[R], I[R] = Z[R], J[R] # rebase if z is closer to zero
M = Z.real ** 2 + Z.imag ** 2 < r ** 2
E[M], I[M] = (2 * S[I[M]] + E[M]) * E[M] + C[M], I[M] + 1
Z[M], dZ[M] = S[I[M]] + E[M], 2 * Z[M] * dZ[M] + 1
305

edits