“…When C 2 splines are considered, the algorithms for choosing free nonlinearity parameters are, in comparison, fairly unsatisfactory; see, e.g., [23] for rational splines and [24] for rational lacunary splines. Recently by Archer and Le Gruyer [1] a direct method for computing C 2 interpolants was given which guarantee convexity preservation on the subinterval [x 1 , x n−1 ] ⊂ I . For reviews on shape preserving interpolation we refer to [9,30].…”