--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/lqg/lqg_visc.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,23 @@
+#!/usr/bin/python
+# 13.6
+# Viscosity control problem of MacGregor
+
+import scipy as sp
+from lqg import lqg
+from cl import cl
+
+A, dA = sp.array([1, -0.44]), 1
+B, dB = sp.array([0.51, 1.21]), 1
+F, dF = sp.array([1, -1]), 1
+C, dC = sp.array([1, -0.44]),1
+k = 1
+V, W = 1, 1
+dV, dW = 0, 0
+rho = 1
+int1 = 1
+
+R1,dR1,Sc,dSc = lqg(A,dA,B,dB,C,dC,k,rho,V,dV,W,dW,F,dF)
+
+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
+print "Nu,dNu,Du,dDu,Ny,dNy,Dy,dDy,yvar,uvar"