A multi-segment subdivision scheme for an arbitrary univariate real polynomial on a finite interval is established. The method uses Newton’s interpolation in combination with the AIFS (Affine invariant Iterated Function Systems) to construct a fractal-type algorithm that products polynomial geometry.