Transcript Document

Ch 5.4: Regular Singular Points
Recall that the point x0 is an ordinary point of the equation
d2y
dy
P ( x) 2  Q( x)  R ( x) y  0
dx
dx
if p(x) = Q(x)/P(x) and q(x)= R(x)/P(x) are analytic at at x0.
Otherwise x0 is a singular point.
Thus, if P, Q and R are polynomials having no common
factors, then the singular points of the differential equation
are the points for which P(x) = 0.
Example 1: Bessel and Legendre Equations
Bessel Equation of order :


x 2 y  xy  x 2  2 y  0
The point x = 0 is a singular point, since P(x) = x2 is zero
there. All other points are ordinary points.
Legendre Equation:
1  x y  2xy     1y  0
2
The points x = 1 are singular points, since P(x) = 1- x2 is
zero there. All other points are ordinary points.
Solution Behavior and Singular Points
If we attempt to use the methods of the preceding two
sections to solve the differential equation in a neighborhood
of a singular point x0, we will find that these methods fail.
This is because the solution may not be analytic at x0, and
hence will not have a Taylor series expansion about x0.
Instead, we must use a more general series expansion.
A differential equation may only have a few singular points,
but solution behavior near these singular points is important.
For example, solutions often become unbounded or
experience rapid changes in magnitude near a singular point.
Also, geometric singularities in a physical problem, such as
corners or sharp edges, may lead to singular points in the
corresponding differential equation.
Numerical Methods and Singular Points
As an alternative to analytical methods, we could consider
using numerical methods, which are discussed in Chapter 8.
However, numerical methods are not well suited for the study
of solutions near singular points.
When a numerical method is used, it helps to combine with it
the analytical methods of this chapter in order to examine the
behavior of solutions near singular points.
Solution Behavior Near Singular Points
Thus without more information about Q/P and R/P in the
neighborhood of a singular point x0, it may be impossible to
describe solution behavior near x0.
It may be that there are two linearly independent solutions
that remain bounded as x  x0; or there may be only one,
with the other becoming unbounded as x  x0; or they may
both become unbounded as x  x0.
If a solution becomes unbounded, then we may want to know
if y   in the same manner as (x - x0)-1 or |x - x0|-½, or in
some other manner.
Example 1
Consider the following equation
x 2 y  2 y  0,
which has a singular point at x = 0.
It can be shown by direct substitution that the following
functions are linearly independent solutions, for x  0:
y1 ( x)  x 2 , y2 ( x)  x 1
Thus, in any interval not containing the origin, the general
solution is y(x) = c1x2 + c2 x -1.
Note that y = c1 x2 is bounded and analytic at the origin, even
though Theorem 5.3.1 is not applicable.
However, y = c2 x -1 does not have a Taylor series expansion
about x = 0, and the methods of Section 5.2 would fail here.
Example 2
Consider the following equation
x 2 y  2 xy  2 y  0
which has a singular point at x = 0.
It can be shown the two functions below are linearly
independent solutions and are analytic at x = 0:
y1 ( x)  x, y2 ( x)  x 2
Hence the general solution is
y( x)  c1 x  c2 x 2
If arbitrary initial conditions were specified at x = 0, then it
would be impossible to determine both c1 and c2.
Example 3
Consider the following equation
x 2 y  5xy  3 y  0,
which has a singular point at x = 0.
It can be shown that the following functions are linearly
independent solutions, neither of which are analytic at x = 0:
y1 ( x)  x 1 , y2 ( x)  x 3
Thus, in any interval not containing the origin, the general
solution is y(x) = c1x -1 + c2 x -3.
It follows that every solution is unbounded near the origin.
Classifying Singular Points
Our goal is to extend the method already developed for solving
P( x) y  Q( x) y  R( x) y  0
near an ordinary point so that it applies to the neighborhood of
a singular point x0.
To do so, we restrict ourselves to cases in which singularities
in Q/P and R/P at x0 are not too severe, that is, to what might
be called “weak singularities.”
It turns out that the appropriate conditions to distinguish weak
singularities are
lim x  x0 
x  x0
Q( x)
2 R( x)
is finite and lim x  x0 
is finite.
x

x
0
P( x)
P( x)
Regular Singular Points
Consider the differential equation
P( x) y  Q( x) y  R( x) y  0
If P and Q are polynomials, then a regular singular point x0
is singular point for which
lim x  x0 
x  x0
Q( x)
2 R( x)
is finite and lim x  x0 
is finite.
x

x
0
P( x)
P( x)
For more general functions than polynomials, x0 is a regular
singular point if it is a singular point with
Q( x)
x  x0 
and
P( x)
x  x0 
2
R( x)
are analytic at x  x0 .
P( x)
Any other singular point x0 is an irregular singular point.
Example 4: Bessel Equation
Consider the Bessel equation of order 


x 2 y  xy  x 2  2 y  0
The point x = 0 is a regular singular point, since both of the
following limits are finite:
Q( x)
 x 
lim x  x0 
 lim x 2   1,
x x
P ( x ) x 0  x 
2
2
2 R( x)
2  x  
2

lim x  x0 
 lim x 



x x
P( x) x0  x 2 
0
0
Example 5: Legendre Equation
Consider the Legendre equation
1  x y  2xy     1y  0
2
The point x = 1 is a regular singular point, since both of the
following limits are finite:
 2x 
 2x 
 lim 
  1,
2 
x

1
1 x 
 x 1
Q( x)



lim x  x0
 lim x  1
x x
x 1
P( x)
0
lim x  x0 
2
x  x0
R( x)
    1 
2     1 
 lim x  1 
 lim  x  1
0
2 
x

1
x

1
P( x)
 1 x 
 x 1 
Similarly, it can be shown that x = -1 is a regular singular
point.
Example 6
Consider the equation
2
2 x x  2  y  3 xy   x  2 y  0
The point x = 0 is a regular singular point:


Q( x)
3x
3x


lim x  x0 
 lim x
 lim
 0  ,
2 
2
x  x0
x

0
x

0
P( x)
2x  2
 2 x x  2  
lim x  x0 
2
x  x0
R( x)
x2 
x
2


 lim x 
 lim
0
2 
x

0
x

0
P( x)
2x  2
 2 x x  2  
The point x = 2, however, is an irregular singular point, since
the following limit does not exist:


Q( x )
3x
3x

lim x  x0 
 lim x  2

lim
2 
x  x0
P( x) x 2


2
x
x

2

 x  2 2 x x  2 
Example 7: Nonpolynomial Coefficients
(1 of 2)
Consider the equation
x   / 22 y  cos x y  sin x y  0
Note that x =  /2 is the only singular point.
We will demonstrate that x =  /2 is a regular singular point by
showing that the following functions are analytic at  /2:
cos x
cos x
x   / 2

,
2
x   / 2 x   / 2
sin x
x   / 2
 sin x
2
x   / 2
2
Example 7: Regular Singular Point
(2 of 2)
Using methods of calculus, we can show that the Taylor series of
cos x about  /2 is
(1) n 1
x   / 22n1
cos x  
n  0 ( 2n  1) !

Thus

cos x
(1) n 1
x   / 22n ,
 1  
x  / 2
n 1 ( 2n  1) !
which converges for all x, and hence is analytic at  /2.
Similarly, sin x analytic at  /2, with Taylor series
(1) n
x   / 22n
sin x  
n  0 ( 2n) !

Thus  /2 is a regular singular point of the differential equation.