; function to return a density of nulls for given corrected dimensionless G and q

function pnull, gcz, q
Pnull=gcz*(q)^3
return, pnull

end