Introduction - If you have any usage issues, please Google them yourself
The basic idea is to take B-spline function and basis functions in [a, b] in the evenly divided intervals, xi = a+ ih, h = (ba)/n, because the cubic spline function is n+3 dimensional space , so the expansion points to the x-1, xn+1, then any cubic spline function can be expressed as a linear combination: