neumann, a FENICS code which uses the finite element method to solve a two dimensional boundary value problem in which homogeneous Neumann boundary conditions are imposed, based on a program by Doug Arnold.
Note that I have installed FENICS using Docker, and so to run this script I issue the commands: