%% Generate an undecorated version of Figure 2a

function []=Figure2a()

%% Read in the data
% 2-loop
data_name=sprintf('./critical2loop_stability.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3); % number of data points
ds2=numbers(3*(1:Nd)-2); % spatial dimensions
lambda1_2loop=numbers(3*(1:Nd)-1); % 1st stability exponent
lambda2_2loop=numbers(3*(1:Nd)-0); % 2nd stability exponent

% 3-loop, resummed
data_name=sprintf('./critical3loop_stability.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3); % number of data points
ds3=numbers(3*(1:Nd)-2); % spatial dimensions
lambda1_3loopR=numbers(3*(1:Nd)-1); % 1st stability exponent
lambda2_3loopR=numbers(3*(1:Nd)-0); % 2nd stability exponent
    
%% Plot them
figure();
hold on
plot(ds2,lambda1_2loop,'LineStyle','-','Color',[0,0,0])
plot(ds2,lambda2_2loop,'LineStyle','-','Color',[0,0,0])
plot(ds3,lambda1_3loopR,'LineStyle','--','Color',[0,0,0])
plot(ds3,lambda2_3loopR,'LineStyle','--','Color',[0,0,0])
xlabel('$d$','Interpreter','latex')
ylabel('$\lambda$','Interpreter','latex')
xlim([2,6])
ylim([-2,10])
axis square
set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 Figure2a.eps
    
end