Function phase_shift 

pub fn phase_shift(l: u32, rho: f64) -> f64

Hard-sphere phase shift φ_l(ρ).

Reference: SAMMY rml/mrml07.f Sinsix function (lines 254-342)