--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/lqg/spec_ex.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,17 @@
+#!/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
+