Wie schon in Abschnitt 3.3 gezeigt wurde, besitzt
jedes Element mit den lokalen Knotennummern
eine fixe Abbildung auf die globalen Gitterknoten.
Jedes Element besitzt einen Abbildungsvektor der Länge
, der einen lokalen
Index
der Elementsmatrix auf einen globalen Index
der Gesamtmatrix
abbildet und
des weiteren als Transformation
bezeichnet wird.
Summiert man alle
Elementsintegrale zu einem Gesamtintegral

auf, geht man von den lokalen Elementsknoten zu einer globalen Knotennumerierung des
Gesamtgitters mit
Knoten

über und führt man die Dirichletschen Randbedingungen derart ein, daß für die
Dirichlet-Knoten
der Bereich von
bis
reserviert ist, so erhält man das Gesamtintegral
Da die Gesamtmatrix symmetrisch ist,
(siehe auch Abschnitt 3.10.1),
vereinfacht sich (3.48) zu
Der letzte Term in (3.49) ist unabhängig von den gesuchten
Potentialen
für
, und wird daher für die weiteren Berechnungen als
Konstante
bezeichnet. Die ebenfalls konstante Spaltenmatrix
besteht aus
Dirichlet-Knotenwerten, multipliziert mit Matrixkoeffizienten
.
Bildet man die partiellen Ableitungen nach den unbekannten Knotenpotentialen
für
von (3.49) und setzt das Gleichungssystem null,
so hat man das gesuchte Minimum gefunden.
Die partiellen Ableitungen nach der Komponente
und
ergeben

bzw.
