|
0
|
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"
|