diff -r 000000000000 -r 0efde00f9229 imc/visc_imc1.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/imc/visc_imc1.py Fri May 27 14:24:59 2011 +0530 @@ -0,0 +1,23 @@ +#!/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