Modified matrix and loading vector for Dirichlet boundary.
Added contribution from Neumann condition to loading vector.