Example 2.  Consider the heat equation where  [Graphics:Images/CrankNicolsonMod_gr_43.gif].   The length of the rod is  [Graphics:Images/CrankNicolsonMod_gr_44.gif].  Assume that the ends of the rod are held at the temperature  [Graphics:Images/CrankNicolsonMod_gr_45.gif].  Assume that the initial temperature distribution is

      [Graphics:Images/CrankNicolsonMod_gr_46.gif].
      
Now investigate what happens when the step size is too large.  This time the t interval is larger  [0.0, 0.35], and the step size is "too large."
Apply the forward difference method with  [Graphics:Images/CrankNicolsonMod_gr_47.gif]  and  obtain temperature distributions for  [Graphics:Images/CrankNicolsonMod_gr_48.gif].
We will use  [Graphics:Images/CrankNicolsonMod_gr_49.gif] .  This forces  [Graphics:Images/CrankNicolsonMod_gr_50.gif].

Solution 2.

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


Now set up the table of solutions.

 

 

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

Next, solve it.

[Graphics:../Images/CrankNicolsonMod_gr_53.gif]
[Graphics:../Images/CrankNicolsonMod_gr_54.gif]


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

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

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

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

[Graphics:../Images/CrankNicolsonMod_gr_59.gif]
[Graphics:../Images/CrankNicolsonMod_gr_60.gif]
[Graphics:../Images/CrankNicolsonMod_gr_61.gif]
[Graphics:../Images/CrankNicolsonMod_gr_62.gif]


To see the numerical values enter the command:

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




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

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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004