%% Generate an undecorated version of Figure 2c function []=Figure2c() %% Read in the data % 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); % number of data points ds2=numbers(3*(1:Nd)-2); % spatial dimensions gIs_2loop=numbers(3*(1:Nd)-1); % critical coupling $g^{I}$ gIIs_2loop=numbers(3*(1:Nd)-0); % critical coupling $g^{II}$ % 3-loop, resummed data_name=sprintf('./critical3loop_couplings.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 gIs_3loopR=numbers(3*(1:Nd)-1); % critical coupling $g^{I}$ gIIs_3loopR=numbers(3*(1:Nd)-0); % critical coupling $g^{II}$ %% Plot them figure(); hold on plot(ds2,gIs_2loop,'LineStyle','-','Color',[1,0,0]) plot(ds2,gIIs_2loop,'LineStyle','-','Color',[1,0.7,0]) plot(ds3,gIs_3loopR,'LineStyle','--','Color',[1,0,0]) plot(ds3,gIIs_3loopR,'LineStyle','--','Color',[1,0.7,0]) xlabel('$d$','Interpreter','latex') ylabel('$g_{\star}^{\mathcal{X}}$','Interpreter','latex') xlim([2,6]) ylim([0,2]) axis square set(gcf, 'PaperPositionMode', 'auto'); print -depsc2 Figure2c.eps end