pub fn gen_harmonic(n: u64, m: f64) -> f64
Computes the generalized harmonic number of order n of m e.g. (1 + 1/2^m + 1/3^m + ... + 1/n^m)
n
m
(1 + 1/2^m + 1/3^m + ... + 1/n^m)
Returns 1 as a special case when n == 0
1
n == 0