In this paper we discuss a nonlocal approximation to the classical heat equation with Neumann boundary conditions. We considerwt∊(x,t)=1∊N+2∫ΩJx-y∊(w∊(y,t)-w∊(x,t))dy+C1∊N∫∂ΩJx-y∊g(y,t)dSy,(x,t)∈Ω‾×(0,T),w(x,0)=u0(x),x∈Ω‾,and we show that the corresponding solutions, w∊, converge to the classical solution of the local heat equation vt=Δv with Neumann boundary conditions, ∂v∂n(x,t)=g(x,t), and initial condition v(0)=u0, as the parameter ∊ goes to zero. The obtained convergence is in the weak star on L∞ topology.