8.5 Series Solutions Near a Regular Singular Point, Part I We now consider the question of solving the general second order linear equation P(x)y'' + Q(x)y' + R(x)y = 0 (1) in the neighborhood of a regular singular point x = x 0 . For convenience we assume that x 0 = 0. If x 0 ≠ 0, we can transform the equation into one for which the regular singular point is at the origin by letting x − x 0 equal t. We assume that where a 0 ≠ 0. In other words, r is the exponent of the first term in the series, and a 0 is its coefficient.
30
Embed
8.5 Series Solutions Near a Regular Singular Point, Part I We now consider the question of solving the general second order linear equation P(x)y'' + Q(x)y'
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
8.5 Series Solutions Near a Regular Singular Point, Part I
We now consider the question of solving the general second order linear equation
P(x)y'' + Q(x)y' + R(x)y = 0 (1)
in the neighborhood of a regular singular point x = x0. For convenience we assume that x0 = 0. If x0 ≠ 0, we can transform the equation into one for which the regular singular point is at the origin by letting x − x0 equal t.
We assume that
where a0 ≠ 0. In other words, r is the exponent of the first term in the series, and a0 is its coefficient.
How to proceedAs part of the solution, we have to
determine:
1. The values of r for which Eq. (1) has a solution of the form (7).
2. The recurrence relation for the coefficients an.
3. The radius of convergence of the series
ExampleSolve the differential equation
2x2 y' − xy' + (1 + x)y = 0.
Answer:
1. Find the indicial equation for diff. eq.
2. Find roots of the indicial equation called the exponents at the singularity for the regular singular point x = 0.
3. Find the general solution.
8.6 Series Solutions Near a Regular Singular Point, Part II
Consider the general problem of determining a solution of the equation (1)
Where (2)
and both series converge in an interval |x| < ρ for some ρ > 0. The point x = 0 is a regular singular point, and the corresponding Cauchy–Euler equation (3) is
Ctd.We seek a solution of Eq. (1) for x > 0 and
assume that it has the form
(4)
where a0 ≠ 0, and we have written y = φ(r, x) to emphasize that φ depends on r as well as x.
We Get
where
Indicial EquationThe equation F(r ) = 0 is called the indicial
equation.The roots r1 and r2 of F(r)=0 are called the
exponents at the singularity; they determine the qualitative nature of the solution in the neighborhood of the singular point.
Setting the coefficient of xr+n in Eq. (6) equal to zero gives the recurrence relation (8)
Unequal RootsWe can always determine one solution of Eq.
(1) in the form (4)
We can also obtain a second solution
Example
Equal RootsWe can always determine one solution of Eq.
(1) in the form (4)
We can also obtain a second solution
Roots Differing by an IntegerTHEOREM 8.6.1
Consider the differential equation (1)
x2y'' + x[xp(x)]y' + [x2q(x)]y = 0,
where x = 0 is a regular singular point. Then xp(x) and x2q(x) are analytic at x = 0 with convergent power series expansions
for |x| < ρ, where ρ > 0 is the minimum of the radii of convergence of the power series for xp(x) and x2q(x). Let r1 and r2 be the roots of the indicial equation F(r ) = r (r − 1) + p0r + q0 = 0,
with r1 ≥ r2 if r1 and r2 are real.
THEOREM 8.6.1 (Ctd.)Then in either the interval −ρ < x < 0 or the interval 0 < x < ρ,
there exists a solution of the form (21)
where the an(r1) are given by the recurrence relation (8) with a0 = 1 and r = r1.
If r1 − r2 is not zero or a positive integer, then in either the interval −ρ < x < 0 or the interval 0 < x < ρ, there exists a second solution of the form (22)
The an(r2) are also determined by the recurrence relation (8) with a0 = 1 and r = r2. The power series in Eqs. (21) and (22) converge at least for |x| < ρ.
THEOREM 8.6.1 (Ctd.) If r1 = r2, then the second solution is (23)
If r1 − r2 = N, a positive integer, then (24)
The coefficients an(r1), bn(r1), cn(r2), and the constant a can be determined by substituting the form of the series solutions for y in Eq. (1). The constant a may turn out to be zero, in which case there is no logarithmic term in the solution (24). Each of the series in Eqs. (23) and (24) converges at least for |x| < ρ and defines a function that is analytic in some neighborhood of x = 0. In all three cases the two solutions y1(x) and y2(x) forma fundamental set of solutions of the given differential equation.
8.7 Bessel’s EquationWe consider three special cases of
where ν is a constant. x = 0 is a regular singular point of Eq. (1).The indicial equation is
F(r ) = r(r−1)+p0r+q0=r(r−1)+r−ν2 = r2 − ν2 = 0,with the roots r = ±ν. We will consider the three cases ν = 0, ν =
1/2, and ν = 1 for the interval x > 0.
Bessel Equation of Order ZeroIn this case ν = 0, Eq. (1) reduces toL[ y] = x2y'' + xy' + x2 y = 0, (2) and the roots of the indicial equation are equal:
r1 = r2 = 0.The first solution is
The function in brackets is known as the Bessel function of the first kind of order zero and is denoted by J0(x).
Bessel Equation of Order ZeroThe second solution of the Bessel equation
of order zero is
Instead of y2, the second solution is usually taken to be a certain linear combination of J0 and y2. It is known as the Bessel function of the second kind of order zero and is denoted by Y0.
FIGURE 8.7.1 The Bessel functions J0 and Y0.
General solution of the Bessel Equation of order 0 for x>0
Since
The general solution of the Bessel equation of order zero for x > 0 is
y = c1 J0(x) + c2Y0(x).
Bessel Equation of Order One-Half Illustrates the situation in which the roots of the
indicial equation differ by a positive integer but there is no logarithmic term in the second solution.
Setting ν = ½ in Eq. (1) gives
L[ y] = x2 y'' + xy' +(x2 − ¼)y = 0.
The roots of the indicial equation are r1 = ½, r2=− ½.
Bessel Equation of Order One-Half
One solution of the Bessel equation of order one-half is y1 = x−1/2 sin x.
The Bessel function of the first kind of order one-half, J1/2, is defined as (2/π)1/2 y1.
Bessel Equation of Order One-Half
The general solution is y = c1 J1/2(x) + c2 J−1/2(x).
where
The Bessel functions J1/2 and J−1/2.
Bessel Equation of Order OneThis case illustrates the situation in which the
roots of the indicial equation differ by a positive integer and the second solution involves a logarithmic term.
Setting ν = 1 in Eq. (1) gives
L[ y] = x2 y'' + xy' + (x2 − 1)y = 0.
Bessel Equation of Order OneWe get
The general solution for x > 0 is
y = c1 J1(x) + c2Y1 (x).
The Bessel functions J1 and Y1.
Numerical Evaluation of Bessel Functions
we have shown how to obtain infinite series solutions of Bessel’s equation of orders zero, one-half, and one.
In applications it is not unusual to require Bessel functions of other orders.
Accurate numerical approximations can be obtained by appropriate function calls in a computer algebra system in the same way that approximations of elementary functions such as sin x and ex are found.
Chapter Summary
Sections 8.2 and 8.3 Series Solutions Near an Ordinary Point, Parts I and II
Section 8.4 Regular Singular Points
Sections 8.5 and 8.6 Series Solutions Near a Regular Singular Point, Parts I and II