During the time that it is estimating Jacobian it deliberately does not pay attention to constraints. You can use an OutputFcn to take a record of the previous boundary conditions, but as outlined above those previous results might deal with the same time or a ...