Write a Matlab function called compute Bode () that does the following: a. Accepts as its input argument a transfer function G(jw). b. Returns two vectors, each as a function of w (in rad/s): (1) |G(jw)| (in dB) and (2) phase angle, (in deg), where phase angle is negative if the output lags behind the input. c. Also returns 4 values: gain crossover frequency (in rad/s), phase crossover frequency (in rad/s), gain margin (in dB), and phase margin (in deg).