Mandelbrot set: Difference between revisions

Content added Content deleted
(Added output. Added a graphical version.)
m (→‎{{header|Python}}: Removed an optimization to make the code more readable.)
Line 7,062: Line 7,062:
u, v = 0.0, 0.0
u, v = 0.0, 0.0
for k in range(n):
for k in range(n):
u, v = u ** 2 - v ** 2 + a, 2 * u * v + b
if u ** 2 + v ** 2 < r ** 2:
if not u ** 2 + v ** 2 < r ** 2:
u, v = u ** 2 - v ** 2 + a, 2 * u * v + b
T[i][j] = k + 1
else:
break
break
T[i][j] = k + 1


plt.imshow(T, cmap=plt.cm.twilight_shifted)
plt.imshow(T, cmap=plt.cm.twilight_shifted)