--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/freq/lead_exp.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,21 @@
+#!/usr/bin/python
+# 7.3
+import os, sys
+sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
+
+import pylab as pl
+from scipy import signal
+from bode import bode
+
+pol1 = [1, -0.9]
+pol2 = [1, -0.8]
+G1 = signal.lti(pol1,[1, 0])
+G2 = signal.lti(pol2,[1, 0])
+w = pl.linspace(0.001,0.5,1000)
+bode(G1, w)
+bode(G2, w)
+pl.figure()
+G = signal.lti(pol1,pol2)
+bode(G,w)
+
+pl.show()