%Enter Directory
load('DIRECTORY\input.mat');

%%%%%%%%%%%%%%%%%
ax1 = subplot(1,4,1);
clr = [0.8  0.5  0    ;
        1    0    0  ;
        0.25 0.25 0.9  ;
        1    0.5  0.5  ;
        0    0    0  ];
     colormap(clr);
bar(ax1,log10(input(1:13,:)),'EdgeColor','none','BarWidth', 1)
set(gca,'XTickLabel',{'Austria','Bulgaria','Denmark','France','Germany','Romania','Paraguay','Combination','','','','',''},'fontweight','bold');
set(gca,'XTickLabelRotation',45,'FontSize',9)
title('PC1','FontSize',40)
ylim([0 log10(18160130)])
%%%%%%%%%%%%%%%%%
ax2 = subplot(1,4,2); 
clr = [0.8  0.5  0    ;
        1    0    0  ;
        0.25 0.25 0.9  ;
        1    0.5  0.5  ;
        0    0    0  ];
bar(ax2,log10(input(14:26,:)),'EdgeColor','none','BarWidth', 1)
set(gca,'XTickLabel',{'Denmark','Turkey','USA','Australia','Iran', 'Yemen','Combination','','','','','',''},'fontweight','bold');
set(gca,'XTickLabelRotation',45,'FontSize',9)
title('PC2','FontSize',40)
ylim([0 log10(18160130)])
%%%%%%%%%%%%%%%%%
ax3 = subplot(1,4,3); 
clr = [0.8  0.5  0    ;
        1    0    0  ;
        0.25 0.25 0.9  ;
        1    0.5  0.5  ;
        0    0    0  ];
bar(ax3,log10(input(27:39,:)),'EdgeColor','none','BarWidth', 1)
set(gca,'XTickLabel',{'Turkey','Algeria','Niger','Canada','USA','Iran','Pakistan','Combination','','','','',''},'fontweight','bold');
set(gca,'XTickLabelRotation',45,'FontSize',9)
title('PC3','FontSize',40)
ylim([0 log10(18160130)])
%%%%%%%%%%%%%%%%%
ax4 = subplot(1,4,4);
clr = [0.8  0.5  0    ;
        1    0    0  ;
        0.25 0.25 0.9  ;
        1    0.5  0.5  ;
        0    0    0  ];
bar(ax4,log10(input(40:52,:)),'EdgeColor','none','BarWidth', 1)
set(gca,'XTickLabel',{'France','Germany','Spain','Canada','Yemen','India','Japan','Combination','','','','',''},'fontweight','bold');
set(gca,'XTickLabelRotation',45,'FontSize',9)
title('PC4','FontSize',40)
ylim([0 log10(18160130)])