Compute the digamma function, also known as the psi function.
∞
/ exp (-t) exp (-xt)
psi (x) = |---------- - -------------- dt
/ t 1 - exp (-t)
0
Accuracy: The result is a tight enclosure.
psi (infsup (1)) ⇒ ans ⊂ [-0.57722, -0.57721]
Package: interval