%% Plot the values of stability exponents at the nontrivial fixed point as a function of spatial dimensions d
function []=figure4a()
%% Load the data
% Stability exponents obtained through minimal 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); 
ds2=numbers(3*(1:Nd)-2); % spatial dimension d
lambda1_loopy=numbers(3*(1:Nd)-1);
lambda2_loopy=numbers(3*(1:Nd)-0);

% Stability exponents obtained through nonperturbative RG
data_name=sprintf('./criticalNRG_stability.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/4);
dsNRG=numbers(4*(1:Nd)-3); % spatial dimension d
lambda1_NRG=numbers(4*(1:Nd)-2);
lambda2_NRG=numbers(4*(1:Nd)-1);
% lambda0_NRG=numbers(4*(1:Nd)-0); discarded from figure 4a, as the same info is contained in figure 4c as ``nu"

%% Plot
close all
figure()
hold on
plot(ds2,lambda1_loopy,'Color',[0,0,0],'LineStyle','-');
plot(ds2,lambda2_loopy,'Color',[0,0,0],'LineStyle','-');
plot(dsNRG,lambda1_NRG,'Color',[0,0,0],'LineStyle','--');
plot(dsNRG,lambda2_NRG,'Color',[0,0,0],'LineStyle','--');
xlabel('$d$','Interpreter','latex')
ylabel('$\lambda$','Interpreter','latex')
xlim([2,6])
ylim([-2,10])
print(gcf, '-dpng', './figure4a.png');
end