--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/lqg/lqg_julien1.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,22 @@
+#!/usr/bin/python
+
+import scipy as sp
+from lqg import lqg
+from cl import cl
+
+A, dA = sp.convolve([1, -0.7], [1, 0.9]), 2
+B, dB = sp.convolve([0.079, 0.221], [1, 0.9]), 2
+C, dC = sp.array([1, -0.7]), 1
+k = 2
+int1 = 1
+F, dF = sp.array([1, -1]), 1
+V, dV = 1, 0
+W, dW = 1, 0
+rho = 0.1
+
+R1, dR1, Sc, dSc = lqg (A, dA, B, dB, C, dC, k, rho, V, dV, W, dW, F, dF)
+#print R1, dR1, Sc, dSc
+
+Nu, dNu, Du, dDu, Ny, dNy, Dy, dDy, yvar, uvar = \
+     cl(A, dA, B, dB, C, dC, k, Sc, dSc, R1, dR1, int1)
+print Nu, dNu, Du, dDu, Ny, dNy, Dy, dDy, yvar, uvar