pub fn p_harmonic(n: u64, p: u32) -> f64
Returns the sum of the first n - 1 terms of the p-harmonic series