Solution 4.
![[Graphics:../Images/NumericalSolutionDEMod_gr_79.gif]](../Images/NumericalSolutionDEMod_gr_79.gif)
![[Graphics:../Images/NumericalSolutionDEMod_gr_83.gif]](../Images/NumericalSolutionDEMod_gr_83.gif)
![[Graphics:../Images/NumericalSolutionDEMod_gr_84.gif]](../Images/NumericalSolutionDEMod_gr_84.gif)
Make sure the initial condition is what we wanted.
![[Graphics:../Images/NumericalSolutionDEMod_gr_89.gif]](../Images/NumericalSolutionDEMod_gr_89.gif)
![[Graphics:../Images/NumericalSolutionDEMod_gr_90.gif]](../Images/NumericalSolutionDEMod_gr_90.gif)
Unfortunately, the particular choice of quotients of Bessel functions introduces a problem at x=0.
![[Graphics:../Images/NumericalSolutionDEMod_gr_92.gif]](../Images/NumericalSolutionDEMod_gr_92.gif)
This is one of the mysteries of differential equations.
Try sneaking up on x = 0.
![[Graphics:../Images/NumericalSolutionDEMod_gr_95.gif]](../Images/NumericalSolutionDEMod_gr_95.gif)
![[Graphics:../Images/NumericalSolutionDEMod_gr_97.gif]](../Images/NumericalSolutionDEMod_gr_97.gif)
Or we can take the limit.
![[Graphics:../Images/NumericalSolutionDEMod_gr_99.gif]](../Images/NumericalSolutionDEMod_gr_99.gif)