spezielle Matrix, die aus dem Finite-Elemente-Ansatz einer parabolischen Differentialgleichung u_t − ∆u = f entsteht.

Sind die φ_k die Basisfunktionen des Finite- Elemente-Raums, so entsteht durch den Lösungsansatz \begin{eqnarray}u(t):=\displaystyle \sum {y}_{k}(t){\phi }_{k}\end{eqnarray} ein Differentialgleichungssystem für die y_k der Form \begin{eqnarray}M{y}_{t}+Ay=b.\end{eqnarray}

Dabei ist A die sogenannte Steifigkeitsmatrix und M = (M_ij) die Massenmatrix, welche definiert ist durch \begin{eqnarray}{M}_{ij}=\displaystyle \int {\phi }_{i}{\phi }_{j}dx.\end{eqnarray}