sqrt(2 \pi) exp(InverseCDF[NormalDistribution[0, 1], p]^2 / 2)