fitting by jacobi polynomials and get the coefficients

Question

I\'m trying to fit a curve with jacobi polynomials Pn(alpha,beta,x) using python3. It is possible to fit by other orthogonal polynomials using numpy such as legendre, laguer