This routine computes the Debye function of order 4, D4(x).
function debye_D4 extends Modelica.Icons.Function; input Real x; output Real y; end debye_D4;