#!/usr/bin/python

# 10.7

import os, sys

sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]

import scipy as sp

from scipy import signal

from imc_stable1 import imc_stable1

from zpowk import zpowk

from myc2d import myc2d

num = sp.array([1])

den = sp.array([250, 35, 1])

G = signal.lti(num, den)

Ts = 3

B, A, k = myc2d(G, Ts)

alpha = 0.9

k, GiN, GiD = imc_stable1(B, A, k, alpha)

print k, GiN, GiD

zk, dzk = zpowk(k)

Bp, Ap = B, A

t0 = 0

tf = 10

st = 1

Nvar = 0