Skip to content

Series Solutions

6.1 Power Series Method

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.

6.2 Ordinary and Regular Singular Points

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.

6.3 Frobenius Method

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.

6.4 Bessel’s Equation

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).

6.4b Worked Example: Higher-Order ODE

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

6.6 Worked Example: Frobenius Method

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

6.7 Worked Example: Bessel Functions

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