La meilleurs façon de s'assurer que le programme donne le bon résultat est de partir d'une fonction connue, par exemple u = x*x + y*y et de choisir les données f,u0,g, pour retrouver cette solution. Ainsi avec a=nu=1, f = x*x+y*y -4, u0=u et sur un cercle intérieur de rayon r, g= grad u .n = (2x,2y).(-x/r, -y/r) = -2.
Faire les simulations du chapitre 1 du livre Lucquin-Pironneau[] avec ce module
quand l'EDP est du bon type.