function [I,J,n] = gibbstablePT
% -----------------------------------------------------------------
% Author: Philippe Daigle
% Date: 10 mai 2002
% Numerical values of the coefficients 
% and exponents of the dimensionless 
% Gibbs fre energy for region 1
% -----------------------------------------------------------------
% i = 1
I(1) = 0;
J(1) = -2;
n(1) = 0.14632971213167;
% --------------------------------
% i = 2
I(2) = 0;
J(2) = -1;
n(2) = -0.84548187169114;
% --------------------------------
% i = 3
I(3) = 0;
J(3) = 0;
n(3) = -0.37563603672040*10^1;
% --------------------------------
% i = 4
I(4) = 0;
J(4) = 1;
n(4) = 0.33855169168385*10^1;
% --------------------------------
% i = 5
I(5) = 0;
J(5) = 2;
n(5) = -0.95791963387872;
% --------------------------------
% i = 6
I(6) = 0;
J(6) = 3;
n(6) = 0.15772038513228;
% --------------------------------
% i = 7
I(7) = 0;
J(7) = 4;
n(7) = -0.16616417199501*10^-1;
% --------------------------------
% i = 8
I(8) = 0;
J(8) = 5;
n(8) = 0.81214629983568*10^-3;
% --------------------------------
% i = 9
I(9) = 1;
J(9) = -9;
n(9) = 0.28319080123801*10^-3;
% --------------------------------
% i = 10
I(10) = 1;
J(10) = -7;
n(10) = -0.60706301565874*10^-3;
% --------------------------------
% i = 11
I(11) = 1;
J(11) = -1;
n(11) = -0.18990068218419*10^-1;
% --------------------------------
% i = 12
I(12) = 1;
J(12) = 0;
n(12) = -0.32529748770505*10^-1;
% --------------------------------
% i = 13
I(13) = 1;
J(13) = 1;
n(13) = -0.21841717175414*10^-1;
% --------------------------------
% i = 14
I(14) = 1;
J(14) = 3;
n(14) = -0.52838357969930*10^-4;
% --------------------------------
% i = 15
I(15) = 2;
J(15) = -3;
n(15) = -0.47184321073267*10^-3;
% --------------------------------
% i = 16
I(16) = 2;
J(16) = 0;
n(16) = -0.30001780793026*10^-3;
% --------------------------------
% i = 17
I(17) = 2;
J(17) = 1;
n(17) = 0.47661393906987*10^-4;
% --------------------------------
% i = 18
I(18) = 2;
J(18) = 3;
n(18) = -0.44141845330846*10^-5;
% --------------------------------
% i = 19
I(19) = 2;
J(19) = 17;
n(19) = -0.72694996297594*10^-15;
% --------------------------------
% i = 20
I(20) = 3;
J(20) = -4;
n(20) = -0.31679644845054*10^-4;
% --------------------------------
% i = 21
I(21) = 3;
J(21) = 0;
n(21) = -0.28270797985312*10^-5;
% --------------------------------
% i = 22
I(22) = 3;
J(22) = 6;
n(22) = -0.85205128120103*10^-9;
% --------------------------------
% i = 23
I(23) = 4;
J(23) = -5;
n(23) = -0.22425281908000*10^-5;
% --------------------------------
% i = 24
I(24) = 4;
J(24) = -2;
n(24) = -0.65171222895601*10^-6;
% --------------------------------
% i = 25
I(25) = 4;
J(25) = 10;
n(25) = -0.14341729937924*10^-12;
% --------------------------------
% i = 26
I(26) = 5;
J(26) = -8;
n(26) = -0.40516996860117*10^-6;
% --------------------------------
% i = 27
I(27) = 8;
J(27) = -11;
n(27) = -0.12734301741641*10^-8;
% --------------------------------
% i = 28
I(28) = 8;
J(28) = -6;
n(28) = -0.17424871230634*10^-9;
% --------------------------------
% i = 29
I(29) = 21;
J(29) = -29;
n(29) = -0.68762131295531*10^-18;
% --------------------------------
% i = 30
I(30) = 23;
J(30) = -31;
n(30) = 0.14478307828521*10^-19;
% --------------------------------
% i = 31
I(31) = 29;
J(31) = -38;
n(31) = 0.26335781662795*10^-22;
% --------------------------------
% i = 32
I(32) = 30;
J(32) = -39;
n(32) = -0.11947622640071*10^-22;
% --------------------------------
% i = 33
I(33) = 31;
J(33) = -40;
n(33) = 0.18228094581404*10^-23;
% --------------------------------
% i = 34
I(34) = 32;
J(34) = -41;
n(34) = -0.93537087292458*10^-25;
% --------------------------------