Equazione del potenziale (D-AN3)

Risolvere con un metodo alle differenze (ad esempio lo schema a 5 punti) il problema di Neumann

$$\cases{-\Delta u+u = f & in $\Omega$, \cr {\partial\over\partial n}u = 0 & su $\partial\Omega$}$$

dove $\Omega=[0,1]$ e $\Omega=[0,1]\times[0,1]$. Rappresentare graficamente la soluzione.