Stiffness matrix ($K_{ij}=a({\varphi }_j,{\varphi }_i)$)

The stiffness matrix represents a variational bilinear form $a$ in a finite element basis. A symmetric coercive bilinear form gives a positive-definite matrix. Local support of the piecewise-linear hat functions produces a sparse matrix, allowing efficient numerical solution.

For $-u^{\prime \prime }+u=f$ with homogeneous Dirichlet boundary conditions, the conforming piecewise-linear hat functions on a uniform grid give a diffusion stiffness matrix $S$ and a mass matrix $M$. Their nonzero entries are $S_{ii}=2/h$, $S_{i,i+1}=-1/h$, $M_{ii}=2h/3$, and $M_{i,i+1}=h/6$, with symmetric counterparts. The sum is a positive-definite matrix because its coefficient quadratic form is $\int ((u_h^{\prime }{)}^2+u_h^2)$. For forcing $f(x)=x$, the load at $x_i$ is $h{x}_i$.