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