from cmath import *
from pylab import *

def M_1(c):
    k=0
    z=0+0j
    while k<=100 and abs(z)<=2:
        z=z**2+c
        k=k+1
    if abs(z)<=2:
        m=True
    else:
        m=False
    return m

def M_2(Xmin,Xmax,Ymin,Ymax,N):
    for u in range(0,N+1):
        for v in range(0,N+1):
            re=Xmin+u*(Xmax-Xmin)/N
            im=Ymin+v*(Ymax-Ymin)/N
            c=re+im*1j
            if M_1(c):
                plot([re],[im],'.')
    show()
    return