next up previous
Next: About this document

Solving Linear ODE Using the Laplace Transform Method

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''](s) and L[y']:

displaymath67

Here we have used the fact that y(0)=1

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 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). Hence, the solution is

displaymath97




next up previous
Next: About this document


Tue Jul 2 09:31:44 PDT 1996