1 #!/usr/bin/python

     2 

     3 import scipy as sp

     4 from lqg import lqg

     5 from cl import cl

     6 

     7 A, dA = sp.convolve([1, -0.7], [1, 0.9]), 2

     8 B, dB = sp.convolve([0.079, 0.221], [1, 0.9]), 2

     9 C, dC = sp.array([1, -0.7]), 1

    10 k = 2

    11 int1 = 1

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

    13 V, dV = 1, 0

    14 W, dW = 1, 0

    15 rho = 0.1

    16 

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

    18 #print R1, dR1, Sc, dSc

    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 print Nu, dNu, Du, dDu, Ny, dNy, Dy, dDy, yvar, uvar