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