Please provide the Matlab code to solve this question, thank you !!!

Consider the homogeneous equation with non-constant coefficients: y" + (a 2q cos(2t))y = 0, y(0) = 1, y(0) = 0 = This is called Mathieu's equation and it arises in optics. If for example q = 0, this IVP has solution y(t) = cos(Vat). We therefore expect that for q small, there are also periodic solutions. (a) Use ode45 to approximate the solution of Mathieu's equation for q = 0.1 and a = 0.2, 0.5, 0.7,0.8 over the interval 0