1 #!/usr/bin/env python

     2 # 12.4

     3 import os, sys

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

     5 

     6 import pylab as pl

     7 from gpc_N import gpc_N

     8 

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

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

    11 k=1

    12 rho = 0.8

    13 N1, N2, Nu = 0, 3, 2

    14 

    15 K, KH1, KH2, Tc, dTc, Sc, dSc, R1, dR1 = gpc_N(A, dA, B, dB, k, N1, N2, Nu, rho)

    16 

    17 print "K =", K

    18 print "KH1 =", KH1

    19 print "KH2 =", KH2

    20 print "Tc =", Tc

    21 print "dTc =", dTc

    22 print "Sc =", Sc

    23 print "dSc =", dSc

    24 print "R1 =", R1

    25 print "dR1 =", dR1

    26