%% Generate an undecorated version of the middle row of Figure 7 function []=Figure7middle() polydispersity_indices=[0,3,6,9,11]; % polydispersity, in unit of percentage. for ii=1:size(polydispersity_indices,2) polydispersity_index=polydispersity_indices(ii); % polydispersity, in unit of percentage. figure(); hold on for packing_index=1:5 % packing_index=1 corresponds to packing fraction phi=phi1 (depends on polydispersity); % packing_index=2 corresponds to packing fraction phi=phi2 (depends on polydispersity); % packing_index=3 corresponds to packing fraction phi=phi3 (depends on polydispersity); % packing_index=4 corresponds to packing fraction phi=phi4 (depends on polydispersity); % packing_index=5 corresponds to packing fraction phi=phi5 (depends on polydispersity); %% Read in the data data_name=sprintf('./hexChi_poly%d_phi%d.dat',polydispersity_index,packing_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all'); Nr=round(size(numbers,1)/3);% number of radial data points numbers=reshape(numbers,3,Nr); R=numbers(1,1:Nr).'; % cavity radii where hexatic point-to-set susceptibilities are measured chiPTS=numbers(2,1:Nr).'; % hexatic point-to-set susceptibilities chiPTS_error=numbers(3,1:Nr).'; % error bars %% Plot them errorbar(R,chiPTS,chiPTS_error) end xlabel('$R$','Interpreter','latex') ylabel('$\chi^{\rm PTS}_{6}$','Interpreter','latex'); xlim([0,20]) ylim([0,0.15]) axis square set(gcf, 'PaperPositionMode', 'auto'); % Print out if polydispersity_index==0 print -depsc2 Figure7middle0.eps end if polydispersity_index==3 print -depsc2 Figure7middle3.eps end if polydispersity_index==6 print -depsc2 Figure7middle6.eps end if polydispersity_index==9 print -depsc2 Figure7middle9.eps end if polydispersity_index==11 print -depsc2 Figure7middle11.eps end end end