%% Generate an undecorated version of Figure 1(b)

function []=Figure1b()

%% Read in the data
data_name=sprintf('./tau_alphas.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ts=numbers(3*(1:8)-2); % temperatures
tau_alphas=numbers(3*(1:8)-1); % tau_alphas, i.e., relaxation time
tau_alphas_error=numbers(3*(1:8)); % error bars

    
%% Plot them
close all
figure(1);
errorbar(Ts,tau_alphas,tau_alphas_error);

%% Label axis and print it out
xlabel('$T$','Interpreter','latex')
ylabel('$\tau_{\alpha}$','Interpreter','latex')
xlim([0.4,1])
ylim([1,10^4])
set(gca,'yscale','log')
axis square
set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 Figure1b.eps

end