Computation of Neumann Localised Boundary Domain Integral Equations