clear all
close all
load('dr.mat')
figure
plot(c(:,3)./c(:,2),c(:,end),'-*','linewidth',2,'markersize',10)
hold on
plot(b(:,3)./b(:,2),b(:,end),'-d','linewidth',2,'markersize',10)
hold on
plot(a(:,3)./a(:,2),a(:,end),'-s','linewidth',2,'markersize',10)
legend('Ca=500','Ca=20','Ca=0.2')
set(gca,'fontsize',20)
xlabel('h/w')
ylabel('occupancy percentage')
saveas(gcf,'occudr.png')
print(gcf,'-depsc','-painters','occudr.eps')

figure
plot(c(:,3)./c(:,2),c(:,end).*c(:,2).*c(:,3),'-*','linewidth',2,'markersize',10)
hold on
plot(b(:,3)./b(:,2),b(:,end).*b(:,2).*b(:,3),'-d','linewidth',2,'markersize',10)
hold on
plot(a(:,3)./a(:,2),a(:,end).*a(:,2).*a(:,3),'-s','linewidth',2,'markersize',10)
legend('Ca=500','Ca=20','Ca=0.2')
set(gca,'fontsize',20)
xlabel('h/w')
ylabel('occupancy volume')
saveas(gcf,'occudrv.png')
print(gcf,'-depsc','-painters','occudrv.eps')