Mandelbrot set: Difference between revisions

Content added Content deleted
m (→‎Mercator zooms of the Mandelbrot set using NumPy: Do not transpose the Mercator map.)
Line 7,155: Line 7,155:


p0 = -0.41031 + 0.61021 * 1j # common center of all circles
p0 = -0.41031 + 0.61021 * 1j # common center of all circles
r0 = 1.0 # radius of the largest circle
r0 = 1.0 # radius of the largest circle


x = np.linspace(0, 2 * np.pi, d, endpoint=False)
x = np.linspace(0, 2 * np.pi, d, endpoint=False)
Line 7,179: Line 7,179:


A, B = C.real, C.imag
A, B = C.real, C.imag
S = (150 * D) ** 2
S = (100 * D) ** 2


fig, ax = plt.subplots(3, 2, figsize=(14, 21), dpi=100)
fig, ax = plt.subplots(3, 2, figsize=(10, 15), dpi=100)
ax[0, 0].scatter(A[0:200], B[0:200], s=S[0:200], c=T[0:200], cmap=plt.cm.twilight_shifted, vmin=0, vmax=n)
ax[0, 0].scatter(A[0:200], B[0:200], s=S[0:200], c=T[0:200], cmap=plt.cm.twilight_shifted, vmin=0, vmax=n)
ax[0, 1].scatter(A[100:300], B[100:300], s=S[0:200], c=T[100:300], cmap=plt.cm.twilight_shifted, vmin=0, vmax=n)
ax[0, 1].scatter(A[100:300], B[100:300], s=S[0:200], c=T[100:300], cmap=plt.cm.twilight_shifted, vmin=0, vmax=n)