author | Puneeth Chaganti <punchagan@fossee.in> |
Fri, 27 May 2011 14:24:59 +0530 | |
changeset 0 | 0efde00f9229 |
permissions | -rw-r--r-- |
#!/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