Added function to determine Nyquist plot

Code based on the lti.bode() function.
n=10000 decided through trial and error. Might be somewhat high for
simple systems.