--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/minv/mv_visc.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,19 @@
+#!/usr/bin/env python
+# 11.11
+# Viscosity control problem of MacGregor
+import os, sys
+sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
+
+import pylab as pl
+import mv_nm as mv_nm
+from cl import cl
+
+A, dA = pl.array([1, -0.44]), 1 
+B, dB = pl.array([0.51, 1.21]), 1
+C, dC = pl.array([1, -0.44]), 1
+k, int1 = 1, 1
+
+Sc, dSc, Rc, dRc = mv_nm(A, dA, B, dB, C, dC, k, int1)
+
+[Nu, dNu, Du, dDu, Ny, dNy, Dy, dDy, yvar, uvar] = \
+cl(A, dA, B, dB, C, dC, k, Sc, dSc, Rc, dRc, int1)