A method is described for constructing an interpolant to a set of arbi
trary data points (x(i), y(i)), i = 1, 2, ..., n. The constructed inte
rpolant is a piecewise parametric cubic polynomial and satisfies C1 co
ntinuity, and it reproduces all parametric polynomials of degree two o
r less exactly. The experiments to compare the new method with Bessel
method and spline method are also shown.