statrs::function::harmonic
pub fn harmonic(t: u64) -> f64
Computes the t-th harmonic number
t
Returns 1 as a special case when t == 0
1
t == 0