# Elementary Differential Equations and Boundary Value Problems - Boyce W.E.

**Download**(direct link)

**:**

**583**> 584 585 586 587 588 589 .. 609 >> Next

Jp (t).

For any nonnegative constant p, the differential equation

t2 X'(t) + tx (t) + (t2 — p2 )x(t) = 0

is known as Bessel’s equation of order p, and its solutions are the Bessel functions of order p. In normalized form, Bessel’s equation becomes

X'(t) + -tX(t) +

t2 - p2

x(t) = 0

From this we can see that tp(t) = 1 and t2q(t) = t2 — p2, so that tp(t) and t2q(t) are analytic at t0 = 0. Therefore zero is a regular singular point and, using equation (8), we find that the indicial equation (with P0 = 1, Q0 = — p2)

is

r(r — 1) + r — p2 = r2 — p2 = 0 Application of Frobenius’ Theorem yields a solution Jp given by the formula

Jp (t) = tpJ2

(—1)n

^ 22nn! (p + 1)( p + 2) ¦¦¦( p + n)

t2 n

n=0

The function Jp(t) is called the Bessel function of order p of the first kind. The series converges and is bounded for all t. If p is not an integer, it can be shown that a second solution of Bessel’s equation is J— p(t) and that the general solution of Bessel’s equation is a linear combination of Jp(t) and

J— p (t).

For the special case p = 0, we get the function J0 (t) used in the aging spring model in the second submodule of Module 11:

n0

©

2n t2 t

= 1 - - +

t6

4 64 2304

+

Note that even though t = 0 is a singular point of the Bessel equation of order zero, the value of J0(0) is finite [J0(0) = 1]. See Figure 11.2.

? Check that J0 (t) is a solution of Bessel’s equation of order 0.

When p is an integer we have to work much harder to get a second solution that is linearly independent of Jp(t). The result is a function Yp(t) called the Bessel function of order p of the second kind. The general formula for Yp (t) is extremely complicated. We show only the special case Y0 (t), used in the aging spring model:

Actually ó is an unending decimal and non-repeating (or so most mathematicians believe), and 0.5772 gives the first four digits.

Y0 (t) =

2

/ t \

(y + ln - +

n=0

(-1)Ë+1ßË (ë!)2

(I)

2n

where Hn = 1 + (1/2) + (1/3) + ¦¦¦ + (1/n) and y is Euler’s constant: ó =

limn

3(Hn — lnn) ú 0.5772.

210 Chapter 11

t

Figure 11.2: The graph of J0(t) [dark] looks like the graph of the decaying sinusoid *J2/ntcos(t— ÿ/4) [light].

The general solution of Bessel’s equation of integer order p is

x(t) = ci Jp (t) + C2Yp (t) (10)

for arbitrary constants c1 and c2. An important thing to note here is that the value of Yp (t) at t = 0 does reflect the singularity at t = 0; in fact, Yp (t) ^ —<x as t ^ 0+, so that a solution having the form given in equation (10) is bounded only if c = 0.

Bessel functions appear frequently in applications involving cylindrical geometry and have been extensively studied. In fact, except for the functions you studied in calculus, Bessel functions are the most widely used functions in science and engineering.

¦ Transforming Bessel's Equation to the Aging Spring Equation

See “Aging Springs” in Bessel’s equation of order zero can be transformed into the aging spring equa-

M°duie n. tion x" + e-atx = 0. To do this, we take

t = (2/a) ln(2/as) (11)

where the new independent variable s is assumed to be positive. Then we can use the chain rule to find the first two derivatives of the displacement x of the

Transforming Bessel's Equation to the Aging Spring Equation

211

We use w in place of x in the aging spring section of Module 11.

References

aging spring with respect to s:

dx dx dt dx

ds dt ds dt

as

d2x d ds2 ds

dx'

~dt

as

+

dx d

dt ds

as

d2 x dt

dx 2

dt2 ds V as/ + Ë as2 dt2 as as dt as2

d2x 4 dx 2

dt2 (as)2 dt as2

Bessel's equation of order p = 0 is given by:

2 d2 x dx 2

s2—-r + s— + Szx=0 ds2 ds

and when we substitute in the derivatives we just found, we obtain

2 ( d2 x 4 dx 2 \ dx ( 2 \ 2

I “TT 7—ã÷ + -----9 I + s~r (---I + sx = 0

dt2 (as)2 dt as2 j dt \ as)

**583**> 584 585 586 587 588 589 .. 609 >> Next