Orthogonal Polynomial Series in Numpy: Playing with Legendre and Chebyshev Series