Solutions of initial boundary value problems for parabolic functional differential equations are approximated by solutions of implicit difference schemes. The existence and uniqueness of approximate solutions is proved. The proof of the stability is based on a comparison technique with nonlinear estimates of the Perron type for given operators. It is shown that the new methods are considerably better than the explicit difference schemes. Numerical examples are presented.