Do the question in python 3.

Starting point of the code: def polyfit_4(x,f): #Add code here

Consider the four polynomials Po (x) 1, P (x)-1+ P2(z) = 1 + z +z?" and hence define the polynomial p(x) aoPo (x)+ap(x)+a2pe (x) +asps (x) Given four coordinates (ro fo). (i,fi). (x2,fa), and (x3, fs), with none of the four r, equal, let the coefficients ao, a1,a2,a be defined such that p(x) is the interpolating polynomial for which p(zi) = fi foriE {0,1,2,3} Write a function named polyfit_4' which accepts two arguments 1. a length 4 NumPy vector 'x' of floating point values, with 'xU containing the value of 2. a length 4 NumPy vector 'P of floating point values, with 1l containing the value of f and returns . a length 4 NumPy vector 'a' of floating point values, with all' containing the value of ai