from random import *
from math import *
from pylab import *
def Julia(P,N):
    LX =[]
    LY =[]
    for i in range(P):
        z0 = random()*4-2+1j*(random()*4-2)
        z=z0
        k=0
        while (abs(z)<2 and k<N):
            k = k+1
            z = z**2-1
        if k == N:
            LX.append(z0.real)
            LY.append(z0.imag)
    scatter(LX,LY)
    show()