Many important differential equations, especially in physics and engineering, don't have solutions you can write down with simple functions like polynomials or exponentials. For these, we often turn to a powerful technique: finding the solution as an infinite series. This approach allows us to solve complex equations like the Legendre and Bessel equations by representing their solutions as a power series.

This method works particularly well when we're looking for a solution around a specific point, $x_0$. If the coefficients of the differential equation are well-behaved (analytic) at that point, we can use a standard power series. But things get more interesting when we have a regular singular point, a type of mild singularity where standard power series methods fail. This is where a more general technique comes into play.

The Frobenius method is a way to find a series solution for a second-order linear ODE around a regular singular point, say $x=0$. Instead of assuming the solution is a simple power series, we assume a slightly more general form.

The process involves substituting this series form into the differential equation. After taking derivatives and aligning the powers of $x$, we can solve for the unknown exponent $r$ and the coefficients $a_n$. The very first step, which involves the lowest power of $x$, is the most crucial.

Legendre's equation is a classic example that appears when solving problems with spherical symmetry. It's written as:

By substituting the series $y = \sum a_n x^n$ into the equation and solving for the coefficients, we find a recurrence relation. This relation links each coefficient to the one two steps before it ($a_{n+2}$ in terms of $a_n$).

A fascinating thing happens if the constant $p$ is a non-negative integer. Look at the numerator of the recurrence relation, $(n-p)$. When $n$ becomes equal to $p$, the numerator becomes zero. This means $a_{p+2}$ will be zero, and so will all subsequent coefficients in that series ($a_{p+4}, a_{p+6}, \dots$). One of the series solutions terminates and becomes a polynomial. These are the famous Legendre polynomials.

Bessel's equation is another cornerstone of mathematical physics, often appearing in problems involving cylindrical symmetry, like the vibrations of a circular drumhead or heat flow in a cylinder.

After substituting the Frobenius series $y = \sum a_n x^{n+r}$ into Bessel's equation, we find the indicial equation by looking at the coefficient of the lowest power of $x$.

The nature of the solutions depends on the value of $p$. When we work through the recurrence relation, we get solutions called Bessel functions. The first solution, corresponding to the root $r=p$, is called a Bessel function of the first kind, denoted $J_p(x)$.

For the series to be a valid solution, it must converge. Using tests like the ratio test, we find that the series solutions for both the Legendre and Bessel equations converge for all finite values of $x$. For Legendre's equation, the interval of convergence is typically $|x|<1$, related to the singular points at $x=\pm 1$.

The Frobenius method provides a systematic way to tackle differential equations with regular singular points, unlocking series solutions for fundamental equations like Bessel's and extending the power of series methods beyond ordinary points.