The above solution is really only appropriate if the reason that c(1) becomes negative is due to round-off or overshooting (that is, a smaller step size would not have gone negative.)
If you are trying to do the equivalent of bounding off a floor, then it is not appropriate to use this solution, and instead you should use an ode event function to terminate the integration, then change the boundary condition, and then restart the integration.
