Consider the following method for solving the ODE y = f(y,t) y=yn + hf (yn, ta) Y+1=Yn+haf (y..t+1)+(1-a) f(un.tu)) where 0 a 1 (a) Apply this method to y'=iwy, where is a real number, and find the optimal value of a for stability. What is the largest time step you can take with this optimal value of a? (b) With the value of a obtained in part (a), we solve the system y'= iwy with y(0) = 1 and step size h = 1/w. What are the amplitude and phase error after 100 steps? (e) Now find the value of a that gives you maximum possible accuracy. (d) For the value of a obtained in part (e), what are the stability characteristics of the method when applied to the ODE (w real).