#!/usr/bin/env python
# 12.7

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

import pylab as pl
from gpc_Nc import gpc_Nc

A, dA = pl.array([1, -0.44]), 1
B, dB = pl.array([0.51, 1.21]), 1
C, dC = pl.array([1, -0.44]), 1
k = 1
rho = 1
N1, N2, Nu = 0, 2, 0
K, KH1, KH2, Tc, dTc, Sc, dSc, R1, dR1 = gpc_Nc(A, dA, B, dB, C, dC, k, N1, N2, Nu, 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