Develop a Matlab function myCubicFit that calculates the coefficients $a_0$ through $a_3$ of the cubic polynomial $f(x)=$

$a_0+a_{1}x+a_2{x}^2+a_3{x}^3$ that best fits given data points $(x_i,y_i).$ The input arguments to the function must be the two

vectors $x$ and $y$ that contain the values of the data points. The output of the function shall be a four$-$element column

vector a that contains the values of $a_0$ through $a_3,$ and the error of the fit $E$ defined as the sum of the squares of the

residuals. Do NOT use any Matlab build$-$in functions to solve the required system of linear equations. Instead use

myGaussJordan from Module $2.$ Use Matlab's 'operator to calculate the transpose of a matrix and $*$ to calculate the

product of matrices. Inside the function check that the length of the vectors $x$ and $y$ is the same. If not, display an

error message, set the vector a and $E$ to realmax$()$ and exit the function.