Example 3.  Consider the wave equation where  [Graphics:Images/FiniteDifferencePDEMod_gr_78.gif].  The string at rest has length  [Graphics:Images/FiniteDifferencePDEMod_gr_79.gif].  Assume that the initial position

        [Graphics:Images/FiniteDifferencePDEMod_gr_80.gif]  

Use the finite difference method to solve the wave equation over the rectangle  [Graphics:Images/FiniteDifferencePDEMod_gr_81.gif].  

Solution 3.

Remark.  We will be using  [Graphics:../Images/FiniteDifferencePDEMod_gr_82.gif] .  This forces  [Graphics:../Images/FiniteDifferencePDEMod_gr_83.gif].

Using the same functions as in Example 2, we increase the interval width ever so slightly to be  [Graphics:../Images/FiniteDifferencePDEMod_gr_84.gif].  Notice the instability this produces which was predicted in advance in our mathematical analysis discussed in the text.

 

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


Now set up the table of solutions.  Solve it.  Plot it.  

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


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

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

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

[Graphics:../Images/FiniteDifferencePDEMod_gr_90.gif]
[Graphics:../Images/FiniteDifferencePDEMod_gr_91.gif]
[Graphics:../Images/FiniteDifferencePDEMod_gr_92.gif]
[Graphics:../Images/FiniteDifferencePDEMod_gr_93.gif]

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004