Cubic spline interpolation function
function Spline input Real[:] x "x coordinations of interpolating points"; input Real[:, 4] a "cubic polynom coefficients of curve segments between interpolating points"; input Real xVal "input value of x to calculate y value"; output Real yVal "y value at xVal"; end Spline;
author: Ondrej Vacek