We propose a method to compute equilibria in dynamic models with several continuous state variables and occasionally binding constraints. These constraints induce non-differentiabilities in policy functions. We develop an interpolation technique that addresses this problem directly: It locates the non-differentiabilities and adds interpolation nodes there. To handle this flexible grid, it uses Delaunay interpolation, a simplicial interpolation technique. Hence, we call this method Adaptive Simplicial Interpolation (ASI). We embed ASI into a time iteration algorithm to compute recursive equilibria in an infinite horizon endowment economy where heterogeneous agents trade in a bond and a stock subject to various trading constraints. We show that this method computes equilibria accurately and outperforms other grid schemes by far.