This routine computes the Digamma function, psi(x) = G'(x)/G(x).
function digamma extends Modelica.Icons.Function; input Real x; output Real y; end digamma;