extendiendo mandelbrot para generar julia

trabajando en un proyecto que requiere que use el mismo código, tenga en cuenta en el mismo archivo para generar el conjunto de mandelbrot y los conjuntos de julia, tengo un conjunto de mandelbrot que funciona pero puedo ver cómo extender al conjunto de julia usando el mismo código. tal vez no entiendo las diferencias entre? alguien puede elaborar

import numpy as np
import matplotlib.pyplot as plt
import math



def Mandelbrot(zmin, zmax, m, n, tmax=256):

    xs = np.linspace(zmin, zmax, n)
    ys = np.linspace(zmin, zmax, m)
    X, Y = np.meshgrid(xs, ys)


    Z = X + 1j * Y
    C = np.copy(Z)
    M = np.ones(Z.shape) * tmax

    for t in xrange(tmax):
        mask = np.abs(Z) <= 2.
        Z[ mask] = Z[mask]**2 + C[mask]
        M[-mask] -= 1.
    return M

list=Mandelbrot(-2,2,500,500)
plt.imshow(list.T, extent=[-2, 1, -1.5, 1.5])
plt.gray()
plt.savefig('mandelbrot.png')

Respuestas a la pregunta(2)

Su respuesta a la pregunta