Checks whether the solution generated by the ODE solver contains negative values
checkSolution(soln, tol = -0.1)
soln
Matrix from ode solver out$solution
tol
tolerance