Fixed Legendre Polynomials - Factorial now returns a double, which can be larger than an int; Fully tested l_max = 10, though l_max as large as 170 may be possible (based on factorial(170)<DBL_MAX<factorial(171))