21.04.2015 Views

Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro

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.

296 LESSON 30. THE METHOD OF FROBENIUS<br />

By the <strong>in</strong>dicial equation, α 2 + α(p 0 − 1) + q 0 = 0, so that<br />

α = 1 2<br />

[1 − p 0 ± √ (1 − p 0 ) 2 − 4q 0<br />

]<br />

(30.126)<br />

Designate the two roots by α + and α − and def<strong>in</strong>e<br />

∆ = α + − α − = √ (1 − p 0 ) 2 − 4q 0 (30.127)<br />

If they are real, then the larger root satisfies 2α + ≥ 1 − p 0 , so that<br />

n − 1 + 2α + + p 0 ≥ n − 1 + 1 − p 0 + p 0 = n > 0 (30.128)<br />

Therefore (30.124) gives a solution for the larger of the two roots.<br />

other root gives<br />

The<br />

n − 1 + 2α − + p 0 = n − 1 + 1 − p 0 − √ (p 0 − 1) 2 − 4q 0 + p 0<br />

(30.129)<br />

This is never zero unless the two roots differ by an <strong>in</strong>teger. Thus the smaller<br />

root also gives a solution, so long as the roots are different and do not differ<br />

by an <strong>in</strong>teger.<br />

If the roots are complex, then<br />

where ∆ ≠ 0 and therefore<br />

α = 1 2 [1 − p 0 ± i∆] (30.130)<br />

n − 1 + 2α + p 0 = n − 1 + 1 − p 0 ± i∆ + p 0 = n ± i∆ (30.131)<br />

which can never be zero. Hence both complex roots lead to solutions.<br />

When the difference between the two roots of the <strong>in</strong>dicial equation is an<br />

<strong>in</strong>teger, a second solution can be found by reduction of order.<br />

Theorem 30.2. (Second Frobenius solution.) Suppose that p(t) and q(t)<br />

are both analytic with some radius of convergence r, and that α 1 ≥ α 2 are<br />

two (possibly dist<strong>in</strong>ct) real roots of the <strong>in</strong>dicial equation<br />

α 2 + (p 0 − 1)α + q 0 = 0 (30.132)<br />

and let ∆ = α 1 − α 2 . Then for some set of constants {c 0 , c 1 , ...}<br />

y 1 (t) = (t − t 0 ) α1<br />

∞ ∑<br />

k=0<br />

c k (t − t 0 ) k (30.133)

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

Saved successfully!

Ooh no, something went wrong!