#!/usr/bin/python

import scipy as sp
from lqg import lqg
from cl import cl

A, dA = sp.convolve([1, -0.7], [1, 0.9]), 2
B, dB = sp.convolve([0.079, 0.221], [1, 0.9]), 2
C, dC = sp.array([1, -0.7]), 1
k = 2
int1 = 1
F, dF = sp.array([1, -1]), 1
V, dV = 1, 0
W, dW = 1, 0
rho = 0.1

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

Nu, dNu, Du, dDu, Ny, dNy, Dy, dDy, yvar, uvar = \
     cl(A, dA, B, dB, C, dC, k, Sc, dSc, R1, dR1, int1)
print Nu, dNu, Du, dDu, Ny, dNy, Dy, dDy, yvar, uvar