%% Generate an undecorated version of Figure 1b

function []=figure1b()
%% Load the data
% N_n=1
x1=10.^((-2):0.01:2).'; % betas
y1=log(2*cosh(x1)); % lambdas

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

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

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

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

y1=-y1./x1;
y2=-y2./x2;
y3=-y3./x3;
y4=-y4./x4;
y5=-y5./x5;
x1=1./x1;
x2=1./x2;
x3=1./x3;
x4=1./x4;
x5=1./x5;

%% Plot
figure();
hold on
plot(x1,y1);
plot(x2,y2);
plot(x3,y3);
plot(x4,y4);
plot(x5,y5);
xlim([0,1])
ylim([-1.2,-0.9])
xlabel('$\tilde{T}$','Interpreter','latex')
ylabel('$\tilde{f}$','Interpreter','latex')

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