Now we must add up all the contributions from each element into a global matrix and global vector.

where Nel is equal to the total number of elements in the discretized solution domain and i represents the node numbers (vertices).