In [5]:
import numpy as np
from matplotlib import pyplot as pp

w = np.arange(0, 3, .001)
A = 1/(1 + .1j*w - w*w)

amplitude = np.absolute(A)
phase = np.angle(A)
In [3]:
pp.plot(w, amplitude)
pp.xlabel('Omega')
pp.ylabel('Amplitude')
pp.axis([0,3,-2,12])
pp.show()
In [4]:
pp.plot(w, phase)
pp.xlabel('Omega')
pp.ylabel('Phase')
pp.axis([0,3,-5,2])
pp.show()
In [ ]: