--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/mpc/gpc_wtc.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,28 @@
+#!/usr/bin/env python
+# 12.7
+
+import os, sys
+sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
+
+import pylab as pl
+from gpc_Nc import gpc_Nc
+
+A, dA = pl.array([1, -0.44]), 1
+B, dB = pl.array([0.51, 1.21]), 1
+C, dC = pl.array([1, -0.44]), 1
+k = 1
+rho = 1
+N1, N2, Nu = 0, 2, 0
+K, KH1, KH2, Tc, dTc, Sc, dSc, R1, dR1 = gpc_Nc(A, dA, B, dB, C, dC, k, N1, N2, Nu, rho)
+
+print "K =", K
+print "KH1 =", KH1
+print "KH2 =", KH2
+print "Tc =", Tc
+print "dTc =", dTc
+print "Sc =", Sc
+print "dSc =", dSc
+print "R1 =", R1
+print "dR1 =", dR1
+
+