In [6]:
import matplotlib.pylab as plt
import numpy as np
%matplotlib inline

In [11]:
def create(positions,rnopoints,polarnopoints):
    radialarray=np.arange(rnopoints)**2
    polararray=np.arange(polarnopoints)/(2*np.pi)

In [12]:
position=np.array([0,0,0])
create(position,5,5)

In [ ]: