Initial commit.
#!/usr/bin/python
# 13.2
import scipy as sp
from spec import spec
A = sp.convolve([-0.5, 1], [-0.9, 1])
dA = 2
B = 0.5*sp.array([-0.9, 1])
dB = 1
rho = 1
r, beta1, sigma = spec(A, dA, B, dB, rho)
print "sigma", sigma
print "beta1", beta1
print "r", r