This algorithm describes the UVIP3P subroutine subprogram that implements the method of univariate interpolation and smooth curve fitting [1]. The subroutine is written in ANSI Standard FORTRAN [2]. Possible use of a higher degree polynomial is also implemented as the user's option. Unless the option of the use of a higher degree polynomial is exercised by the user, the method has the accuracy of a third-degree (cubic) polynomial even when the input data points are given at unequal intervals.