Dirichlet boundary condition fixed (row modification v.0)
Added Dirichlet second order
Dirichlet condition via penalization method