Write MATLAB to implement natural cubic spline

Write a MATLAB function to implement the natural cubic spline in Algorithm 3.4. Note that the index should run from 1 instead of 0 that is used in the textbook. This is because MATLAB does not allow index 0. Your function should use the same order of input below. Note that we do not return a,bi,ci,dy as in the algorithm. Instead, we calculate the spline at x. Since the cubic spline is a piecewise polynomial, one needs to locate the subinterval [xi x+1l that contains x before calculating S(x) - StCx). function y = cubic spline (X,Y,x) %X=numbers x 1 x 2 x n % Y-values f (x 1) f (x 2) f (x n) % x-interpolation point x % y- spline value at x, i.e. s(x) % your implementation below end Save your function in cubic_spline.m. Then, compute the natural cubic spline for the data below. -2.4061 | 1.0830 | 0.6440 | 0.4068 | -0.2448 | -0.1158 0 f (x) | -0.3984 | -0.76 1 1 | -0.9688 | -0.9791 | -0.7899 | -0.4397 | 0 (continued) 0.11580.2448 0.4068 0.6440 1.08302.4061 f(x) 0.43y/ e1.7) 0.7899 ga97 0.9791 gaa3313 0.:/611. 0.:Was 0.3984 Natural Cubic Spline ALGORITHM 3.4 To construct the cubic spline interpolant S for the function f, defined at the numbers X0