This routine computes the Legendre polynomial of order l, Pl(x).
function legendre_Pl extends Modelica.Icons.Function; input Integer l; input Real x; output Real y; end legendre_Pl;