diff -r 000000000000 -r 0efde00f9229 mpc/gpc_wt.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/mpc/gpc_wt.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,26 @@
+#!/usr/bin/env python
+# 12.4
+import os, sys
+sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
+
+import pylab as pl
+from gpc_N import gpc_N
+
+A, dA = pl.array([1, -0.8]), 1
+B, dB = pl.array([0.4, 0.6]), 1
+k=1
+rho = 0.8
+N1, N2, Nu = 0, 3, 2
+
+K, KH1, KH2, Tc, dTc, Sc, dSc, R1, dR1 = gpc_N(A, dA, B, dB, 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
+