1 #!/usr/bin/python

     2 # 13.2

     3 

     4 import scipy as sp

     5 from spec import spec

     6 

     7 A = sp.convolve([-0.5, 1], [-0.9, 1])

     8 dA = 2

     9 B = 0.5*sp.array([-0.9, 1])

    10 dB = 1

    11 rho = 1

    12 

    13 r, beta1, sigma = spec(A, dA, B, dB, rho)

    14 print "sigma", sigma

    15 print "beta1", beta1

    16 print "r", r

    17