Solution 4.

[Graphics:../Images/NumericalSolutionDEMod_gr_79.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_80.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_81.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_82.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_83.gif]

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

[Graphics:../Images/NumericalSolutionDEMod_gr_85.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_86.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_87.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_88.gif]

Make sure the initial condition is what we wanted.

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

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

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

Unfortunately, the particular choice of quotients of Bessel functions introduces a problem at x=0.

[Graphics:../Images/NumericalSolutionDEMod_gr_92.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_93.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_94.gif]

This is one of the mysteries of differential equations.
Try sneaking up on x = 0.

[Graphics:../Images/NumericalSolutionDEMod_gr_95.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_96.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_97.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_98.gif]

Or we can take the limit.

[Graphics:../Images/NumericalSolutionDEMod_gr_99.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_100.gif]
[Graphics:../Images/NumericalSolutionDEMod_gr_101.gif]


Converted by Mathematica      January 28, 2003