write a general function, in MATLAB, to fit a sinusoidal curve of the type y = a0 + a1 sin(t) + a2 cos(t) to a data set (t,y) with the following format: function [A,r] = fourier(t,y,) where A is a vector that contains the coefficients a0, a1, and a2, and r is the correlation coefficient, and is the angular frequency of the oscillations (if you are not familiar with this parameter, please look it up, and, in particular, see how it relates to the oscillations period). This function should implement the general linear least-squares theory studied in class and therefore use matrix operations (only) to compute the coefficients; a loop should NOT be used. The function should plot the data as well as the sinusoidal fit on the same plot.