; 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