What is the MATLAB code for this?

The world population in selected years is given in the table below. Find a polynomial of the form P = a_n t^n + a_n-1 t^n-1 + .. + a_0 that best fit the given data (t is the number of years after 1850 and P is the population in millions). Make a plot of the world population versus the year that shows the given data points and the polynomial model. Using "spline" method and the polynomial model, interpolate the population in 2017. If the actual world population in 2017 is 7515millions, what is the percentage error between the actual value and the interpolated value from part c