clear all

filename0 = 'Fig05'  

load Fig05_Pabs den0 den1 den2 PSOL0 PSOL1 PSOL2

load('ehtable','whitejet','greyjet','redblue', 'color7');
set(groot,'defaultAxesColorOrder',color7)

hh=figure(1);
hh.Position = [3 300 420 300];
set(gcf,'color','w');

subplot(111)
plot(den0, PSOL0,'o-', 'color', color7(2,:), 'linewidth', 1.2)
hold on
    plot([2.2e18,2.2e18],[30,100],'k--')
hold off
set(gca,'fontsize', 15)
set(gca, 'ylim', [30,100]);
xlabel('n_{ant} (m^{-3})','fontsize', 17)
ylabel('P_{abs} (%)', 'fontsize', 17)
text(2.3e18, 36, 'n_{FW cut}', 'fontsize', 17)

saveas(gca, filename0,'epsc')
saveas(gca, filename0,'png')
