clear all
close all
load('cycle.mat')
figure
plot(a(:,3)./a(:,2),a(:,end-1),'-s','linewidth',2,'markersize',10)
hold on
plot(b(:,3)./b(:,2),b(:,end-1),'-*','linewidth',2,'markersize',10)
hold on
plot(c(:,3)./c(:,2),c(:,end-1),'-d','linewidth',2,'markersize',10)
ylim([-0.05,1])
legend('Ca/Re=5','Ca/Re=0.2','Ca/Re=0.05')
set(gca,'fontsize',20)
saveas(gcf,'occu.png')
print(gcf,'-depsc','-painters','occuimdr.eps')


figure
plot(a(:,3)./a(:,2),a(:,end-1)./drp,'-s','linewidth',2,'markersize',10)
hold on
plot(b(:,3)./b(:,2),b(:,end-1)./drp,'-*','linewidth',2,'markersize',10)
hold on
plot(c(:,3)./c(:,2),c(:,end-1)./drp,'-d','linewidth',2,'markersize',10)
ylim([-0.05,1])
legend('Ca/Re=5','Ca/Re=0.2','Ca/Re=0.05')
set(gca,'fontsize',20)
saveas(gcf,'occu1.png')
print(gcf,'-depsc','-painters','occuimbydr.eps')