#!/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