This routine computes the Polygamma function, psi_n(x).
function polygamma extends Modelica.Icons.Function; input Integer n; input Real x; output Real y; end polygamma;