Fritch-Carlson interpolation

Syntax: vout = INTERPOLATE(x,y,xout,'FC')

If the FC keyword is entered, the interpolant is calculated using the Fritsch and Carlson method of monotone piecewise cubic interpolation. This algorithm produces a visually pleasing interpolant, that is, the interpolating curve has no extraneous "bumps" or "wiggles".

For an explanation of this method, see:

SIAM Journal of Numerical Analysis, volume 17, number 2, April 1980.

  Lagrange interpolation
  SPLINTERP function