Checks whether the solution generated by the ODE solver contains negative values
checkSolution(soln, tol = -0.1)
soln
tol