In [1]:

    

def julia(z, c, maxiter=200):
    for n in xrange(maxiter):
        if abs(z) > 4:
            return n-1
        
        z = z*z + c
    
    return maxiter


    

In [8]:

    

import numpy as np


    

In [ ]:

    

#@time m = [ uint8(julia(complex(r,i), complex(-.06,.67))) for i=1:-.002:-1, r=-1.5:.002:1.5 ];


    

In [20]:

    

#X, Y = np.ogrid[1:-1:-.002, -1.5:1.5:.002]
#m = [[ int(julia(complex(r,i), complex(-.06,.67))) for i in np.arange(1, 1,-.002)] 
#     for r in np.arange(-1.5, 1.5, 0.002)];


    

In [23]:

    

%%time
for i in np.r_[1:-1:-.002]:
    for r in np.r_[-1.5:1.5:0.002]:
        julia(complex(r,i), complex(-.06, .67))


    

    




CPU times: user 16.4 s, sys: 99.4 ms, total: 16.5 s
Wall time: 16.5 s

In [ ]: