Write a code to compute a natural spline s(x) which interpolates a collection of given points (x0; y0); (x1; y1); : : : (xn; yn) where x0 < x1 < x2 < : : : < xn (do not assume they are equidistributed). 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).