%% Generate an undecorated version of Figure 1d

function []=figure1d()
%% Load the data
% N_n=1
x1=10.^((-2):0.01:2).'; % betas
y1=1./(log(cosh(x1)./sinh(x1))); % xis

% N_n=2
data_name=sprintf('xis%d.dat',2);
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/2);
x2=numbers(2*(1:Ndata)-1); % betas
y2=numbers(2*(1:Ndata)-0); % xis
% Restriction to xi<10^5
x2=x2(1:(Ndata-5)); y2=y2(1:(Ndata-5));

% N_n=3
data_name=sprintf('xis%d.dat',3);
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/2);
x3=numbers(2*(1:Ndata)-1); % betas
y3=numbers(2*(1:Ndata)-0); % xis

% N_n=4
data_name=sprintf('xis%d.dat',4);
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/2);
x4=numbers(2*(1:Ndata)-1); % betas
y4=numbers(2*(1:Ndata)-0); % xis

% N_n=5
data_name=sprintf('xis%d.dat',5);
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/2);
x5=numbers(2*(1:Ndata)-1); % betas
y5=numbers(2*(1:Ndata)-0); % xis

%% Plot
figure();
hold on
plot(x1,y1);
plot(x2,y2);
plot(x3,y3);
plot(x4,y4);
plot(x5,y5);
set(gca,'XScale','log');
set(gca,'YScale','log');
xlim([10^(-1),10^(1)])
ylim([10^(-0),10^(3)])
xlabel('$\tilde{\beta}$','Interpreter','latex')
ylabel('$\xi$','Interpreter','latex')

set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 Figure1d.eps
end