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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

252 LESSON 27. HIGHER ORDER EQUATIONS<br />

Variation of Parameters<br />

Theorem 27.12 (Variation of Parameters.). Suppose that y 1 , ..., y n are a<br />

fundamental set of solutions to<br />

L n y = a n (t)y (n) + a n−1 (t)y (n−1) + · · · + a 0 (t)y = 0 (27.194)<br />

Then a particular solution to<br />

is given by<br />

L n y = a n (t)y (n) + a n−1 (t)y (n−1) + · · · + a 0 (t)y = f(t) (27.195)<br />

∫<br />

∫<br />

W 1 (s)f(s)<br />

y P = y 1<br />

t W (s)a n (s) ds + y W 2 (s)f(s)<br />

2<br />

t W (s)a n (s) ds<br />

∫<br />

W n (s)f(s)<br />

+ · · · + y n ds (27.196)<br />

W (s)a n (s)<br />

t<br />

t<br />

where W j (t) is the determ<strong>in</strong>ant of W[y 1 , ..., y n ](t) with the jth column<br />

replaced by a vector with all zeroes except for a 1 <strong>in</strong> the last row. In<br />

particular, for n = 2, a particular solution to a(t)y ′′ + b(t)y ′ + c(t)y = f(t)<br />

is<br />

∫<br />

∫<br />

y 2 (s)f(s)<br />

y p = −y 1 (t)<br />

W (s)a(s) ds + y y 1 (s)f(s)<br />

2(t)<br />

ds (27.197)<br />

W (s)a(s)<br />

Proof for general case. Look for a solution of the form<br />

t<br />

y = u 1 y 1 + · · · + u n y n (27.198)<br />

This is under-determ<strong>in</strong>ed so we can make additional assumptions; <strong>in</strong> particular,<br />

we are free to assume that<br />

Then<br />

u ′ 1y 1 + · · · + u ′ ny n = 0<br />

u ′ 1y ′ 1 + · · · + u ′ ny ′ n = 0<br />

.<br />

u ′ 1y (n−2)<br />

1 + · · · + u ′ ny n (n−2) = 0<br />

(27.199)<br />

y ′ = u 1 y 1 ′ + · · · + u n y n<br />

′<br />

y ′′ = u 1 y 1 ′′ + · · · + u n y n<br />

′′<br />

.<br />

y (n−1) = u 1 y (n−1)<br />

1 + · · · + u n y n<br />

(n−1)<br />

y (n) = u 1 y (n)<br />

1 + · · · + u n y n<br />

(n)<br />

+ u ′ 1y (n−1)<br />

1 + · · · + u ′ ny n<br />

(n−1)<br />

(27.200)

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

Saved successfully!

Ooh no, something went wrong!