--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/lqg/lqg_as1.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,23 @@
+#!/usr/bin/python
+# 13.8
+# Solves Example 3.1 of Ahlen and Sternad in Hunt's book
+
+import scipy as sp
+from lqg import lqg
+from cl import cl
+
+A, dA = sp.array([1, -0.9]), 1
+B, dB = sp.array([0.1, 0.08]), 1
+F, dF = 1, 0
+C, dC = 1, 0
+k = 2
+V, dV = 1, 0
+W, dW = sp.array([1, -1]), 1
+rho = 0.1
+
+R1,dR1,Sc,dSc = lqg_simple(A,dA,B,dB,C,dC,k,rho,V,dV,W,dW,F,dF)
+
+print "R1", R1
+print "dR1", dR1
+print "Sc", Sc
+print "dSc", dSc