def _particle_func(self, coords, pos, wid):
dx, dy, dz = [c - p for c,p in zip(coords, pos)]<EOL>dr2 = dx*dx + dy*dy + dz*dz<EOL>return np.exp(-dr2/(<NUM_LIT:2>*wid*wid))<EOL>
Draws a gaussian, range is (0,1]. Coords = [3,n]