Added NewtonCotes formula and neumann (not working)
Added neumann boundary condition on the whole boundary
Dirichlet boundary testing day 1
Merge branch 'sparse', remote-tracking branch 'origin/sparse' into sparse