1 #!/usr/bin/python

     2 # 10.6

     3 import os, sys

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

     5 

     6 import scipy as sp

     7 from imc_stable1 import imc_stable1

     8 from zpowk import zpowk

     9 

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

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

    12 k = 1

    13 alpha = 0.5

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

    15 print "k =", k

    16 print "GiN =", GiN

    17 print "GiD =", GiD

    18 

    19 Bp, Ap = B, A

    20 Ts = 0.1

    21 t0 = 0

    22 tf = 20

    23 Nvar = 0.01