#!/usr/bin/env python

# 12.4

import os, sys

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

import pylab as pl

from gpc_N import gpc_N

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

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

k=1

rho = 0.8

N1, N2, Nu = 0, 3, 2

K, KH1, KH2, Tc, dTc, Sc, dSc, R1, dR1 = gpc_N(A, dA, B, dB, 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