This thesis presents a solver to handle constrained-state-space ODEs. This solver locates the points where any of the states go outside of the constraining set, and then, transfers the control from the continuous time ODE to the function governing the behavior of the system on the boundaries of the constrained set. The main contribution of this solver is found simulating complex right-hand-sided ODEs for long periods of time.