How to model transfer functions using Python (without using higher level packages) and plot the magnitude and phase response?

Question

I\'m trying to model a transfer function in Python and thought I could do it by simply plotting the transfer function at many frequencies. This seemed to work for a 2nd order LP