A MATLAB implementation of a 2D particle physics simulation based on solving MLCP problem using Projected Gauss-Seidel method.