Model order reduction for delay systems by iterative interpolation