|
0
|
1 |
#!/usr/bin/python
|
|
|
2 |
# 13.5
|
|
|
3 |
|
|
|
4 |
import scipy as sp
|
|
|
5 |
from lqg import lqg
|
|
|
6 |
from cl import cl
|
|
|
7 |
|
|
|
8 |
A, dA = sp.array([1, -1.4, 0.45]), 2
|
|
|
9 |
C, dC = sp.array([1, -0.5]), 1
|
|
|
10 |
B, dB = 0.5 * sp.array([1, -0.9]), 1
|
|
|
11 |
k = 1
|
|
|
12 |
int1 = 0
|
|
|
13 |
F, dF = 1, 0
|
|
|
14 |
V, dV = 1, 0
|
|
|
15 |
W, dW = 1, 0
|
|
|
16 |
rho = 1
|
|
|
17 |
|
|
|
18 |
R1, dR1, Sc, dSc = lqg (A, dA, B, dB, C, dC, k, rho, V, dV, W, dW, F, dF)
|
|
|
19 |
|
|
|
20 |
Nu, dNu, Du, dDu, Ny, dNy, Dy, dDy, yvar, uvar = \
|
|
|
21 |
cl(A, dA, B, dB, C, dC, k, Sc, dSc, R1, dR1, int1)
|
|
|
22 |
|
|
|
23 |
print R1, dR1, Sc, dSc
|
|
|
24 |
print Nu, dNu, Du, dDu, Ny, dNy, Dy, dDy, yvar, uvar
|