Ordinary Points Lecture Slides Ordinary Points
Power Series Solutions
Lecture I
The Variable Coefficient Challenge
Beyond constant coefficients
Consider the linear homogeneous second-order ODE:
\[ P(x)y'' + Q(x)y' + R(x)y = 0 \]
If \(P(x)\), \(Q(x)\), and \(R(x)\) are polynomials, we often lack elementary solutions. We seek an analytic solution near a point \(x_0\).
Definition: Ordinary Point
A point \(x_0\) is an ordinary point if \(P(x_0) \neq 0\).
Near an ordinary point, the functions \(p(x) = Q(x)/P(x)\) and \(q(x) = R(x)/P(x)\) are analytic (have Taylor series expansions).
The Ansatz and Derivatives
Assume a solution of the form:
\[ y(x) = \sum_{n=0}^{\infty} a_n (x-x_0)^n \]
First Derivative
\[ y' = \sum_{n=1}^{\infty} n a_n (x-x_0)^{n-1} \]
Second Derivative
\[ y'' = \sum_{n=2}^{\infty} n(n-1) a_n (x-x_0)^{n-2} \]
"The shift of indices is our primary tool for aligning powers of \(x\) to collect coefficients."
Deriving Recurrence Relations
1. Substitute series into ODE.
This yields a recurrence relation , expressing \(a_{n+k}\) in terms of previous coefficients.
Example: Airy's Equation
\(y'' - xy = 0\)
Leads to: \(a_{n+2} = \frac{a_{n-1}}{(n+2)(n+1)}\)
Fundamental Results
There exist two linearly independent series solutions centered at \(x_0\).
The radius of convergence \(\rho\) is at least as large as the distance from \(x_0\) to the nearest singularity of \(p(x)\) or \(q(x)\) in the complex plane.
General Solution: \(y = c_1 y_1(x) + c_2 y_2(x)\)
Recurrence Relation Workshop Ordinary Points
Differential Equations | Workshop 01
Student:
Date:
"In this workshop, we transition from the theory of analytic functions to the mechanics of recurrence. We will treat ordinary points—those where the leading coefficient is non-zero—and construct two fundamental solutions."
1
Analysis of Airy's Equation
Consider the differential equation \(y'' - xy = 0\) near the ordinary point \(x_0 = 0\).
A. Substitution and Index Shifting
Substitute \(y = \sum_{n=0}^{\infty} a_n x^n\) and shift the indices to express the equation as a single power series \(\sum_{n=0}^{\infty} [\dots] x^n = 0\).
B. The Recurrence Relation
Extract the recurrence relation for \(a_n\). For what values of \(n\) is this relation valid?
C. Fundamental Solutions
Derive the first four non-zero terms for two linearly independent solutions, \(y_1(x)\) and \(y_2(x)\), such that \(y_1(0)=1, y_1'(0)=0\) and \(y_2(0)=0, y_2'(0)=1\).
2
Hermite Equation & Polynomials
The Hermite differential equation is given by \(y'' - 2xy' + 2\lambda y = 0\), where \(\lambda\) is a constant.
A. Recurrence Analysis
Derive the recurrence relation relating \(a_{n+2}\) to \(a_n\).
Theoretical Inquiry: Polynomial Solutions
Under what condition on \(\lambda\) does the series solution terminate, resulting in a polynomial solution? If \(\lambda = 2\), find the corresponding Hermite polynomial \(H_2(x)\) normalized such that the leading coefficient is 4.
B. Radius of Convergence
Determine the radius of convergence for the Hermite series solutions when \(\lambda\) is not an integer. Justify your answer using the distance to the nearest singularity.
Ordinary Points Teaching Notes Ordinary Points
Teaching Notes & Facilitation Guide
Core Objectives
Students should grasp that power series are not just approximations, but exact representations of functions near ordinary points.
Mastery of index shifting (\(n \to n+k\)) is the technical bottleneck of this lesson; ensure students visualize the shift on a number line.
Understanding the radius of convergence in the complex plane provides deep insight into why solutions fail.
Common Student Errors
Technical Errors
Forgetting to pull powers of \(x\) outside the summation sign before shifting.
Inconsistent starting indices for the combined sum (missing the \(a_0, a_1\) standalone terms).
Conceptual Errors
Assuming the series always converges everywhere if the point is ordinary.
Mistaking a recurrence relation for the final solution.
Deep Inquiry Questions
"Why does Fuchs' Theorem guarantee that the radius of convergence is at least the distance to the nearest singularity?"
Discussion guide: Connect back to the Taylor series of \(p(z)\) and \(q(z)\) and the properties of analytic continuation.
"If we find a polynomial solution (like in Hermite's equation), what does that imply about the second linearly independent solution?"
Discussion guide: It must be an infinite series, as the polynomial behavior only occurs when the recurrence 'breaks' in one branch.
Workshop Key Points
Airy Equation: Recurrence is \(a_{n+2} = a_{n-1} / [(n+2)(n+1)]\). Note that \(a_2 = 0\). This leads to terms in groups of 3.
Hermite Equation: Recurrence is \(a_{n+2} = \frac{2n - 2\lambda}{(n+2)(n+1)} a_n\). Termination occurs when \(2n = 2\lambda\) for some even or odd branch.
Frobenius Method Lecture Slides Singular Points
The Method of Frobenius
Lecture II
Classification of Singularities
Standard Form: \(y'' + p(x)y' + q(x)y = 0\)
Regular Singular Point
\(x_0\) is Regular if:
\((x-x_0)p(x)\) is analytic at \(x_0\)
\((x-x_0)^2 q(x)\) is analytic at \(x_0\)
Irregular Singular Point
If either limit fails to exist or the functions are not analytic, the point is Irregular . Series solutions are significantly more complex or may not exist.
Limit \( (x \to x_0) \) check is essential.
The Frobenius Ansatz
Generalizing the power series to handle poles:
\[ y(x) = x^r \sum_{n=0}^{\infty} a_n x^n = \sum_{n=0}^{\infty} a_n x^{n+r} \]
Step 1
Substitute \(y, y', y''\) into the ODE. \(a_0 \neq 0\).
Step 2
Solve the Indicial Equation for the possible values of \(r\).
Step 3
Derive recurrence relations for coefficients \(a_n\).
The Indicial Equation
The lowest power of \(x\) (coefficient of \(a_0\)) gives the indicial equation :
\[ r(r-1) + p_0 r + q_0 = 0 \]
where \(p_0 = \lim_{x\to 0} xp(x)\) and \(q_0 = \lim_{x\to 0} x^2q(x)\).
Case Analysis
Case 1: \(r_1 - r_2 \neq \text{integer}\).Case 2: \(r_1 = r_2\).Case 3: \(r_1 - r_2 = \text{integer}\). Indicial Equation Practice Singularities & Frobenius
Advanced Differential Equations | Problem Set 02
Scholar:
01
Singularity Classification
Classify the singular points (at finite \(x\)) for the following differential equations. For each, determine if the point is Regular or Irregular .
A. \(x^2(1-x)y'' + xy' + y = 0\)
Work Area: Points \(x=0\) and \(x=1\)
B. \(x^2 y'' + e^x y' + y = 0\)
Work Area: Point \(x=0\)
02
The Indicial Equation
Consider the equation: \( 4xy'' + 2y' + y = 0 \) near \(x=0\).
Part A: Substitution
Assume \( y = \sum_{n=0}^{\infty} a_n x^{n+r} \). Substitute into the ODE and extract the coefficients for the lowest power of \(x\).
Part B: Roots
Determine the roots \(r_1\) and \(r_2\) of the indicial equation. Which Frobenius case is this?
Part C: Recurrence
Find the general recurrence relation for \(a_n\) in terms of \(a_{n-1}\).
03
Series Solution Construction
Using the larger root \(r_1\) from Problem 2, construct the first three non-zero terms of the series solution \(y_1(x)\).
Challenge Reflection
Observe the second root \(r_2\). In this specific case, does the second solution \(y_2(x)\) require a logarithmic term? Explain your reasoning based on the difference \(r_1 - r_2\).
Frobenius Answer Key Answer Key
Frobenius Method & Singularities
PROFESSOR USE ONLY
1. Singularity Classification
A. \(x^2(1-x)y'' + xy' + y = 0\)
x = 0: \(xp(x) = \frac{x}{x(1-x)} = \frac{1}{1-x} \to 1\) as \(x\to 0\). \(x^2q(x) = \frac{x^2}{x^2(1-x)} = \frac{1}{1-x} \to 1\). Regular Singular Point.x = 1: \((x-1)p(x) = \frac{(x-1)x}{x^2(1-x)} = \frac{-1}{x} \to -1\). \((x-1)^2q(x) = \frac{(x-1)^2}{x^2(1-x)} = \frac{1-x}{x^2} \to 0\). Regular Singular Point.
B. \(x^2 y'' + e^x y' + y = 0\)
x = 0: \(xp(x) = \frac{e^x}{x}\). As \(x \to 0\), this limit is \(\infty\). Irregular Singular Point.
2. The Indicial Equation
ODE: \(4xy'' + 2y' + y = 0\). Divide by \(4x\): \(y'' + \frac{1}{2x}y' + \frac{1}{4x}y = 0\).
\(p_0 = \lim_{x\to 0} x(\frac{1}{2x}) = 1/2\). \(q_0 = \lim_{x\to 0} x^2(\frac{1}{4x}) = 0\).
Indicial Eq: r(r-1) + (1/2)r + 0 = 0
Roots: r₁ = 1/2, r₂ = 0 . Case: \(r_1 - r_2 = 1/2\) (Not an integer). Case 1.
Recurrence: \(a_n = -\frac{a_{n-1}}{4(n+r)(n+r-1) + 2(n+r)} = -\frac{a_{n-1}}{2(n+r)(2n+2r-1)}\).
3. Series Solution (r = 1/2)
Using \(r = 1/2\), the recurrence becomes: \(a_n = -\frac{a_{n-1}}{2(n+1/2)(2n)} = -\frac{a_{n-1}}{n(2n+1)}\).
n=1: \(a_1 = -a_0 / (1 \cdot 3) = -a_0 / 3\)
n=2: \(a_2 = -a_1 / (2 \cdot 5) = a_0 / (3 \cdot 10) = a_0 / 30\)
Solution: \(y_1(x) = a_0 \sqrt{x} [ 1 - \frac{1}{3}x + \frac{1}{30}x^2 - \dots ]\)
Bessel Functions Lecture Slides Bessel Functions
Cylindrical Symmetry
Lecture III
Bessel's Equation
Arising from Laplacian in cylindrical coordinates
\[ x^2y'' + xy' + (x^2 - \nu^2)y = 0 \]
The parameter \(\nu\) (the order) can be any real or complex number. The point \(x=0\) is a regular singular point .
Classification
Applying the Frobenius Method at \(x=0\):
Indicial Equation: \(r^2 - \nu^2 = 0 \implies r = \pm \nu\)
First Solution: \(J_\nu(x)\) (Bessel function of the first kind)
Second Solution: \(Y_\nu(x)\) (Bessel function of the second kind)
The First Kind: \(J_\nu(x)\)
\[ J_\nu(x) = \sum_{m=0}^{\infty} \frac{(-1)^m}{m! \Gamma(m + \nu + 1)} \left( \frac{x}{2} \right)^{2m+\nu} \]
Properties
Analytic at \(x=0\) (for \(\nu \geq 0\)).
Oscillatory behavior for large \(x\).
Amplitude decays as \(x^{-1/2}\).
Convergence
The series converges for all \(x\) in the complex plane (entire function for integer \(\nu\)). For large \(x\), \(J_\nu(x) \approx \sqrt{\frac{2}{\pi x}} \cos(x - \frac{\nu\pi}{2} - \frac{\pi}{4})\).
The Second Kind: \(Y_\nu(x)\)
Also known as Weber functions or Neumann functions .
\[ Y_\nu(x) = \frac{J_\nu(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)} \]
For integer \(n\), we take the limit \(\nu \to n\).
Divergence at Origin
\(Y_\nu(x)\) is singular at \(x=0\). It approaches \(-\infty\) as \(x \to 0\).
"In physical problems requiring finite solutions at the origin (e.g., solid cylinders), the coefficient of \(Y_\nu\) must be zero."
Bessel Property Cheat Sheet Bessel Compendium
Essential Relations & Properties | Reference Sheet 03
Generating Function
\[ e^{\frac{x}{2}(t - 1/t)} = \sum_{n=-\infty}^{\infty} J_n(x) t^n \]
Recurrence Relations
\[ J_{\nu+1}(x) = \frac{2\nu}{x}J_\nu(x) - J_{\nu-1}(x) \]
\[ J_{\nu}'(x) = \frac{1}{2}[J_{\nu-1}(x) - J_{\nu+1}(x)] \]
\[ \frac{d}{dx}[x^\nu J_\nu(x)] = x^\nu J_{\nu-1}(x) \]
\[ \frac{d}{dx}[x^{-\nu} J_\nu(x)] = -x^{-\nu} J_{\nu+1}(x) \]
Integral Representations
\[ J_n(x) = \frac{1}{\pi} \int_{0}^{\pi} \cos(n\theta - x\sin\theta) d\theta \]
(For integer n)
Asymptotic Behavior
Small Argument (x &to; 0):
\[ J_\nu(x) \approx \frac{1}{\Gamma(\nu+1)} \left( \frac{x}{2} \right)^\nu \]
Large Argument (x &to; ∞):
\[ J_\nu(x) \approx \sqrt{\frac{2}{\pi x}} \cos\left( x - \frac{\nu\pi}{2} - \frac{\pi}{4} \right) \]
Orthogonality
Bessel functions of fixed order \(\nu\) satisfy an orthogonality condition on the interval \([0, a]\) with respect to the weight function \(w(x) = x\).
\[ \int_{0}^{a} x J_\nu(\alpha_i x) J_\nu(\alpha_j x) dx = 0 \quad (i \neq j) \]
where \(\alpha_i, \alpha_j\) are roots of \(J_\nu(\alpha a) = 0\).
Cylindrical Harmonics
Solutions to Helmholtz Eq in polar coords: \( \psi = [A J_n(kr) + B Y_n(kr)] e^{in\theta} \)
Circular Drumhead Problem Set Vibrating Drumheads
Differential Equations | Physics Case Study
Name:
Problem Context
Consider a circular membrane of radius \(a\) stretched across a frame. The displacement \(u(r, \theta, t)\) satisfies the wave equation \(\nabla^2 u = \frac{1}{c^2} \frac{\partial^2 u}{\partial t^2}\). In polar coordinates, assuming a separation of variables \(u = R(r)\Theta(\theta)T(t)\), the radial part leads to Bessel's equation.
01
The Radial Component
Given the boundary condition \(u(a, \theta, t) = 0\) (clamped edge), explain why the radial solution must be of the form \(R(r) = A J_n(k r)\) and why \(B Y_n(k r)\) is physically excluded.
02
Natural Frequencies
Let \(\alpha_{n,m}\) denote the \(m\)-th positive zero of \(J_n(x)\). Show that the natural frequencies of vibration for the drumhead are given by:
\[ f_{n,m} = \frac{c \alpha_{n,m}}{2\pi a} \]
03
Fundamental Mode vs. First Overtone
The fundamental mode corresponds to \(n=0, m=1\) (\(\alpha_{0,1} \approx 2.405\)). The first overtone with a nodal line corresponds to \(n=1, m=1\) (\(\alpha_{1,1} \approx 3.832\)).
A. Frequency Ratio
Calculate the ratio of the first overtone frequency to the fundamental frequency. Is it a harmonic (integer multiple)?
B. Nodal Visualization
Sketch the nodal lines for the (0,1) and (1,1) modes in the space provided. Indicate regions of opposing phase.
04
Total Displacement
Express the total displacement \(u(r, \theta, t)\) as a generalized Fourier-Bessel series for an arbitrary initial displacement \(f(r, \theta)\) at \(t=0\).
Legendre Polynomials Lecture Slides Legendre Polynomials
Spherical Symmetry
Lecture IV
Legendre's Differential Equation
Arising from Potential Theory in spherical coordinates
The standard form of Legendre's equation:
\[ (1-x^2)y'' - 2xy' + \ell(\ell+1)y = 0 \]
This equation appears when solving Laplace's equation \(\nabla^2 V = 0\) for systems with azimuthal symmetry, where \(x = \cos \theta\).
Key Observations
Points \(x = \pm 1\) are regular singular points .
Polynomial solutions \(P_\ell(x)\) exist only when \(\ell\) is a non-negative integer.
The second solution \(Q_\ell(x)\) diverges at \(x = \pm 1\) (singular at the poles).
Rodrigues' Formula
A compact operator representation for \(P_\ell(x)\):
\[ P_\ell(x) = \frac{1}{2^\ell \ell!} \frac{d^\ell}{dx^\ell} (x^2 - 1)^\ell \]
Order 0
\[ P_0(x) = 1 \]
Order 1
\[ P_1(x) = x \]
Order 2
\[ P_2(x) = \frac{1}{2}(3x^2 - 1) \]
Order 3
\[ P_3(x) = \frac{1}{2}(5x^3 - 3x) \]
Generating Function
\[ \frac{1}{\sqrt{1 - 2xt + t^2}} = \sum_{\ell=0}^{\infty} P_\ell(x) t^\ell \]
This relates the Legendre polynomials to the multipole expansion of the gravitational or electrostatic potential of a point charge.
Recurrence Relations
Three-Term Recurrence
\[ (\ell+1)P_{\ell+1}(x) = (2\ell+1)xP_\ell(x) - \ell P_{\ell-1}(x) \]
Derivative Relation
\[ P_{\ell+1}'(x) - P_{\ell-1}'(x) = (2\ell+1)P_\ell(x) \]
Rodrigues Formula Proof Sheet Rodrigues Analysis
Differential Equations | Proof Sheet 04
Researcher:
Theorem: Rodrigues' Formula
The Legendre polynomials \(P_n(x)\) which are solutions to \((1-x^2)y'' - 2xy' + n(n+1)y = 0\) can be generated by the formula: \[ P_n(x) = \frac{1}{2^n n!} \frac{d^n}{dx^n} (x^2 - 1)^n \]
01
Direct Verification
Use Rodrigues' formula to explicitly derive the Legendre polynomial \(P_3(x)\). Show every step of the differentiation process.
02
Proof of Orthogonality
Using Rodrigues' formula, prove that \(\int_{-1}^{1} P_n(x) P_m(x) dx = 0\) for \(n \neq m\). Hint: Assume \(m < n\) and use integration by parts \(n\) times.
03
The Normalization Constant
Evaluate the integral \(I = \int_{-1}^{1} [P_n(x)]^2 dx\) using Rodrigues' formula and integration by parts. Show that the result is \(\frac{2}{2n+1}\).
Potential Theory Discussion Guide Potential Theory
Legendre Polynomials Facilitation Guide
The Physical Motivation
Legendre polynomials are the natural language of azimuthal symmetry . In electrostatics, the potential \(V\) due to a point charge \(q\) at position \(d\) on the \(z\)-axis is: \[ V(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \frac{q}{|\mathbf{r} - \mathbf{d}|} \] Expanding this in terms of \(r\) and \(\theta\) (where \(\cos\theta = x\)) leads directly to the generating function: \[ \frac{1}{\sqrt{r^2 + d^2 - 2rd\cos\theta}} = \frac{1}{d} \sum_{\ell=0}^{\infty} \left(\frac{r}{d}\right)^\ell P_\ell(\cos\theta) \]
Key Instructional Pivots
Singularity Check
Emphasize that the series for \(P_\ell(x)\) generally diverges at \(x=\pm 1\) unless the series terminates. This is why we only care about integer \(\ell\) in physical space (the poles must be finite).
Parity Property
Students should observe that \(P_\ell(-x) = (-1)^\ell P_\ell(x)\). This is crucial for problems with hemispherical symmetry or integrating over half-spaces.
Seminar Questions
"Why does the second solution \(Q_\ell(x)\) always diverge at \(x=\pm 1\)? What does this imply about the physics of a solid sphere?"
Response: It implies that for any region containing the axis of symmetry (\(\theta=0\) or \(\pi\)), the coefficient of \(Q_\ell\) must vanish to maintain a finite potential.
"How does the generating function approach provide an advantage over the recurrence relation method?"
Response: The generating function captures all polynomials at once and reveals the deep connection to inverse-distance laws (Green's functions).
Note on Zeros
Remind students that \(P_\ell(x)\) has exactly \(\ell\) distinct zeros in the interval \((-1, 1)\). This inter-spacing of zeros is a precursor to the discussion on Sturm-Liouville eigenfunctions in the final lesson.
Sturm Liouville Lecture Slides Sturm-Liouville
Theory & Orthogonality
Lecture V
The Sturm-Liouville Problem
A unifying framework for special functions
The canonical equation for \(\mathcal{L}[y] + \lambda w(x)y = 0\):
\[ \frac{d}{dx} \left[ p(x) \frac{dy}{dx} \right] + [q(x) + \lambda w(x)]y = 0 \]
By defining the operator \(\mathcal{L}\), we can treat differential equations as eigenvalue problems in Hilbert space.
Weight Functions
Every special function has a characteristic weight \(w(x)\):
Bessel: \(w(x) = x\)
Legendre: \(w(x) = 1\)
Hermite: \(w(x) = e^{-x^2}\)
Chebyshev: \(w(x) = \frac{1}{\sqrt{1-x^2}}\)
Fundamental Theorems
Theorem 1: Real Eigenvalues
For a regular Sturm-Liouville problem, all eigenvalues \(\lambda_n\) are real and can be ordered: \(\lambda_1 < \lambda_2 < \dots < \lambda_n \to \infty\).
Theorem 2: Orthogonality
Eigenfunctions \(y_n, y_m\) corresponding to distinct eigenvalues \(\lambda_n, \lambda_m\) are orthogonal with respect to the weight function \(w(x)\).
\[ \int_a^b y_n(x) y_m(x) w(x) dx = 0 \]
Generalized Fourier Series
The eigenfunctions form a complete basis for \(L^2([a,b], w)\):
\[ f(x) = \sum_{n=1}^{\infty} c_n y_n(x) \]
\[ c_n = \frac{\int_a^b f(x) y_n(x) w(x) dx}{\int_a^b [y_n(x)]^2 w(x) dx} \]
This result justifies the expansion of any physical state (potential, displacement, temperature) into the special functions appropriate for the geometry of the problem.
Orthogonality Proof Challenge Orthogonality Challenge
Sturm-Liouville Theory | Proof Workshop 05
Scholar:
01
Operator Transformation
Every second-order linear ODE \(A(x)y'' + B(x)y' + C(x)y + \lambda D(x)y = 0\) can be put into the self-adjoint Sturm-Liouville form by multiplying by an integrating factor \(\mu(x)\). Show that for Bessel's equation of order \(\nu\):
\[ x^2y'' + xy' + (\lambda^2 x^2 - \nu^2)y = 0 \]
the Sturm-Liouville form is \(\frac{d}{dx}[x y'] + (-\frac{\nu^2}{x} + \lambda^2 x)y = 0\). Identify \(p(x), q(x), \lambda^2,\) and \(w(x)\).
02
Lagrange's Identity
Let \(\mathcal{L}y = (p y')' + q y\). Prove Lagrange's Identity :
\[ u\mathcal{L}v - v\mathcal{L}u = \frac{d}{dx} [p(x)(u v' - v u')] \]
03
General Proof of Orthogonality
Integrate Lagrange's identity from \(x=a\) to \(x=b\). Using the Sturm-Liouville boundary conditions (where \(p(x)(u v' - v u')|_a^b = 0\)), prove that if \(\lambda_n \neq \lambda_m\), then:
\[ (\lambda_n - \lambda_m) \int_a^b y_n y_m w(x) dx = 0 \]
Eigenfunction Mastery Exit Ticket Mastery Check
Sturm-Liouville & Special Functions
1. Identification
Identify the weight function \(w(x)\) required for the orthogonality of Legendre polynomials \(P_n(x)\) on the interval \([-1, 1]\).
2. Conceptual Symmetry
Why must we include the factor \(x\) in the integral \(\int_0^a x J_n(\alpha_i x) J_n(\alpha_j x) dx = 0\)? Connect this to the geometry of the Laplacian.
3. Completion
State the completeness property of Sturm-Liouville eigenfunctions in your own words. Why is this property the "holy grail" for mathematical physicists?
Unit Exit Ticket Sequence Final
