Legendre Polynomials Python

Question

I am trying to calculate the 720th Legendre polynomial with scipy:

>>>from scipy.special import legendre 
>>>print(legendre(720))

it prints the powers correctly but returns NAN for the coefficients how can I fix this?

Answer 1

It's probably a float overflow error . If you use the monic keyword of the legendre function you get non- nan coefficients:

>>> from scipy.special import legendre  
>>> print(legendre(720, monic=True))

720 sounds like a high number, maybe it's more stable to use the eval_legendre function.