1 #!/usr/bin/python

     2 # 9.4

     3 import scipy as sp

     4 def desired(Ts,rise,epsilon):

     5     Nr = rise/Ts

     6     omega = sp.pi/2/Nr

     7     rho = epsilon**(omega/sp.pi)

     8     phi = sp.array([1 -2*rho*sp.cos(omega), rho**2])

     9     dphi = len(phi)-1

    10     return phi, dphi