pub fn radial_probability(n: u32, l: u32, r: f64) -> f64
Radial probability density r^2 |R_nl(r)|^2 (probability of finding electron at radius r).