1 #!/usr/bin/env python

     2 # 12.3

     3 # Camacho and Bordon's GPC example; model formation

     4 import os, sys

     5 sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]

     6 

     7 import pylab as pl

     8 from gpc_bas import gpc_bas

     9 

    10 A, dA = pl.array([1, -0.8]), 1

    11 B, dB = pl.array([0.4, 0.6]), 1

    12 N=3

    13 k=1

    14 rho = 0.8

    15 

    16 K, KH1, KH2, Tc, dTc, Sc, dSc, R1, dR1 = gpc_bas(A, dA, B, dB, N, k, rho)

    17 

    18 print "K =", K

    19 print "KH1 =", KH1

    20 print "KH2 =", KH2

    21 print "Tc =", Tc

    22 print "dTc =", dTc

    23 print "Sc =", Sc

    24 print "dSc =", dSc

    25 print "R1 =", R1

    26 print "dR1 =", dR1

    27