FORTRAN subroutine for plotting a cubic spline function