clear all
close all
load('im.mat')
figure
plot(a(:,3)./a(:,2),a(:,end-1),'-s','linewidth',1)
hold on
plot(b(:,3)./b(:,2),b(:,end-1),'-d','linewidth',1)
hold on
plot(c(:,3)./c(:,2),c(:,end-1),'-o','linewidth',1)
hold on
plot(d(:,3)./d(:,2),d(:,end-1),'-+','linewidth',1)
% legend('Ca=20','Ca=6.67','Ca=1.25','Ca=0.2','Ca=0.05')
ylim([-0.05,1])
set(gca,'fontsize',12)
set(gcf,'Position',[10 10 300 230])
xlabel('h/w')
ylabel('occupancy percentage')
saveas(gcf,'occu.png')
print(gcf,'-depsc','-painters','occutestnew.eps')