pub fn f_trigamma(x: f64) -> f64
Computes the trigamma function ψ₁(x).
The trigamma function is the second derivative of the logarithm of the gamma function.