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.