--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/mpc/gpc_ex12.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,27 @@
+#!/usr/bin/env python
+# 12.3
+# Camacho and Bordon's GPC example; model formation
+import os, sys
+sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
+
+import pylab as pl
+from gpc_bas import gpc_bas
+
+A, dA = pl.array([1, -0.8]), 1
+B, dB = pl.array([0.4, 0.6]), 1
+N=3
+k=1
+rho = 0.8
+
+K, KH1, KH2, Tc, dTc, Sc, dSc, R1, dR1 = gpc_bas(A, dA, B, dB, N, k, 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
+