The Lagrange polynomials form the normal basis of most Finite Elements Method (FEM) software. There are other polynomials which are used as well, but these are the workhorse of most solvers.
(proof-reading through HN is a mildly embarrassing process, sorry about that! I do go over these posts and proof-read them several times myself before publishing)
the last matrix before the appendix is not the identity matrix, right now the matrix is: \begin{bmatrix}
1 & 0 & 0 & \dots & 0\\
1 & 0 & 0 & \dots & 0\\
1 & 0 & 0 & \dots & 0\\
\vdots & \vdots & \vdots & \ddots &\vdots \\
1 & 0 & 0 & \dots & 1
\end{bmatrix}
(proof-reading through HN is a mildly embarrassing process, sorry about that! I do go over these posts and proof-read them several times myself before publishing)