allMre0 = [];
for i = 1:20
   allMre0 = [allMre0,mre0(i,:)]; 
end

allMreGac2 = [];
for i = 1:20
   allMreGac2 = [allMreGac2,mreGac2(i,:)]; 
end

allMre1 = [];
for i = 1:20
   allMre1 = [allMre1,mre1(i,:)]; 
end

allMre2 = [];
for i = 1:20
   allMre2 = [allMre2,mre2(i,:)]; 
end

allMre4 = [];
for i = 1:20
   allMre4 = [allMre4,mre4(i,:)]; 
end

allMre8 = [];
for i = 1:20
   allMre8 = [allMre8,mre8(i,:)]; 
end

allMre16 = [];
for i = 1:20
   allMre16 = [allMre16,mre16(i,:)]; 
end

allMrex = [];
for i = 1:20
   allMrex = [allMrex,mrex(i,:)]; 
end

figure;
hold all;
h=gca;
set(h,'FontSize',20);
set(h,'FontName','Times New Roman');

% get size
tempSize = size(allMreGac2,2);
% 
% % get tick label
% yTempTickLabel = sort(log(allMreGac2));
% yTempTickLabel = exp(yTempTickLabel);
% yTickLabelStr = num2str(yTempTickLabel(1), '%10.3f');
% % convert tick labels to string
% for i = 2:tempSize
%     yTickLabelStr = [yTickLabelStr '|' num2str(yTempTickLabel(i), '%10.3f')];
% end
% 
% set(gca,'YTickLabel',yTickLabelStr);

% axis([0 tempSize -10 10])
% plot(1:tempSize,sort(log(allMreGac2)),'LineWidth',1.8);
% plot(1:tempSize,sort(log(allMre1)),'LineWidth',1.8);
% plot(1:tempSize,sort(log(allMre2)),'LineWidth',1.8);
% plot(1:tempSize,sort(log(allMre4)),'LineWidth',1.8);
% plot(1:tempSize,sort(log(allMre8)),'LineWidth',1.8);
% plot(1:tempSize,sort(log(allMre16)),'LineWidth',1.8);
% plot(1:tempSize,sort(log(allMrex)),'LineWidth',1.8);

semilogy(1:tempSize,sort(allMreGac2));
semilogy(,1:tempSizesort(allMre1));
semilogy(,1:tempSizesort(allMre2));
semilogy(,1:tempSizesort(allMre4));
semilogy(,1:tempSizesort(allMre8));
semilogy(sort(allMre16));
semilogy(sort(allMrex));


ylabel('MRE','fontsize',36);
h_legend=legend('TEAK','k=1','k=2','k=4','k=8','k=16','k=best K','Location','NorthWest');
set(h_legend,'FontSize',22);

hold off;