In this work we prove the existence and uniqueness of solutions as well as the validation of a comparison principle for a discrete model associated to a nonlocal diffusion problem with Neumann conditions. We show that the solutions for the discrete model converge to the solutions of the continuous model when the mesh parameter goes to zero. Finally, we perform some numerical experiments.