Example 12.33. Use convolution to solve the initial value problem

            [Graphics:Images/LaplaceConvolutionMod_gr_152.gif]  
            with
            [Graphics:Images/LaplaceConvolutionMod_gr_153.gif].  

                          [Graphics:Images/LaplaceConvolutionMod_gr_154.gif]
                                    A graph of the solution.

Explore Solution 12.33.

The solution is accomplished in two steps. (a) First solve  [Graphics:../Images/LaplaceConvolutionMod_gr_170.gif]  with  [Graphics:../Images/LaplaceConvolutionMod_gr_171.gif].  
(b)  Second solve for  [Graphics:../Images/LaplaceConvolutionMod_gr_172.gif]  with  [Graphics:../Images/LaplaceConvolutionMod_gr_173.gif].  Then the desired solution is  [Graphics:../Images/LaplaceConvolutionMod_gr_174.gif].  

(a) First solve  [Graphics:../Images/LaplaceConvolutionMod_gr_175.gif]  with  [Graphics:../Images/LaplaceConvolutionMod_gr_176.gif].  

[Graphics:../Images/LaplaceConvolutionMod_gr_177.gif]




[Graphics:../Images/LaplaceConvolutionMod_gr_178.gif]

 

 

 

 

(b)  Second solve for  [Graphics:../Images/LaplaceConvolutionMod_gr_179.gif]  with  [Graphics:../Images/LaplaceConvolutionMod_gr_180.gif].
The convolution theorem is used to make the construction of  [Graphics:../Images/LaplaceConvolutionMod_gr_181.gif]  where  [Graphics:../Images/LaplaceConvolutionMod_gr_182.gif]  and  [Graphics:../Images/LaplaceConvolutionMod_gr_183.gif].  

[Graphics:../Images/LaplaceConvolutionMod_gr_184.gif]




[Graphics:../Images/LaplaceConvolutionMod_gr_185.gif]

 

 

 

 

Aside. We can substitute v[t] into the D.E.'s and verify that it is the solution.

[Graphics:../Images/LaplaceConvolutionMod_gr_186.gif]





[Graphics:../Images/LaplaceConvolutionMod_gr_187.gif]

[Graphics:../Images/LaplaceConvolutionMod_gr_188.gif]
[Graphics:../Images/LaplaceConvolutionMod_gr_189.gif]




Remark. Mathematica 4 obtained a different form of the solution and the following comment applies when using it.
Now guide Mathematica using  trigonometric substitutions to simplify the expression for  v[t].
Notice.  This cannot be done automatically.

[Graphics:../Images/LaplaceConvolutionMod_gr_190.gif]




[Graphics:../Images/LaplaceConvolutionMod_gr_191.gif]

 

 

 

 

Aside. We can substitute y[t] into the D.E.'s and verify that it is the solution.

[Graphics:../Images/LaplaceConvolutionMod_gr_192.gif]





[Graphics:../Images/LaplaceConvolutionMod_gr_193.gif]

[Graphics:../Images/LaplaceConvolutionMod_gr_194.gif]
[Graphics:../Images/LaplaceConvolutionMod_gr_195.gif]




Aside. We use Mathematica/s DSolve package to solve the D.E. and plot the solution.

[Graphics:../Images/LaplaceConvolutionMod_gr_196.gif]





[Graphics:../Images/LaplaceConvolutionMod_gr_197.gif]



Remark. This involves the version unsimplified of v[t].  So we will plot our previously simplified function y[t].

[Graphics:../Images/LaplaceConvolutionMod_gr_198.gif]




[Graphics:../Images/LaplaceConvolutionMod_gr_199.gif]

[Graphics:../Images/LaplaceConvolutionMod_gr_200.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) 2006 John H. Mathews, Russell W. Howell