#!/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)