This PR adds some type conversations for algebraic operations involving I/O systems, motivated by an example I was putting together where I wanted to generate a Nyquist plot for a linearized nonlinear I/O system (typeLinearIOSystem) and a simple lead compensator (typeTransferFunction). UsingP * C generated a type conversation error butC * P was OK as wasct.tf(P) * C. It seemed annoying to have to do the explicit typecasting, so I decided to implement some of the operations that were needed to makeP * C work (generates aLinearICSystem, which is a subclass ofStateSpace).

The following functionality has been added:

• Added__sub__, so that you can write G1 - G2 for I/O systems
• Added conversions for scalars, arrays, state space, and transfer function systems so that they can be combined with I/O systems on either the left or right.
• Set the Numpy__array_priority__ value so that operations involving ndarrays are handled as I/O system operations and I/O system operations take precedence over LTI operations involving ndarrays.

I also did some PEP8 cleanup oniosys.py while I was at it, and put in unit tests for the code.