1 #!/usr/bin/env python

     2 # 11.11

     3 # Viscosity control problem of MacGregor

     4 import os, sys

     5 sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]

     6 

     7 import pylab as pl

     8 import mv_nm as mv_nm

     9 from cl import cl

    10 

    11 A, dA = pl.array([1, -0.44]), 1 

    12 B, dB = pl.array([0.51, 1.21]), 1

    13 C, dC = pl.array([1, -0.44]), 1

    14 k, int1 = 1, 1

    15 

    16 Sc, dSc, Rc, dRc = mv_nm(A, dA, B, dB, C, dC, k, int1)

    17 

    18 [Nu, dNu, Du, dDu, Ny, dNy, Dy, dDy, yvar, uvar] = \

    19 cl(A, dA, B, dB, C, dC, k, Sc, dSc, Rc, dRc, int1)