How to use Lagrange Multipliers in FreeFEM

f is the source term, vb discretizes the constraint int(u) = b1.