Add mixed Laplacian on embedded square test

Solve the mixed Laplacian with inhomogeneous Neumann boundary conditions
on a unit square embedded in 3D.