Development and validation of a P1 FEM solver for Poisson's problem (Neumann conditions): variational modeling, linear finite element discretization, calculation and assembly of mass and stiffness matrices, construction of the right-hand side, and solution of the linear system. Use of GMSH.