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