Neumann Boundary Condition Numerical Method