Solve an obstacle problem (i.e, a partial differential equation with constraints) - View it on GitHub