Cubic spline interpolation coefficients
function SplineCoefficients input Real[:] x "x coordinations of interpolating points"; input Real[:] y "y coordinations of interpolating points"; input Real[:] slope "slopes at interpolating points"; output Real[size(x, 1) + 1, 4] a "cubic polynom coefficients of curve segments between interpolating points"; end SplineCoefficients;
author: Ondrej Vacek