`scipy.special._orthogonal:gegenbauer`

Signature

def gegenbauer ( n , alpha , monic = False )

Gegenbauer (ultraspherical) polynomial.

Defined to be the solution of

(1 - x^{2}) \frac{d ^{2}}{d x ^{2}} C_{n}^{(α)} - (2 α + 1) x \frac{d}{d x} C_{n}^{(α)} + n (n + 2 α) C_{n}^{(α)} = 0

for $α > - 1/2$ ; $C_{n}^{(α)}$ is a polynomial of degree $n$ .

n : int: Degree of the polynomial.
alpha : float: Parameter, must be greater than -0.5.
monic : bool, optional: If True, scale the leading coefficient to be 1. Default is False.

The polynomials $C_{n}^{(α)}$ are orthogonal over $[- 1, 1]$ with weight function $(1 - x^{2})^{(α - 1/2)}$ .

import numpy as np
from scipy import special
import matplotlib.pyplot as plt

✓

We can initialize a variable ``p`` as a Gegenbauer polynomial using the `gegenbauer` function and evaluate at a point ``x = 1``.

p = special.gegenbauer(3, 0.5, monic=False)
p

✓

p(1)

✗

To evaluate ``p`` at various points ``x`` in the interval ``(-3, 3)``, simply pass an array ``x`` to ``p`` as follows:

x = np.linspace(-3, 3, 400)
y = p(x)

✓

We can then visualize ``x, y`` using `matplotlib.pyplot`.

fig, ax = plt.subplots()

✓

ax.plot(x, y)
ax.set_title("Gegenbauer (ultraspherical) polynomial of degree 3")
ax.set_xlabel("x")
ax.set_ylabel("G_3(x)")

✗

plt.show()

✓