SciPy

scipy.special.gegenbauer

scipy.special.gegenbauer(n, alpha, monic=0)[source]

Return the nth order Gegenbauer (ultraspherical) polynomial, C^(alpha)_n(x), orthogonal over [-1,1] with weighting function (1-x**2)**(alpha-1/2) with alpha > -1/2