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.

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

Proof. (of theorem 30.1). Suppose that<br />

is a solution of (30.70), where<br />

y = (t − t 0 ) α S(x) (30.91)<br />

S =<br />

∞∑<br />

c n (t − t 0 ) n . (30.92)<br />

n=0<br />

We will derive formulas for the c k .<br />

For n = 1, we start by differentiat<strong>in</strong>g equation (30.38)<br />

0=2(t − t 0 )S ′′ + (t − t 0 ) 2 S ′′′<br />

+ [p(t) + 2α] S ′ + (t − t 0 )p ′ (t)S ′ + (t − t 0 ) [p(t) + 2α] S ′′′<br />

+ [q ′ (t) + αp ′ (t)] S + [q(t) + α(α − 1) + αp(t)] S ′ (30.93)<br />

Everyth<strong>in</strong>g on the right hand side of this equation is analytic, and hence<br />

cont<strong>in</strong>uous. Tak<strong>in</strong>g the limit as t → t 0 , we have,<br />

0= [p(t 0 ) + 2α] S ′ (t 0 )<br />

+ [q ′ (t 0 ) + αp ′ (t 0 )] S + [q(t 0 ) + α(α − 1) + αp(t 0 )] S ′ (t 0 ) (30.94)<br />

By Taylor’s Theorem c n = S (n) (t 0 )/n! so that c 0 = S(t 0 ) and c 1 = S ′ (t 0 );<br />

similarly, p(t 0 ) = p 0 , p ′ (t 0 ) = p 1 , q(t 0 ) = q 0 , and q ′ (t 0 ) = q 1 <strong>in</strong> (30.66), and<br />

therefore<br />

0 = (p 0 + 2α)c 1 + [q 1 + αp 1 ]c 0 + [q 0 + α(α − 1) + αp 0 ]c 1 (30.95)<br />

The third term is zero (this follows from (30.68)) and hence<br />

c 1 = − q 1 + αp 1<br />

p 0 + 2α c 0. (30.96)<br />

For n > 1, we will need to differentiate equation (30.38) n times, to obta<strong>in</strong><br />

a recursion relationship between the rema<strong>in</strong><strong>in</strong>g coefficients, start<strong>in</strong>g with<br />

0 = D (n) [ (t − t 0 ) 2 S ′′] + D (n) {[2α + p(t)] (t − t 0 )S ′ }<br />

We then apply the identity<br />

+ D (n) {[α(α − 1) + αp(t) + q(t)] S} (30.97)<br />

D n (uv) =<br />

n∑<br />

k=0<br />

( n<br />

k)<br />

(D k u)(D n−k v) (30.98)

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

Saved successfully!

Ooh no, something went wrong!