6.1 Power Series Method
For an ODE y"′+p(x)y′+q(x)y=0 where p and q are analytic near x0Substitute the Power series y=∑n=0∞an(x−x0)n and solve for the coefficients.
6.2 Ordinary and Regular Singular Points
x0 is an ordinary point if p and q are analytic at x0. It is a regular singular Point if (x−x0)p(x) and (x−x0)2q(x) are analytic at x0.
6.3 Frobenius Method
At a regular singular point x0=0Substitute y=∑n=0∞anxn+r. The indicial equation determines the possible values of r.
Theorem 6.1. If the roots r1≥r2 of the indicial equation differ by a non-integer, there Are two linearly independent solutions of the form xr1∑anxn and xr2∑bnxn.
6.4 Bessel’s Equation
Bessel’s equation of order ν:
x2y′′+xy′+(x2−ν2)y=0
For ν∈/ZThe solutions are Jν(x) and J−ν(x) (Bessel functions of the First kind). For ν=n∈NThe second solution is the Weber function Yn(x).
6.4b Worked Example: Higher-Order ODE
Problem. Solve y′′′−6y′′+11y′−6y=0.
Solution
Solution. Characteristic equation: r3−6r2+11r−6=0.
Trying r=1: 1−6+11−6=0. Factor: (r−1)(r2−5r+6)=(r−1)(r−2)(r−3)=0.
Roots: r=1,2,3 (three distinct real roots).
y=c1ex+c2e2x+c3e3x. ■
### 6.5 Worked Example: Power Series Method
Problem. Solve y′′−xy=0 (Airy’s equation) using power series about x0=0.
Solution
Solution. Since p(x)=0 and q(x)=−x are both analytic everywhere, x0=0 is an ordinary Point. Substitute y=∑n=0∞anxn:
y′=∑n=1∞nanxn−1, y′′=∑n=2∞n(n−1)anxn−2.
y′′−xy=∑n=2∞n(n−1)anxn−2−∑n=0∞anxn+1=0.
Shift indices: first sum ∑m=0∞(m+2)(m+1)am+2xmSecond sum ∑m=1∞am−1xm.
For m=0: 2⋅1⋅a2=0⟹a2=0.
For m≥1: (m+2)(m+1)am+2−am−1=0⟹am+2=(m+2)(m+1)am−1.
This gives: a3=6a0, a4=12a1, a5=20a2=0 a6=30a3=180a0Etc.
Since a2=0All a3k+2=0.
y(x)=a0(1+6x3+180x6+⋯)+a1(x+12x4+504x7+⋯).
These are the Airy functions Ai(x) and Bi(x) (up to normalization). ■
6.6 Worked Example: Frobenius Method
Problem. Solve 2xy′′+y′+xy=0 near x=0 using the Frobenius method.
Solution
Solution. Rewrite in standard form: y′′+2x1y′+21y=0.
x=0 is a regular singular point since xp(x)=1/2 and x2q(x)=x2/2 are analytic at 0.
Substitute y=∑n=0∞anxn+r, a0=0:
y′=∑n=0∞(n+r)anxn+r−1
y′′=∑n=0∞(n+r)(n+r−1)anxn+r−2
Substituting into 2xy′′+y′+xy=0:
∑n=0∞2(n+r)(n+r−1)anxn+r−1+∑n=0∞(n+r)anxn+r−1+∑n=0∞anxn+r+1=0
For n=0: [2r(r−1)+r]a0=0. Since a0=0: r(2r−2+1)=0⟹r(2r−1)=0.
Indicial equation: r=0 or r=1/2.
For general n≥1: [2(n+r)(n+r−1)+(n+r)]an+an−2=0
(n+r)(2n+2r−1)an=−an−2
an=−(n+r)(2n+2r−1)an−2
For r=0: an=−n(2n−1)an−2. Odd coefficients vanish (a1=0). Even: a2=−6a0 a4=120a0Etc.
For r=1/2: an=−(n+1/2)(2n)an−2=−n(2n+1)an−2.
y=C1∑k=0∞a2k(0)x2k+C2x1/2∑k=0∞a2k(1/2)x2k. ■
6.7 Worked Example: Bessel Functions
Problem. Find the first three nonzero terms of J0(x)The Bessel function of the first kind Of order zero, which satisfies x2y′′+xy′+x2y=0.
Solution
Solution. Here ν=0. The indicial equation gives r2=0 (repeated root r=0).
Substituting y=∑n=0∞anx2n (we can show only even powers appear):
y′=∑n=1∞2nanx2n−1, y′′=∑n=1∞2n(2n−1)anx2n−2.
x2y′′+xy′+x2y=∑n=1∞2n(2n−1)anx2n+∑n=1∞2nanx2n+∑n=0∞anx2n+2=0.
For n=0: a0 is free.
For the recurrence: 4n2an+an−1=0⟹an=−4n2an−1 for n≥1.
a1=−4a0, a2=64a0, a3=−2304a0.
Setting a0=1: J0(x)=1−4x2+64x4−2304x6+⋯. ■