Write a program in MATLAB that will do the following:

a. Make a Nichols plot of an open-loop transfer function

b. Allow the user to read the Nichols plot display and enter the value of M_{p}

c. Make closed-loop magnitude and phase plots

d. Display the expected values of percent overshoot, settling time, and peak time

e. Plot the closed-loop step response Test your program on a unity feedback system with the forward-path transfer function

and explain any discrepancies.