1D cubic spline interpolator with irregular monotonic X spacing
No error checking, assumes x>= x at point point llo, updating llo