The publication of this function was inspired by the poor job done by several implementations of the 1d cubic spline publications in Matlab's file exchange repos.

This function was part of one of my homework assignments when I was a master student in NTU back in 2011.

Although it is not a perfect implementation, it is simple to read (the formulation is close the math equations), it can be parallelized (if parfor is available), the boundary conditions can be easily expanded/extended, and it is written in a single script with less than 125 lines of code (including a short intro and example).

Enjoy ! : )