fn digamma(x: f64) -> f64
Digamma for x >= 1: DLMF 5.5.2 up to x >= 10, then DLMF 5.11.2.
x >= 1
x >= 10