On effectue la même procédure de Galerkin qu'en I.2.1.1 p. . On a donc le même théorème d'existence et d'unicité du problème variationnel dans Vh (I.2.1) et l'on définit toujours les matrices du problème linéaire (I.5.36).
La matrice D est la matrice du produit scalaire dans Vh, C est la matrice de la forme bilinéaire
, le second membre, par abus de langage, est toujours noté b.
Nous présentons comment calculer les termes du système linéaire I.5.36 puis comment calculer une approchée discrète de la solution du problème I.5.21. L'implémentation numérique est proposée par le choix d'un espace Vh qui permet une inversibilité inconditionnelle du
système I.5.36.