Tutorial/Laplace equation with Neumann boundary condition

We propose here to solve a first very simple academic example with GMSH and GetDP. The problem is the following: \begin{equation} \begin{cases} \Delta u + u = f & \text{in } \Omega)\\ \displaystyle{\frac{\partial u}{\partial \mathbf{n}} = 0} & \text{on }\partial\Omega \end{cases} \end{equation}