--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/freq/lead_lag.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,23 @@
+#!/usr/bin/python
+# 7.4
+import os, sys
+sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
+
+import pylab as pl
+from scipy import signal
+from bode import bode
+
+w = pl.linspace(0.001,pl.pi,10000)
+a = pl.linspace(0.001,0.999,100)
+
+zero = a
+pole = 0.9*a
+A = pl.vstack((pl.ones_like(zero), -zero))
+A = A.T
+B = pl.vstack((pl.ones_like(pole), -pole))
+B = B.T
+sys = map(signal.lti, A, B)
+
+
+
+#pl.show()