function plotEOS()

close all
load finalRho_1.5.dat
load finalRho_2.0.dat
load finalRho_2.2.dat
load finalRho_2.4.dat
load eos1.5.dat
load eos2.2.dat
plot(finalRho_1_5(:,1), finalRho_1_5(:,2), '-ko', 'MarkerSize', 8, 'LineWidth', 1)
hold on;
plot(finalRho_2_0(:,1), finalRho_2_0(:,2), '-kx', 'MarkerSize', 8, 'LineWidth', 1)
plot(finalRho_2_2(:,1), finalRho_2_2(:,2), '-k^', 'MarkerSize', 8, 'LineWidth', 1)
plot(finalRho_2_4(:,1), finalRho_2_4(:,2), '-k>', 'MarkerSize', 8, 'LineWidth', 1)
set(gca, 'FontSize', 20)
xlabel('P^*')
ylabel('\rho^*')
leg = legend('D=1.50\sigma','D=2.00\sigma','D=2.20\sigma','D=2.40\sigma', 'Location', 'NorthWest');
set(leg, 'FontSize', 14);
legend boxoff;
axes('Position',[.55 .2 .33 .31])
box on
plot(eos2_2(:,1), eos2_2(:,2) / (pi * 1.1^2), '-rv', 'MarkerSize', 5, 'LineWidth', 1);
hold on;
plot(finalRho_2_2(19:44,1), finalRho_2_2(19:44,2), '-k^', 'MarkerSize', 5, 'LineWidth', 1);
plot(eos1_5(:, 2), eos1_5(:, 1), 'rs', 'MarkerSize', 5, 'LineWidth', 1);
plot(finalRho_1_5(19:44,1), finalRho_1_5(19:44,2), '-ko', 'MarkerSize', 5, 'LineWidth', 1);
xlim([7 20]);
set(gca, 'FontSize', 16)
h = figure(1);
saveas(h, 'eos.png');