function [I,J,n] = region1_gfe
% Description:      Loads the numerical values of the coefficients and
%                   exponents of the dimensionless Gibbs free energy for
%                   Region 1 (Table 2, page 7)
% -------------------------------------------------------------------------
% Syntax:           [I,J,n] = region1_gfe
% -------------------------------------------------------------------------
% Output:           I           Array for the exponent parameter I             
%                   J           Array for the exponent parameter J
%                   n           Array for the coefficient parameter n

I(1:8)   = 0;
I(9:14)  = 1;
I(15:19) = 2;
I(20:22) = 3;
I(23:25) = 4;
I(26)    = 5;
I(27:28) = 8;
I(29)    = 21;
I(30)    = 23;
I(31)    = 29;
I(32)    = 30;
I(33)    = 31;
I(34)    = 32;
%----------------------------------------------------------
J(1)  = -2;
J(2)  = -1;
J(3)  =  0;
J(4)  =  1;
J(5)  =  2;
J(6)  =  3;
J(7)  =  4;
J(8)  =  5;
J(9)  = -9;
J(10) = -7;
J(11) = -1;
J(12) =  0;
J(13) =  1;
J(14) =  3;
J(15) = -3;
J(16) =  0;
J(17) =  1;
J(18) =  3;
J(19) =  17;
J(20) = -4;
J(21) =  0;
J(22) =  6;
J(23) = -5;
J(24) = -2;
J(25) =  10;
J(26) = -8;
J(27) = -11;
J(28) = -6;
J(29) = -29;
J(30) = -31;
J(31) = -38;
J(32) = -39;
J(33) = -40;
J(34) = -41;
%----------------------------------------------------------
n(1)  =  0.14632971213167;
n(2)  = -0.84548187169114;
n(3)  = -0.37563603672040*10^1;
n(4)  =  0.33855169168385*10^1;
n(5)  = -0.95791963387872;
n(6)  =  0.15772038513228;
n(7)  = -0.16616417199501*10^-1;
n(8)  =  0.81214629983568*10^-3;
n(9)  =  0.28319080123801*10^-3;
n(10) = -0.60706301565874*10^-3;
n(11) = -0.18990068218419*10^-1;
n(12) = -0.32529748770505*10^-1;
n(13) = -0.21841717175414*10^-1;
n(14) = -0.52838357969930*10^-4;
n(15) = -0.47184321073267*10^-3;
n(16) = -0.30001780793026*10^-3;
n(17) =  0.47661393906987*10^-4;
n(18) = -0.44141845330846*10^-5;
n(19) = -0.72694996297594*10^-15;
n(20) = -0.31679644845054*10^-4;
n(21) = -0.28270797985312*10^-5;
n(22) = -0.85205128120103*10^-9;
n(23) = -0.22425281908000*10^-5;
n(24) = -0.65171222895601*10^-6;
n(25) = -0.14341729937924*10^-12;
n(26) = -0.40516996860117*10^-6;
n(27) = -0.12734301741641*10^-8;
n(28) = -0.17424871230634*10^-9;
n(29) = -0.68762131295531*10^-18;
n(30) =  0.14478307828521*10^-19;
n(31) =  0.26335781662795*10^-22;
n(32) = -0.11947622640071*10^-22;
n(33) =  0.18228094581404*10^-23;
n(34) = -0.93537087292458*10^-25;