Solving Linear ODE Using Laplace Transforms

How can we use Laplace transforms to solve ode? The procedure is best
illustrated with an example. Consider the ode

displaymath61

This is a linear homogeneous ode and can be solved using standard methods.

Let Y(s)=L[y(t)](s). Instead of solving directly for y(t), we derive a
new equation for Y(s). Once we find Y(s), we inverse transform
to determine y(t).

The first step is to take the Laplace transform of both sides of the
original differential equation. We have

displaymath63

Obviously, the Laplace transform of the function 0 is 0. If we look at
the left-hand side, we have

displaymath65

Now use the formulas for the L[y'']and L[y']:

displaymath67

Here we have used the fact that y(0)=2. And,

displaymath69

Hence, we have

displaymath71

The Laplace-transformed differential equation is

displaymath73

This is a linear algebraic equation for Y(s)! We have converted a
differential equation into a algebraic equation! Solving for Y(s), we have

displaymath75

We can simplify this expression using the method of partial fractions:

displaymath77

Recall the inverse transforms:

displaymath79

Using linearity of the inverse transform, we have

displaymath81

Another Example

Consider the ode:

displaymath83

This problem has an inhomogeneous term. In the direct approach
one solves for the homogeneous solution and the particular solution
separately. For this problem the particular solution can be determined
using variation of parameters or the method of undetermined
coefficients. Using the Laplace transform technique we can solve for
the homogeneous and particular solutions at the same time.

Let Y(s) be the Laplace transform of y(t). Taking the Laplace transform
of the differential equation we have:

displaymath85

The Laplace transform of the LHS  L[y''+4y'+5y] is

displaymath87

The Laplace transform of the RHS is

displaymath89

Equating the LHS and RHS and using the fact that y(0)=1 y'(0)=2, we
obtain

displaymath91

Solving for Y(s), we obtain:

displaymath93

Using the method of partial fractions it can be shown that

displaymath95

Using the fact that the inverse of 1/(s-1) is e^t and that the inverse of
1/[(s+2)^2+1] is exp(2t)sin(t), it follows that

displaymath97


[ODE Home] [1st-Order Home] [2nd-Order Home] [Laplace Transform Home] [Notation] [References]

Copyright © 1996 Department of Mathematics, Oregon State University

If you have questions or comments, don't hestitate to contact us.