17.02.2014 Views

Mathematical Methods for Physicists: A concise introduction - Site Map

Mathematical Methods for Physicists: A concise introduction - Site Map

Mathematical Methods for Physicists: A concise introduction - Site Map

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

ORDINARY DIFFERENTIAL EQUATIONS<br />

Case 1 Distinct roots not di€ering by an integer<br />

This is the simplest case. Let r 1 and r 2 be the roots of the indicial equation (2.29).<br />

If we insert r ˆ r 1 into the recurrence relation and determine the coecients<br />

a 1 , a 2 ; ... successively, as be<strong>for</strong>e, then we obtain a solution<br />

y 1 …x† ˆx r 1<br />

…a 0 ‡ a 1 x ‡ a 2 x 2 ‡†:<br />

Similarly, by inserting the second root r ˆ r 2 into the recurrence relation, we will<br />

obtain a second solution<br />

y 2 …x† ˆx r 2<br />

…a 0<br />

* ‡ a 1<br />

*x ‡ a 2<br />

*x 2 ‡†:<br />

Linear independence of y 1 and y 2 follows from the fact that y 1 =y 2 is not constant<br />

because r 1 r 2 is not an integer.<br />

Case 2 Double roots<br />

The indicial equation (2.29) has a double root r if, and only if,<br />

…g 0 1† 2 4h 0 ˆ 0, and then r ˆ…1 g 0 †=2. We may determine a ®rst solution<br />

<br />

y 1 …x† ˆx r …a 0 ‡ a 1 x ‡ a 2 x 2 ‡† r ˆ 1g <br />

0<br />

…2:30†<br />

2<br />

as be<strong>for</strong>e. To ®nd another solution we may apply the method of variation of<br />

parameters, that is, we replace constant c in the solution cy 1 …x† by a function<br />

u…x† to be determined, such that<br />

y 2 …x† ˆu…x†y 1 …x†<br />

…2:31†<br />

is a solution of Eq. (2.27). Inserting y 2 and the derivatives<br />

y2 0 ˆ u 0 y 1 ‡ uy1 0 y2 00 ˆ u 00 y 1 ‡ 2u 0 y1 0 ‡ uy1<br />

00<br />

into the di€erential equation (2.27) we obtain<br />

x 2 …u 00 y 1 ‡ 2u 0 y1 0 ‡ uy1 00 †‡xg…u 0 y 1 ‡ uy1†‡huy 0 1 ˆ 0<br />

or<br />

x 2 y 1 u 00 ‡ 2x 2 y1u 0 0 ‡ xgy 1 u 0 ‡…x 2 y1 00 ‡ xgy1 0 ‡ hy 1 †u ˆ 0:<br />

Since y 1 is a solution of Eq. (2.27), the quantity inside the bracket vanishes; and<br />

the last equation reduces to<br />

x 2 y 1 u 00 ‡ 2x 2 y1u 0 0 ‡ xgy 1 u 0 ˆ 0:<br />

Dividing by x 2 y 1 and inserting the power series <strong>for</strong> g we obtain<br />

<br />

u 00 ‡ 2 y 1<br />

0 ‡ g <br />

0<br />

y 1 x ‡ u 0 ˆ 0:<br />

Here and in the following the dots designate terms which are constants or involve<br />

positive powers of x. Now from Eq. (2.30) it follows that<br />

y1<br />

0 ˆ xr1 ‰ra 0 ‡…r ‡ 1†a 1 x ‡Š<br />

y 1 x r ˆ 1 ra 0 ‡…r ‡ 1†a 1 x ‡<br />

ˆ r<br />

‰a 0 ‡ a 1 x ‡Š x a 0 ‡ a 1 x ‡ x ‡:<br />

90

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!