By Evalf, and other Nutils contributors
Solves Poisson's equation
Δu = 1 subject to boundary constraints, using
the fact that the solution to the strong form minimizes the functional
∫ .5 ‖∇u‖² - u. The domain is a unit square, and the solution is constrained
to zero along the entire boundary.