You can avoid the numerical glitch at about x=7.7 seen in the pic either by setting TOL to 10^-6 or lower or by adding 10^5 as the third argument (intervals) of the odesolve function
You can avoid the numerical glitch at about x=7.7 seen in the pic either by setting TOL to 10^-6 or lower or by adding 10^5 as the third argument (intervals) of the odesolve function