1 #!/usr/bin/env python

     2 # 12.7

     3 

     4 import os, sys

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

     6 

     7 import pylab as pl

     8 from gpc_col import gpc_col

     9 

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

    11 B, dB = pl.array([0.51, 1.21]), 1

    12 N = 2

    13 k = 1

    14 C, dC = pl.array([1, -0.44]), 1

    15 rho = 1

    16 

    17 K, KH1, KH2, Tc, dTc, Sc, dSc, R1, dR1 = gpc_col(A, dA, B, dB, C, dC, N, k, rho)

    18 

    19 print "K =", K

    20 print "KH1 =", KH1

    21 print "KH2 =", KH2

    22 print "Tc =", Tc

    23 print "dTc =", dTc

    24 print "Sc =", Sc

    25 print "dSc =", dSc

    26 print "R1 =", R1

    27 print "dR1 =", dR1

    28