scipy.special.jacobi

scipy.special.jacobi(n, alpha, beta, monic=False)

Returns the nth order Jacobi polynomial, P^(alpha,beta)_n(x) orthogonal over [-1,1] with weighting function (1-x)**alpha (1+x)**beta with alpha,beta > -1.