Two algorithms for establishing PACSIC (parametric cubic spline
interpolation curve) for any given dataset as an open or closed PLF
(piecewise linear function) were established. Both algorithms do not
require any prescribed tangent directions or any boundary conditions,
though these directions can also be pre-defined if necessary. The
direction at each vertex was simply determined by the two neighboring
vertices. The LS (least square) G1 PACSIC, in particular, can be applied
to approximate any 2d or 3d dataset. Simple examples were used to
illustrate the flexibility of the algorithms, and various illustrative
examples were also shown to demonstrate the efficiency of both
algorithms. For completeness, the G2 PACSIC was briefly investigated,
where the G2 conditions were given in terms of elastic parameters.
Certainly, similar algorithms for generic PACSIC and LS PACSIC can also
be developed if the tangent directions at vertices are evaluated by
using three or more neighboring vertices; or even by raising the degree
from cubic to quartic.