Sturm--Liouville theory is actually a generalization for infinite dimensional case the famous eigenvalue/eigenvector problems for finite square matrices that we discussed in Part I of this tutorial. Although a Sturm--Liouville problem can be formulated in operator form as L[ y ] = λy similarly to matrix eigenvalue problem Ax = λx, the operator L is usually an unbounded differential operator and y is a smooth function.
The corresponding theory that is known as Sturm--Liouville theory originated in two articles by the French mathematician of German decent Jacques Charles François Sturm (1803--1855).
Sturm, J.C.F.,
"Mémoire sur les Équations différentielles linéaires du second ordre",
Journal de Mathématiques Pures et Appliquées, 1836, Vol. 1,
pp. 106-186. [sept. 28, 1833]
Sturm, , J.C.F.,
"Mémoire sur une classe des d'Équations à différences partielles".
Journal de Mathématiques Pures et Appliquées, 1836, Vol. 1,
pp. 373--444.
Following the same 1836 year, Sturm together with his friend, the
French
mathematician Joseph
Liouville (1809--1882), published very influential articles that provided crucial groundwork for the theory.
Sturm & Liouville,
"Démonstration d'un théorème de M. Cauchy relatif aux racines imaginaires des équations".
Journal de Mathématiques Pures et Appliquées, 1836, Vol. 1,
pp. 278--289.
Sturm & Liouville,
"Note sur un théorème de M. Cauchy relatif aux racines des équations simultanées", Comptes rendus de l'Académie des Sciences (English: Proceedings of the Academy of Sciences), 1837, Vol. 4, pp. 720--739.
Sturm & Liouville,
"Extrait d'un Mémoire sur le développement des fonctions en séries dont les différents termes sont assujettis à satisfaire à une même équation différentielle lindaire, contenant un paramètre variable",
Journal de Mathématiques Pures et Appliquées, 1837, Vol 2,
pp. 220--233.
Comptes rendus de l'Académie des Sciences ((English: Proceedings of the Academy of sciences), 1837, Vol. 4, pp. 675--677.
Return to computing page for the first course APMA0330
Return to computing page for the second course APMA0340
Return to Mathematica tutorial for the first course APMA0330
Return to Mathematica tutorial for the second course APMA0340
Return to the main page for the first course APMA0330
Return to the main page for the second course APMA0340
Return to Part V of the course APMA0340
Introduction to Linear Algebra with Mathematica
Here y is a nontrivial (meaning not identically zero) function of the free variable x, and the functions p(x), q(x), and w(x) > 0 are specified at the outset. In the simplest of cases all coefficients are continuous on the finite closed interval [a,b], and p has continuous derivative. In this simplest of all cases, this nontrivial function y is called a solution if it is continuously differentiable on (a,b) and satisfies the given equation at every point in (a,b). The function w(x) is called the "weight" or "density" function.
\[
L\left[ x, \texttt{D} \right] y = w\,\lambda y , \qquad B_0 y = 0, \quad B_{\ell} y = 0.
\]
This boundary value problem has an obvious solution---identically zero function. Since we are not after such a trivial solution, we need something more. The Sturm--Liouville problem (S-L, for short) consists of two parts: the first part is about finding such values of parameter λ for which the problem has a nontrivial solution (not identically zero); such values are called eigenvalues. The second part includes determination of nontrivial solutions that are called eigenfunctions. Note that a Sturm--Liouville problem may have other constraints, not necessarily formulated above.
There are two kinds of Sturm--Liouville problems. One is called classical or regular when p(x) > 0 and w(x) > 0 for all points from the closed interval x∈[0,ℓ]. These assumptions are necessary to render the theory as simple as possible while retaining considerable generality. It turns out that these conditions are valid in many problems that we will consider shortly. If p(0) = 0, or p(ℓ) = 0, or p(0) = 0 = p(ℓ), the Sturm--Liouville problem is said to be singular.
Jacques Charles François Sturm (1803--1855) was a French mathematician of Switzerland descend. Charles spent his adult life in Paris. His primary interests were fluid mechanics and differential equations. Sturm with the Swiss engineer Daniel Colladon was the first to accurately determine the speed of sound in water. In mathematics, he won the coveted Grand prix des Scienes Mathematiques for his work in differential equations. Sturm held the chair of mechanics at the Sorbonne and was elected a member of the French Academy of Sciences. The asteroid 31043 Sturm is named for him. Sturm's name is one of the 72 names engraved at the Eiffel Tower.
Joseph Liouville was a French mathematician known for his work in analysis, differential geometry, and number theory and for his discovery of transcendental numbers---i.e., numbers that are not the roots of algebraic equations having rational coefficients. He was also influential as a journal editor and teacher. Joseph founded the Journal de Mathématiques Pures et Appliquées which retains its high reputation up to today.
In physics and other applications, many problems arise in the form of boundary
value problems involving second order ordinary differential equations, written
in self-adjoint form:
where \( \texttt{D} = {\text d}/{\text d}x \) is
the derivative operator. Sturm was the first who generalized the well known
matrix problem for eigenvalues/eigenvectors to the linear differential
operators by considering the eigenvalue problem:
\[
L\left[ x, \texttt{D} \right] y + \lambda\,\rho (x)\, y =0
\]
for functions y(x) subject to some boundary conditions. Now this problem is referred to as the Sturm--Liouville problem. Here p(x), q(x), and ρ(x) > 0 are specified at the outset, which is usually is some interval (finite or not) of real axis ℝ.
Theorem: Properties of the regular Sturm--Liouville problem.
There exists an infinite number of real eigenvalues that can be arranged in increasing order λ_{1} < λ_{2} < λ_{3} < … λ_{n} < … such that λ_{n} →∞ as n→∞.
For each eigenvalue there is only one eigenfunction (up until nonzero
multiple).
Eigenfunctions corresponding to different eigenvalues are linearly
independent.
The set of eigenfunctions { y_{n} }_{n≥0} corresponding to the set of eigenvalues is
orthogonal with respect to the weight function w(x) on the interval x∈
[0, ℓ]:
The Fourier expansion of a square integrable generic function f ∈ 𝔏²[𝑎,b] on the
interval (𝑎 , b) in terms of an orthogonal system of
real-valued functions
{ φ_{k} } in 𝔏² is given by
\[
f(x) = \sum_{k\ge 1} c_k \varphi_k (x) ,
\]
where the coefficients are given by the inner product
The term complete was introduced by the famous Russian mathematician
Vladimir Steklov (1864--1926). The above identity was stated by the famous French mathematician Marc-Antoine Parseval (1755--1836) in 1799.
The main reason to study Sturn--Liouville problems is that their eigenfunctions
provide a basis for expansion for certain class of functions.
This tells us that basis functions have close relationships with
linear operators, and this basis is orthogonal for self-adjoint second order
differential operators. In other words, solutions of differential equations
can be approximated more efficiently by means of better basis functions. For
example, a periodic function or solution is expressed more efficiently
by periodic basis functions than by polynomials like the Taylor series.
For example, the linear operator
\( L\left[ y \right] = y'' + \lambda \,y \)
determines the basis functions. When λ < 0, the basis functions
will be
is more suitable for approximation of periodic functions than Taylor's series.
Example 1:
Consider the Sturm--Liouville problem
\[
y'' + \lambda \,y =0 , \qquad y(0) = y(T) .
\]
Solutions of the differential equation \( y'' + \lambda \,y =0 \) depend on the sign of parameter λ. If λ = -μ^{2} is negative, the equation has two exponential linearly independent solutions \( y_1 = e^{\mu x} \quad\mbox{and}\quad y_2 = e^{-\mu x} \) that also sometimes can be written through hyperbolic sine and cosine functions. These functions cannot be periodic, so negative λ is not an eigenvalue.
If λ = 0, the general solution of the differential equation \( y'' =0 \) is a linear function \( y = a + b\,x , \) which could be periodic only when b = 0. So a zero is an eigenvalue to which corresponds an eigenfunction---a constant function.
If λ = μ^{2} is positive, then the general solution becomes
\[
y = a\,\cos (\mu x) + b\,\sin (\mu x) ,
\]
with some constants 𝑎 and b. This function is periodic with period T if and only if μ is proportional to 2π/T; so we get eigenvalues
When a function is specified at the boundary, then the corresponding boundary conditions are called after G. Dirichlet.
Solutions of the differential equation \( y'' + \lambda \,y =0 \) depend on the sign of parameter λ. If λ = -μ^{2} is negative, the equation has two exponential linearly independent solutions \( y_1 = e^{\mu x} \quad\mbox{and}\quad y_2 = e^{-\mu x} . \)
Then the general solution becomes
\[
y(x) = c_1 e^{\mu x} + c_2 e^{-\mu x} ,
\]
with some constants c_{1} and c_{2}. From the first boundary condition, we have
It could be shown similarly to the previous example that negative λ do not contribute. So we need to consider only two cases: λ = 0 and λ = μ² > 0. The former gives us an eigenvalue λ = 0 to which corresponds a constant as an eigenfunction. To the latter corresponds the general solution
Example 4:
There are known two Sturm--Liouville problems with mixed boundary conditions when on one end we have the Dirichlet condition while on the other end we have the Neumann condition. So we start with one of them:
with some constant c. It is clear from our previous discussion that eigenvalues must be positive in this case. The Neumann condition at the right end dictates
When a product is zero then at least one multiple must be zero. Obviously, c cannot be zero, otherwise we will have a trivial solution. Since λ > 0 according to assumption, we get the condition
Choosing c_{2} ≠ 0, we see that the last equation is equivalent to
\[
\tan \sqrt{\lambda} = - \sqrt{\lambda} .
\]
The above transcendent equation could be solved only by a numerical procedure such as Newton's method. Nevertheless, if we let \( t = \sqrt{\lambda} , \) then we see from their graphs that there exists an infinite number of roots of \( \tan t = -t . \)
Since tangent function has vertical asymptotes at \( t = \frac{\pi}{2} + n\pi , \quad n=0,1,2,\ldots ; \) the roots t_{n} of the equation\( \tan t = -t \) approach
Return to Mathematica page)
Return to the main page (APMA0340)
Return to the Part 1 Matrix Algebra
Return to the Part 2 Linear Systems of Ordinary Differential Equations
Return to the Part 3 Non-linear Systems of Ordinary Differential Equations
Return to the Part 4 Numerical Methods
Return to the Part 5 Fourier Series
Return to the Part 6 Partial Differential Equations
Return to the Part 7 Special Functions