1 #!/usr/bin/python

     2 # 13.6

     3 # Viscosity control problem of MacGregor

     4 

     5 import scipy as sp

     6 from lqg import lqg

     7 from cl import cl

     8 

     9 A, dA = sp.array([1, -0.44]), 1

    10 B, dB = sp.array([0.51, 1.21]), 1

    11 F, dF = sp.array([1, -1]), 1

    12 C, dC = sp.array([1, -0.44]),1

    13 k = 1

    14 V, W = 1, 1

    15 dV, dW = 0, 0

    16 rho = 1

    17 int1 = 1

    18 

    19 R1,dR1,Sc,dSc = lqg(A,dA,B,dB,C,dC,k,rho,V,dV,W,dW,F,dF)

    20 

    21 Nu,dNu,Du,dDu,Ny,dNy,Dy,dDy,yvar,uvar = cl(A,dA,B,dB,C,dC,k,Sc,dSc,R1,dR1,int1)

    22 print Nu,dNu,Du,dDu,Ny,dNy,Dy,dDy,yvar,uvar

    23 print "Nu,dNu,Du,dDu,Ny,dNy,Dy,dDy,yvar,uvar"