function [Ir,Jr,nr] = region5_gfe_r
% Description:      Loads the numerical values of the coefficients and
%                   exponents of the residual part of the dimensionless 
%                   Gibbs free energy for Region 5 (Table 38, page 37)                 
% -------------------------------------------------------------------------
% Syntax:           [Ir,Jr,nr] = region5_gfe_r
% -------------------------------------------------------------------------
% Output:           Ir          Array for the exponent parameter Ir 
%                   Jr          Array for the exponent parameter Jr  
%                   nr          Array for the coefficient parameter nr

Ir(1:3) = 1;
Ir(4)   = 2;
Ir(5)   = 3;
%----------------------------------------------------------
Jr(1)  = 0;
Jr(2)  = 1;
Jr(3)  = 3;
Jr(4)  = 9;
Jr(5)  = 3;
%----------------------------------------------------------
nr(1) = -0.12563183589592*10^-3;
nr(2) =  0.21774678714571*10^-2;
nr(3) = -0.45942820899910*10^-2;
nr(4) = -0.39724828359569*10^-5;
nr(5) =  0.12919228289784*10^-6;