%% Generate an undecorated version of Figure 2b

function []=Figure2b()

%% Read in the data
% 2-loop
data_name=sprintf('./2loop_critical.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
nus_2loop=numbers(3*(1:Nd)-1); % critical exponent $\nu$
etas_2loop=numbers(3*(1:Nd)-0); % critical exponent $\eta$

% 3-loop, resummed
data_name=sprintf('./3loop_resummed_critical.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
nus_3loopR=numbers(3*(1:Nd)-1); % critical exponent $\nu$
etas_3loopR=numbers(3*(1:Nd)-0); % critical exponent $\eta$
    
%% Plot them
figure();
hold on
plot(ds2,nus_2loop,'LineStyle','-','Color',0.8*[0,1,1])
plot(ds2,etas_2loop,'LineStyle','-','Color',0.5*[0,0,1])
plot(ds3,nus_3loopR,'LineStyle','--','Color',0.8*[0,1,1])
plot(ds3,etas_3loopR,'LineStyle','--','Color',0.5*[0,0,1])
xlabel('$d$','Interpreter','latex')
ylabel('$\nu,\eta$','Interpreter','latex')
xlim([2,6])
ylim([-0.5,1.0])
axis square
set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 Figure2b.eps
    
end