This paper deals with the long-time behaviour of numerical
solutions of neutral delay differential equations that have stable
hyperbolic periodic orbits. It is shown that Runge--Kutta
discretizations of such equations have attractive invariant closed
curves which approximate the periodic orbit with the full order of
the method, in spite of the lack of a finite-time smoothing
property of the flow.