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_Nc import gpc_Nc

     9 

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

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

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

    13 k = 1

    14 rho = 1

    15 N1, N2, Nu = 0, 2, 0

    16 K, KH1, KH2, Tc, dTc, Sc, dSc, R1, dR1 = gpc_Nc(A, dA, B, dB, C, dC, k, N1, N2, Nu, 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 

    28