#!/usr/bin/python

# 10.6

import os, sys

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

import scipy as sp

from imc_stable1 import imc_stable1

from zpowk import zpowk

B = sp.array([0.51, 1.21])

A = sp.array([1, -0.44])

k = 1

alpha = 0.5

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

print "k =", k

print "GiN =", GiN

print "GiD =", GiD

Bp, Ap = B, A

Ts = 0.1

t0 = 0

tf = 20

Nvar = 0.01