USING MATLAB (or you can write just the code)

Write aMatlabfunction [x, t] = fs_approx (a, omega _0, tmin, tmax, n). It is assumed that vector a = [a_- N a_- N + 1 ellipsis a_- 1 a_0 a_1 ellipsis a_N - 1 a_N]. Vector x contains n samples of the Fourier series approximation x_N(t) = sigma _k = - N^N a^k middot e^jk omega _0t between t = tmin(sec) andt=tmax(sec). Note that parameter omega _0 is the angular frequency of the original periodic signal from which Fourier coefficients in a were derived. Output vectort, to be used for plotting, containstime reference for corresponding elements inx. Your program should check whether input vector a is of size2N + 1for some N greaterthanorequalto 0.It should also check whethertmax > fmin. Empty vectors, i.e., x= []andt= []should be returned if input parameters are not valid