Function probability_density_3d 

pub fn probability_density_3d(
    n: u32,
    l: u32,
    m: i32,
    r: f64,
    theta: f64,
    phi: f64,
) -> f64

Probability density |psi|^2 at a point.