def cos_julia_ex():
c = pi/2*(1.0 + 0.6j)
xB = (-1.5, 1.5)
yB = (-1.5, 1.5)
start_time = time.time()
jul = julia(c, xB, yB, cosine, args=0, width=5, height=5, maxiter=1000, dpi=300)
image(jul, cmap=plt.cm.plasma, filename='cos_julia_ex', gamma=1.0)
print('calculation took %s seconds ' % np.round((time.time() - start_time), 3))
def julia_ex0():
c = 1j
xB = (-1.5, 1.5)
yB = (-1.5, 1.5)
start_time = time.time()
jul = julia(c, xB, yB, power, args=2, width=5, height=5, maxiter=500, dpi=300)
image(jul, cmap=plt.cm.viridis, filename='julia_ex0', gamma=0.6)
print('calculation took %s seconds ' % np.round((time.time() - start_time), 3))
def julia_ex1():
c = -0.391 - 0.587j
xB = (-1.5, 1.5)
yB = (-1.5, 1.5)
start_time = time.time()
jul = julia(c, xB, yB, power, args=2, width=5, height=5, maxiter=1000, dpi=300)
image(jul, cmap=plt.cm.gist_ncar, filename='julia_ex1', gamma=0.6)
print('calculation took %s seconds ' % np.round((time.time() - start_time), 3))
def mag_julia_ex1():
c = 2.0 + 0.80j
xB = (-1.0, 1.0)
yB = (-1.0, 1.0)
start_time = time.time()
jul = julia(c, xB, yB, magnetic_2, args=0, width=5, height=5, maxiter=1000, dpi=300)
image(jul, cmap=plt.cm.gist_ncar, filename='mag_julia_ex1', gamma=1.0)
print('calculation took %s seconds ' % np.round((time.time() - start_time), 3))