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