function []=figureS15a() close all figure() hold on %% Data parameters conversion_factor=1; sample_indices=[0,1,2,3,4]; phis=[0.700, 0.740,0.800,0.820,0.840]; % area fraction Zs=[11.01665,14.5729,24.50337,30.35702,39.50089]; % reduced pressure xi_PTSs_guess=[1.2,1.5,2.3,2.7,3.0]; Ns=size(sample_indices,2); %% for sample_index_dummy=1:Ns sample_index=sample_indices(sample_index_dummy); %% Color colorB=(sample_index_dummy-1)/(Ns-1); colorR=1-colorB; color=[colorR,0,colorB]; %% Loading data % Structural PTS data_name=sprintf('./PTS%d.dat',sample_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all'); NCP=round(size(numbers,1)/3); cavity_sizes=numbers(3*(1:NCP)-2); PTSs=numbers(3*(1:NCP)-1); PTS_errors=numbers(3*(1:NCP)); bx=conversion_factor*cavity_sizes; by=PTSs; be=PTS_errors; % Bulk value data_name=sprintf('./AsymptoticBulk%d.dat',sample_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all'); qbulk=numbers; %% Fitting x=0:0.001:4; weight=diag(ones(size(be))); Outliers=false(size(bx)); by_tilde=by-qbulk; Init=[0.8 xi_PTSs_guess(sample_index_dummy) 2]; [fit_values,confidence_intervals]=stretched_exponential_decay_fit(bx,by_tilde,weight,Init,Outliers); a=fit_values(1); b=fit_values(2); c=fit_values(3); f=qbulk+a*exp(-(x/b).^c); % % % % % % % % xi_fit=b; % % % % % % % % xi_error=(confidence_intervals(2,2)-confidence_intervals(1,2))/2; % % % % % % % % xi_error=1.96*xi_error; %68.2689492137086% confidence * 1.96->95% %% Plotting plot(x,f,'Color',color) errorbar(bx,by,be,'Color',color,'Marker','.','MarkerSize',20,'LineStyle','none') xlim([0,4]) ylim([0,1]) xlabel('$R$','Interpreter','latex'); ylabel('$Q$','Interpreter','latex') end print(gcf, '-dpdf', 'figureS15a.pdf'); end