% draw equation of state
close all;

widthpix = 375;
heightpix = 300;

leftmargin = 60;
bottommargin = 55;
topmargin = 10;
rightmargin = 20;

widthtotal = widthpix+leftmargin+rightmargin;
heighttotal = heightpix+topmargin+bottommargin;
left = leftmargin/widthtotal;
bottom = bottommargin/heighttotal;
width = widthpix/widthtotal;
height = heightpix/heighttotal;

psize = 18; % font size
msize = 9; % marker size
fig = figure('rend','painters','pos',[100 100 widthtotal heighttotal]);


datadir = '../data/';
Dstr = '1.5';
Ddble = str2double(Dstr);
dataname = strcat('rhoxiD', Dstr, '_NN.dat');

data = dlmread(strcat(datadir,dataname), '', 1, 0);

%% panel (a)
panelA = axes('Position', [left, bottom, width, height]);

plot(data(:,2), data(:,1), 'linewidth', 1.5);
hold on;
plot([2/sqrt(3), 2/sqrt(3)], [0, 200], '--k');
plot([0.9091,1.0101,1.0526,1.1111,1.14], [7.6938,20.7612,29.8353,61.3047,180.9920], 'b*', 'markersize', msize);
xlabel('$\rho$', 'interpreter', 'latex');
ylabel('$\beta F d$', 'interpreter', 'latex');
xlim([0.6, 1.2]);

%% create table

x0 = 0.63;
y0 = 80;
tbrange = [x0,x0+0.17];
plot(tbrange, [y0+110,y0+110],'-k');
plot(tbrange, [y0+93,y0+93],'-k');
plot(tbrange, [y0,y0],'-k');

text(x0+0.02, y0+102, '$\rho$', 'interpreter', 'latex', 'fontsize', 18);
text(x0+0.1, y0+102, '$\beta F d$', 'interpreter', 'latex', 'fontsize', 18);

text(x0, y0+82, '0.9091', 'interpreter', 'latex', 'fontsize', 18);
text(x0, y0+64, '1.0101', 'interpreter', 'latex', 'fontsize', 18);
text(x0, y0+46, '1.0526', 'interpreter', 'latex', 'fontsize', 18);
text(x0, y0+28, '1.1111', 'interpreter', 'latex', 'fontsize', 18);
text(x0, y0+10, '1.1400', 'interpreter', 'latex', 'fontsize', 18);

text(x0+0.1, y0+82, '7.694', 'interpreter', 'latex', 'fontsize', 18);
text(x0+0.1, y0+64, '20.76', 'interpreter', 'latex', 'fontsize', 18);
text(x0+0.1, y0+46, '29.84', 'interpreter', 'latex', 'fontsize', 18);
text(x0+0.1, y0+28, '61.30', 'interpreter', 'latex', 'fontsize', 18);
text(x0+0.1, y0+10, '181.9', 'interpreter', 'latex', 'fontsize', 18);

set(gca, 'fontname', 'times new roman', 'fontsize', psize);

print('../figures/eos.eps', '-depsc');