%% Generate an undecorated version of Figure 5

function []=Figure5()
%% Data points
liq_FP=[0.694,0; 0.697,0.01573; 0.698,0.03149; 0.716,0.06575; 0.730,0.08450; 0.745,0.09654]; % weakly first-order transition line, (phi_star,Delta)
Deltas=[0,0.03,0.06,0.09,0.11]; % Polydispersity
phi1s=[0.600,0.600,0.620,0.640,0.660]; % Lowest packing fraction data points at each polydispersity 
phi2s=[0.650,0.650,0.670,0.680,0.700];
phi3s=[0.680,0.680,0.690,0.705,0.730];
phi4s=[0.690,0.690,0.700,0.720,0.750];
phi5s=[0.695,0.695,0.710,0.735,0.770]; % Highest packing fraction data points at each polydispersity 

%% Set figure parameters
color1=[1,0,0]; marker1='x'; markersize1=15;
color2=0.5*[0,1,0]; marker2='o'; markersize2=7;
color3=0.7*[0,1,1]; marker3='s'; markersize3=8;
color4=[0,0,1]; marker4='d'; markersize4=8;
color5=[0,0,0]; marker5='+'; markersize5=12;

%% Draw phase diagram
figure()
hFP=line(liq_FP(:,2),liq_FP(:,1));
set(hFP,'Color',[1,0,1]);
hold on
for ii=1:5
    eval(sprintf('h%d=line(Deltas,phi%ds);',ii,ii));
    eval(sprintf('set(h%d,''LineStyle'',''none'',''Marker'',''.'',''Color'',color%d);',ii,ii));
    eval(sprintf('set(h%d,''LineWidth'',1,''Marker'',marker%d,''MarkerSize'',markersize%d,''MarkerFaceColor'',color%d);',ii,ii,ii,ii));
end
xlim([0,0.12])
ylim([0.59,0.8])
xlabel('$\Delta$','Interpreter','latex')
ylabel('$\phi$','Interpreter','latex')
axis square
set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 Figure5.eps
end