This routine computes the Debye function of order 3, D3(x).
function debye_D3 extends Modelica.Icons.Function; input Real x; output Real y; end debye_D3;