We use a complementarity approach to pose the Coulomb friction model, combine it with a non-penetration condition, and append to a differential algebraic problem to characterize the dynamics of multibody systems with friction and contact. The resulting problem is relaxed to a Cone Complementarity Problem, whose solution is shown to represent the first order optimality condition of a quadratic p...