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