Skip to content

Series Solutions | Mathematics

For an ODE y"+p(x)y+q(x)y=0y"' + p(x)y' + q(x)y = 0 where pp and qq are analytic near x0x_0Substitute the Power series y=n=0an(xx0)ny = \sum_{n=0}^{\infty} a_n (x - x_0)^n and solve for the coefficients.

x0x_0 is an ordinary point if pp and qq are analytic at x0x_0. It is a regular singular Point if (xx0)p(x)(x - x_0)p(x) and (xx0)2q(x)(x - x_0)^2 q(x) are analytic at x0x_0.

At a regular singular point x0=0x_0 = 0Substitute y=n=0anxn+ry = \sum_{n=0}^{\infty} a_n x^{n + r}. The indicial equation determines the possible values of rr.

Theorem 6.1. If the roots r1r2r_1 \geq r_2 of the indicial equation differ by a non-integer, there Are two linearly independent solutions of the form xr1anxnx^{r_1}\sum a_n x^n and xr2bnxnx^{r_2}\sum b_n x^n.

Bessel’s equation of order ν\nu:

x2y+xy+(x2ν2)y=0x^2 y'' + xy' + (x^2 - \nu^2)y = 0

For νZ\nu \notin \mathbb{Z}The solutions are Jν(x)J_\nu(x) and Jν(x)J_{-\nu}(x) (Bessel functions of the First kind). For ν=nN\nu = n \in \mathbb{N}The second solution is the Weber function Yn(x)Y_n(x).

Problem. Solve y6y+11y6y=0y''' - 6y'' + 11y' - 6y = 0.

Solution

Solution. Characteristic equation: r36r2+11r6=0r^3 - 6r^2 + 11r - 6 = 0.

Trying r=1r = 1: 16+116=01 - 6 + 11 - 6 = 0. Factor: (r1)(r25r+6)=(r1)(r2)(r3)=0(r - 1)(r^2 - 5r + 6) = (r - 1)(r - 2)(r - 3) = 0.

Roots: r=1,2,3r = 1, 2, 3 (three distinct real roots).

y=c1ex+c2e2x+c3e3xy = c_1 e^x + c_2 e^{2x} + c_3 e^{3x}. \blacksquare

### 6.5 Worked Example: Power Series Method

Problem. Solve yxy=0y'' - xy = 0 (Airy’s equation) using power series about x0=0x_0 = 0.

Solution

Solution. Since p(x)=0p(x) = 0 and q(x)=xq(x) = -x are both analytic everywhere, x0=0x_0 = 0 is an ordinary Point. Substitute y=n=0anxny = \sum_{n=0}^{\infty} a_n x^n:

y=n=1nanxn1y' = \sum_{n=1}^{\infty} na_n x^{n-1}, y=n=2n(n1)anxn2y'' = \sum_{n=2}^{\infty} n(n-1)a_n x^{n-2}.

yxy=n=2n(n1)anxn2n=0anxn+1=0y'' - xy = \sum_{n=2}^{\infty} n(n-1)a_n x^{n-2} - \sum_{n=0}^{\infty} a_n x^{n+1} = 0.

Shift indices: first sum m=0(m+2)(m+1)am+2xm\sum_{m=0}^{\infty} (m+2)(m+1)a_{m+2} x^mSecond sum m=1am1xm\sum_{m=1}^{\infty} a_{m-1} x^m.

For m=0m = 0: 21a2=0    a2=02 \cdot 1 \cdot a_2 = 0 \implies a_2 = 0.

For m1m \geq 1: (m+2)(m+1)am+2am1=0    am+2=am1(m+2)(m+1)(m+2)(m+1)a_{m+2} - a_{m-1} = 0 \implies a_{m+2} = \frac{a_{m-1}}{(m+2)(m+1)}.

This gives: a3=a06a_3 = \frac{a_0}{6}, a4=a112a_4 = \frac{a_1}{12}, a5=a220=0a_5 = \frac{a_2}{20} = 0 a6=a330=a0180a_6 = \frac{a_3}{30} = \frac{a_0}{180}Etc.

Since a2=0a_2 = 0All a3k+2=0a_{3k+2} = 0.

y(x)=a0(1+x36+x6180+)+a1(x+x412+x7504+)y(x) = a_0\left(1 + \frac{x^3}{6} + \frac{x^6}{180} + \cdots\right) + a_1\left(x + \frac{x^4}{12} + \frac{x^7}{504} + \cdots\right).

These are the Airy functions Ai(x)\mathrm{Ai}(x) and Bi(x)\mathrm{Bi}(x) (up to normalization). \blacksquare

Problem. Solve 2xy+y+xy=02xy'' + y' + xy = 0 near x=0x = 0 using the Frobenius method.

Solution

Solution. Rewrite in standard form: y+12xy+12y=0y'' + \frac{1}{2x}y' + \frac{1}{2}y = 0.

x=0x = 0 is a regular singular point since xp(x)=1/2xp(x) = 1/2 and x2q(x)=x2/2x^2 q(x) = x^2/2 are analytic at 00.

Substitute y=n=0anxn+ry = \sum_{n=0}^{\infty} a_n x^{n+r}, a00a_0 \neq 0:

y=n=0(n+r)anxn+r1y' = \sum_{n=0}^{\infty} (n+r)a_n x^{n+r-1}

y=n=0(n+r)(n+r1)anxn+r2y'' = \sum_{n=0}^{\infty} (n+r)(n+r-1)a_n x^{n+r-2}

Substituting into 2xy+y+xy=02xy'' + y' + xy = 0:

n=02(n+r)(n+r1)anxn+r1+n=0(n+r)anxn+r1+n=0anxn+r+1=0\sum_{n=0}^{\infty} 2(n+r)(n+r-1)a_n x^{n+r-1} + \sum_{n=0}^{\infty} (n+r)a_n x^{n+r-1} + \sum_{n=0}^{\infty} a_n x^{n+r+1} = 0

For n=0n = 0: [2r(r1)+r]a0=0[2r(r-1) + r]a_0 = 0. Since a00a_0 \neq 0: r(2r2+1)=0    r(2r1)=0r(2r - 2 + 1) = 0 \implies r(2r - 1) = 0.

Indicial equation: r=0r = 0 or r=1/2r = 1/2.

For general n1n \geq 1: [2(n+r)(n+r1)+(n+r)]an+an2=0[2(n+r)(n+r-1) + (n+r)]a_n + a_{n-2} = 0

(n+r)(2n+2r1)an=an2(n+r)(2n + 2r - 1)a_n = -a_{n-2}

an=an2(n+r)(2n+2r1)a_n = -\frac{a_{n-2}}{(n+r)(2n + 2r - 1)}

For r=0r = 0: an=an2n(2n1)a_n = -\frac{a_{n-2}}{n(2n-1)}. Odd coefficients vanish (a1=0a_1 = 0). Even: a2=a06a_2 = -\frac{a_0}{6} a4=a0120a_4 = \frac{a_0}{120}Etc.

For r=1/2r = 1/2: an=an2(n+1/2)(2n)=an2n(2n+1)a_n = -\frac{a_{n-2}}{(n+1/2)(2n)} = -\frac{a_{n-2}}{n(2n+1)}.

y=C1k=0a2k(0)x2k+C2x1/2k=0a2k(1/2)x2ky = C_1 \sum_{k=0}^{\infty} a_{2k}^{(0)} x^{2k} + C_2 x^{1/2} \sum_{k=0}^{\infty} a_{2k}^{(1/2)} x^{2k}. \blacksquare

Problem. Find the first three nonzero terms of J0(x)J_0(x)The Bessel function of the first kind Of order zero, which satisfies x2y+xy+x2y=0x^2 y'' + xy' + x^2 y = 0.

Solution

Solution. Here ν=0\nu = 0. The indicial equation gives r2=0r^2 = 0 (repeated root r=0r = 0).

Substituting y=n=0anx2ny = \sum_{n=0}^{\infty} a_n x^{2n} (we can show only even powers appear):

y=n=12nanx2n1y' = \sum_{n=1}^{\infty} 2n a_n x^{2n-1}, y=n=12n(2n1)anx2n2y'' = \sum_{n=1}^{\infty} 2n(2n-1) a_n x^{2n-2}.

x2y+xy+x2y=n=12n(2n1)anx2n+n=12nanx2n+n=0anx2n+2=0x^2 y'' + xy' + x^2 y = \sum_{n=1}^{\infty} 2n(2n-1)a_n x^{2n} + \sum_{n=1}^{\infty} 2n a_n x^{2n} + \sum_{n=0}^{\infty} a_n x^{2n+2} = 0.

For n=0n = 0: a0a_0 is free.

For the recurrence: 4n2an+an1=0    an=an14n24n^2 a_n + a_{n-1} = 0 \implies a_n = -\frac{a_{n-1}}{4n^2} for n1n \geq 1.

a1=a04a_1 = -\frac{a_0}{4}, a2=a064a_2 = \frac{a_0}{64}, a3=a02304a_3 = -\frac{a_0}{2304}.

Setting a0=1a_0 = 1: J0(x)=1x24+x464x62304+J_0(x) = 1 - \frac{x^2}{4} + \frac{x^4}{64} - \frac{x^6}{2304} + \cdots. \blacksquare

flowchart TD
A[6_Series Solutions] --> B[Key Concepts]
A --> C[Core Principles]
A --> D[Practical Applications]
B --> E[Fundamental definitions]
C --> F[Design patterns]
D --> G[Real-world usage]

When an ODE has no elementary closed-form solution, power series methods let you build a solution term by term. At an ordinary point, you can expand the solution as a Taylor series and substitute into the ODE to find recurrence relations for the coefficients. At a regular singular point, the solution may have a factor of xrx^r multiplied by a power series — the Frobenius method. The indicial equation determines the possible exponents rr. Bessel functions, Legendre polynomials, and other special functions all arise from this process. The radius of convergence is determined by the distance to the nearest singular point of the ODE, linking local analyticity to global solution behaviour.

Mistake 1: Confusing ordinary and regular singular points. An ordinary point is where p(x)p(x) and q(x)q(x) are analytic. A regular singular point is where (xx0)p(x)(x - x_0)p(x) and (xx0)2q(x)(x - x_0)^2 q(x) are analytic. An irregular singular point is where neither condition holds. Misclassifying the type of singular point leads to incorrect solution methods.

Mistake 2: Forgetting to check the indicial equation roots. When using the Frobenius method, the indicial equation determines the possible values of rr. If the roots differ by an integer, the second solution may involve a logarithmic term. Forgetting to check the roots can lead to missing solutions.

Mistake 3: Assuming that power series solutions always converge. Power series solutions converge only within the radius of convergence, which is determined by the distance to the nearest singular point. Do not assume that a power series solution converges everywhere; always check the radius of convergence.

Mistake 4: Forgetting that the Frobenius method requires regular singular points. The Frobenius method is applicable only at regular singular points. At irregular singular points, the method fails and other techniques (such as asymptotic expansions) are needed. Always verify that the singular point is regular before applying the Frobenius method.

Mistake 5: Confusing Bessel functions of the first and second kind. Bessel functions of the first kind Jν(x)J_\nu(x) are analytic at x=0x = 0, while Bessel functions of the second kind Yν(x)Y_\nu(x) are singular at x=0x = 0. Do not confuse the two; they serve different purposes in solving Bessel’s equation.

This section provides detailed coverage of advanced concepts, including full derivations, proofs, and extended examples.

Complete mathematical derivations and proofs are provided where appropriate. Each step is explained to ensure understanding of the underlying reasoning.

Advanced examples demonstrate the application of concepts to complex problems. These examples go beyond standard exam questions to develop deeper understanding.

This material connects to current research and advanced applications in the field. Understanding these connections provides context for the study material.

Ensure you have mastered the prerequisite material before attempting this advanced content.