Bessel Functions

import numpy as np
import matplotlib.pyplot as plt

The Bessel equation of order \(p\) is

\[ x^2 y'' + xy' + (x^2 - p^2)y = 0 \]

where \(p \in \mathbb{R}\) is any constant however for simplicity, we are only going to consider \(p = n \geq 0\) for \(n=0,1,2,3,\dots\).