Solution of continuous-time dynamic models with inequality constraints

I propose a simple method to compute the transition process in continuous-time dynamic equilibrium models with occasionally binding inequality constraints. By augmenting the dynamic system with an auxiliary variable standard algorithms are able to solve the extended dynamic system taking the inequality into account. The procedure can handle any sequence of regimes with binding and non-binding constraints and determines the exact pattern of regimes endogenously.