This routine computes the Debye function of order 1, D1(x).
function debye_D1 extends Modelica.Icons.Function; input Real x; output Real y; end debye_D1;