The aim of this paper is to present a numerical approximation for quasilinear parabolic differential functional equations with initial boundary conditions of the Neumann type. The convergence result is proved for a difference scheme with the property that the difference operators approximating mixed derivatives depend on the local properties of the coefficients of the differential equation. A numerical example is given.