Numerical Solution Of Mixed Linear Complementarity Problems In Multibody Dynamics With Contact