%% Plot the values of critical couplings at the nontrivial fixed point as a function of spatial dimensions d
function []=figure4c()
%% Load the data
% Critical couplings obtained through minimal 2-loop
data_name=sprintf('./critical2loop_couplings.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
gI_stars2=numbers(3*(1:Nd)-1);
gII_stars2=numbers(3*(1:Nd)-0);

% Critical couplings obtained through nonperturbative RG
data_name=sprintf('./criticalNRG_couplings.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);
dsN=numbers(3*(1:Nd)-2); % spatial dimension d
gI_starsN=numbers(3*(1:Nd)-1);
gII_starsN=numbers(3*(1:Nd)-0);

%% Plot
close all
figure()
hold on
plot(ds2,gI_stars2,'Color',[1,0,0],'LineStyle','-');
plot(ds2,gII_stars2,'Color',[1,0.7,0],'LineStyle','-');
plot(dsN,gI_starsN,'Color',[1,0,0],'LineStyle','--');
plot(dsN,gII_starsN,'Color',[1,0.7,0],'LineStyle','--');
xlabel('$d$','Interpreter','latex')
ylabel('$g_{\star}^{\mathcal{X}}$','Interpreter','latex')
xlim([2,6])
ylim([0,2])
print(gcf, '-dpng', './figure4c.png');
end