Issues with Convergence in Differential Equation Solver Using Dual Numbers

Continuation of this post here: Resolving Julia Warning: dt(2.27e-13) <= dtmin(2.27e-13) in implementation of dual numbers