Numerical periodic orbits of neutral delay differential equations