#!/usr/bin/python
# 13.8
# Solves Example 3.1 of Ahlen and Sternad in Hunt's book

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

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

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

print "R1", R1
print "dR1", dR1
print "Sc", Sc
print "dSc", dSc