function []=figure2b() close all figure() hold on %% Data parameters ceiling_T_VFT=0.0907; ceiling_T_A=0.0738; T_g=mean([ceiling_T_VFT,ceiling_T_A]); %% Loading data x_grain=0.001*(T_g); x_grain_plot=0.03*(T_g); Ts_fit=0:x_grain:T_g; Ts_fit_plot=0:x_grain_plot:T_g; % PEL 1 data_name=sprintf('./sconf_PEL1.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all'); Ndata=round(size(numbers,1)/2); Ts1=numbers(2*(1:Ndata)-1); Sconfs1=numbers(2*(1:Ndata)); where=logical(Ts1