--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/python/desired.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,10 @@
+#!/usr/bin/python
+# 9.4
+import scipy as sp
+def desired(Ts,rise,epsilon):
+    Nr = rise/Ts
+    omega = sp.pi/2/Nr
+    rho = epsilon**(omega/sp.pi)
+    phi = sp.array([1 -2*rho*sp.cos(omega), rho**2])
+    dphi = len(phi)-1
+    return phi, dphi