 Phase-mixed radial probability of a Kepler orbit

For a Kepler orbit with pericentre distance $q$, apocentre distance $Q$, and semi-major axis $a=(q+Q)/2$, a uniform mean anomaly givesThe two radial passages per orbital period are included, and ${\int }_q^Q{p}_r(r)dr=1$.