This routine computes the Debye function of order 2, D2(x).
function debye_D2 extends Modelica.Icons.Function; input Real x; output Real y; end debye_D2;