Write down a code to calculate a natural spline S(x) which interpolates a collection of given points (x0, y0),(x1, y1), , , ,(xn, yn) where x0 < x1 < x2 < . . . < xn. Your code should have a triadiagonal solver for the resulting linear system of equations (you're not allowed to use Matlab's n operator to solve the linear system).