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.

190 LESSON 22. NON-HOMOG. EQS.W/ CONST. COEFF.<br />

Algorithm for 2nd Order, L<strong>in</strong>ear, Constant Coefficients IVP<br />

To solve the <strong>in</strong>itial value problem<br />

ay ′′ + by ′ + cy = f(t)<br />

y(t 0 ) = y 0<br />

y ′ (t 0 ) = y 1<br />

⎫<br />

⎪⎬<br />

⎪ ⎭<br />

(22.21)<br />

1. F<strong>in</strong>d the roots of the characteristic polynomial ar 2 + br + c = 0.<br />

2. Factor the differential equation as a(D − r 1 ) (D − r 2 )y = f(t).<br />

} {{ }<br />

z(t)<br />

3. Substitute z = (D − r 2 )y = y ′ − r 2 y to get a(D − r 1 )z = f(t).<br />

4. Solve the result<strong>in</strong>g differential equation z ′ − r 1 z = 1 f(t) for z(t).<br />

a<br />

5. Solve the first order differential equation y ′ − r 2 y = z(t) where z(t) is<br />

the solution you found <strong>in</strong> step (4).<br />

6. Solve for arbitrary constants us<strong>in</strong>g the <strong>in</strong>itial conditions.<br />

Example 22.1. Solve the <strong>in</strong>itial value problem<br />

⎫<br />

y ′′ − 10y ′ + 21y = 3 s<strong>in</strong> t⎪⎬<br />

y(0) = 1<br />

⎪⎭<br />

y ′ (0) = 0<br />

The characteristic polynomial is<br />

(22.22)<br />

r 2 − 10r + 21 = (r − 3)(r − 7) (22.23)<br />

Thus the roots are r = 3, 7 and s<strong>in</strong>ce a = 1 (the coefficient of y ′′ ), the<br />

differential equation can be factored as<br />

(D − 3) (D − 7)y = 3 s<strong>in</strong> t (22.24)<br />

} {{ }<br />

z<br />

To solve equation (22.24) we make the substitution<br />

Then (22.24) becomes<br />

z = (D − 7)y = y ′ − 7y (22.25)<br />

(D − 3)z = 3 s<strong>in</strong> t (22.26)<br />

z ′ − 3z = 3 s<strong>in</strong> t (22.27)

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

Saved successfully!

Ooh no, something went wrong!