In MATLAB

Newton's divided method of interpolation involves finding the coefficients of the formula, b_0 and b_1, given (x_0, y_0), (x_1, y_1), so that they may formulate f_1(x) = b_0 + b_1 (x - x_0), where b_0 = f(x_0) b_1 = f(x_1) - f(x_0)/x_1 - x_0. This is a linear case. For the third order (or the cubic) interpolation, we need to find four coefficients, b_0, b_1, b_2 and b_3. The third order polynomial, given (x_0, y_0), (x_1, y_1), (x_2, y_2), and (x_3, y_3), is f_3(x) = f[x_0] + f[x_1, x_0](x - x_0) + f[x_2, x_1, x_0](x - x_0)(x - x_1) + f[x_3, x_2, x_1, x_0](x - x_0)(x - x_1)(x - x_2) Using the table above, build the MATLAB code (Newton's divided method) for the Newton's divided method of cubic interpolation. Then, find the missing value of velocity at t = 16