Code must be in MATLAB

3. The Duffing Equation is a nonlinear, second order differential equation of the form The equation can be rewritten as the following system of first order differential equations y' =-5y-az-823 + cos(wt) (a) Let --, -1, -0.3, -0.37, and w-1.2. Write an m-file function for the system of differential equations. (b) Solve the system of equations over the interval t E 10, 38) with a step size of 0.01 for the initial conditions (0)1, y(0)-0 (c) Plot the phase-plane solution found in (b). Label the axes of the plot. (d) Plot the component-curve solutions found in (b). The z(t) and the y(t) curves should be plotted together on one graph. Use a legend to distinguish the curves and be sure to label the axes of the plot