Investigating Neumann BC using ghost points

[shivin9]

Ghost point approach gives lesser accuracy on the boundary but a more accurate solution elsewhere as compared to the one-sided Fornberg coefficient approach.
So we can think of implementing that later.