diff -r 000000000000 -r 0efde00f9229 minv/mv_visc.py --- /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)