Differential Equations Project

Computer Lab Modules

 

Solution of D.E.'s using Laplace Transforms

 

 

Background. Laplace transforms are useful in solving initial value problems for ordinary differential equations and systems of ordinary differential equations. In this exercise we explore the method of using Laplace transforms. First load Mathematica's built in "LaplaceTransform" subroutine package.

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr1.gif]
 

Computer Lab Work.

 

Exercise 1. Use Laplace transforms to solve the initial value problem
[Graphics:e18.txtgr3.gif]
and plot the solution over the interval 0 <= t <= 7.

First set up the initial conditions and function f(t).

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr4.gif]

Second, set up the D.E. in Laplace transformation format and solve for Y[s].

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr5.gif]

Since Y[s] is in its expanded form, we can find the inverse transform of each term in the above sum.
We could look them up in a table or let the following loop do it for us.

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr6.gif]

Find the inverse transform and plot the solution.

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr7.gif]


Exercise 2. Use Laplace transforms to solve the initial value problem
[Graphics:e18.txtgr8.gif]
and plot the solution over the interval 0 <= t <= 7.

First set up the initial conditions and function f(t).

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr9.gif]

Second, set up the D.E. in Laplace transformation format and solve for Y[s].

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr10.gif]

Find the inverse transform and plot the solution.

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr11.gif]


Exercise 3. Use Laplace transforms to solve the initial value problem for the system of D.E.'s
[Graphics:e18.txtgr12.gif]
and plot the solution for 0 <= t <= 5.5

First set up the initial conditions and functions f(t) and g(t)

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr13.gif]

Second, set up the D.E. in Laplace transformation format and solve for X[s] and Y[s].

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr14.gif]

Find the inverse transforms and plot the solution.

[Graphics:e18.txtgr2.gif][Graphics:e18.txtgr15.gif]
 

Solutions.

 

 

Return to the Differential Equations Project

Return to the Numerical Analysis Project

Return to the Complex Analysis Project

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews, 1998