Problem 3 needs to be done before completing 4.

Write a MATLAB function, called Newton_poly that inputs a set of data points (x; y) =(datx, daty), a set x of numbers at which to interpolate, and outputs the polynomial interpolant, y, evaluated at x using Newton polynomial interpolation. Your function header should look something like: function y = Newton_poly(x, datx, daty) Use the code you developed in Problem 3 to interpolate the function f(x) = 1/a + x^2 using the data points datx=-5:1:5. Interpolate at the points x=-5:0.001:5. Plot the results and compare them to what you got in problem 2(b). Explain why you get what you get