A square sliding/residing on a slope under gravity is simulated with mass-spring elasticity potential and implicit Euler time integration. Each time step is solved by minimizing the Incremental Potential with the projected Newton method.
pip install numpy scipy pygame
python simulator.py