function [I,J,n] = region3_hfe
% Description:      Loads the numerical values of the coefficients and
%                   exponents of the dimensionless Helmholtz free energy
%                   for Region 3 (Table 30, page 30)                 
% -------------------------------------------------------------------------
% Syntax:           [I,J,n] = region3_hfe
% -------------------------------------------------------------------------
% 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:12)  = 1;
I(13:18) = 2;
I(19:23) = 3;
I(24:27) = 4;
I(28:30) = 5;
I(31:33) = 6;
I(34)    = 7;
I(35)    = 8;
I(36:37) = 9;
I(38:39) = 10;
I(40)    = 11;
%----------------------------------------------------------
J(1)  = 0;
J(2)  = 0;
J(3)  = 1;
J(4)  = 2;
J(5)  = 7;
J(6)  = 10;
J(7)  = 12;
J(8)  = 23;
J(9)  = 2;
J(10) = 6;
J(11) = 15;
J(12) = 17;
J(13) = 0;
J(14) = 2;
J(15) = 6;
J(16) = 7;
J(17) = 22;
J(18) = 26;
J(19) = 0;
J(20) = 2;
J(21) = 4;
J(22) = 16;
J(23) = 26;
J(24) = 0;
J(25) = 2;
J(26) = 4;
J(27) = 26;
J(28) = 1;
J(29) = 3;
J(30) = 26;
J(31) = 0;
J(32) = 2;
J(33) = 26;
J(34) = 2;
J(35) = 26;
J(36) = 2;
J(37) = 26;
J(38) = 0;
J(39) = 1;
J(40) = 26;
%----------------------------------------------------------
n(1)  =  0.10658070028513*10^1;
n(2)  = -0.15732845290239*10^2;
n(3)  =  0.20944396974307*10^2;
n(4)  = -0.76867707878716*10^1;
n(5)  =  0.26185947787954*10^1;
n(6)  = -0.28080781148620*10^1;
n(7)  =  0.12053369696517*10^1;
n(8)  = -0.84566812812502*10^-2;
n(9)  = -0.12654315477714*10^1;
n(10) = -0.11524407806681*10^1;
n(11) =  0.88521043984318;
n(12) = -0.64207765181607;
n(13) =  0.38493460186671;
n(14) = -0.85214708824206;
n(15) =  0.48972281541877*10^1;
n(16) = -0.30502617256965*10^1;
n(17) =  0.39420536879154*10^-1;
n(18) =  0.12558408424308;
n(19) = -0.27999329698710;
n(20) =  0.13899799569460*10^1;
n(21) = -0.20189915023570*10^1;
n(22) = -0.82147637173963*10^-2;
n(23) = -0.47596035734923;
n(24) =  0.43984074473500*10^-1;
n(25) = -0.44476435428739;
n(26) =  0.90572070719733;
n(27) =  0.70522450087967;
n(28) =  0.10770512626332;
n(29) = -0.32913623258954;
n(30) = -0.50871062041158;
n(31) = -0.22175400873096*10^-1;
n(32) =  0.94260751665092*10^-1;
n(33) =  0.16436278447961;
n(34) = -0.13503372241348*10^-1;
n(35) = -0.14834345352472*10^-1;
n(36) =  0.57922953628084*10^-3;
n(37) =  0.32308904703711*10^-2;
n(38) =  0.80964802996215*10^-4;
n(39) = -0.16557679795037*10^-3;
n(40) = -0.44923899061815*10^-4;