How to use Lagrange Multipliers in FreeFEM

Maybe the documentation from FEniCS is cleaner? It solves exactly the same problem, so the variational formulations are the same. If you want to switch to 3D, switch mesh to mesh3 and int2d to int3d.