%% Generate an undecorated version of Figure S7C function []=figureS7C() volume_fraction_index=5; % corresponds to volume fraction phi=0.618 cavity_radius_indices=[45,55,7,110]; % corresponds to cavity radii R=[2.25,2.59,3.11,4.15] for ii=1:4 %% Read in the data cavity_radius_index=cavity_radius_indices(ii); data_name=sprintf('./PDF%d_%d.dat',volume_fraction_index,cavity_radius_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all'); Ng=round(size(numbers,1)/2); % number of grid points for the histogram q_c=numbers(2*(1:Ng)-1); % core overlap P=numbers(2*(1:Ng)-0); % probability distribution function %% Plot them plot(q_c,P) hold on end xlim([0,1]) ylim([0,6]) xlabel('$q_{\rm c}$','Interpreter','latex') ylabel('$P$','Interpreter','latex') axis square set(gcf, 'PaperPositionMode', 'auto'); print -depsc2 FigureS7C.eps end