m = 500
g = fltarr(m)
for i = 0,(m-1) do begin $
z = (1./10.)*findgen(m)& $
g[i] = Pn(z[i],0) & $
endfor

;This program evaluates the density of null points for different heights and produces an array of this. This program is /very/ slow.

end