Add automatic computation of initial phase #1093

JuliaControl added a nice function to start the phase on a Bode plot based on the number of integrators at the origin. It might be nice to add that to python-control.

There's a pretty good description in the code comments:

https://github.com/JuliaControl/ControlSystems.jl/pull/950/files