x=(0:20)';
px = x(find(x>=1));
clf;
plot(px,ternary(px),'--',px,binary(px),'-');
hold on;
plot(x,ones(size(x))*(log(2)/log(3)),':');
legend('ternary supercells','binary supercells',4);
title('information density');
set(gca,'ylim',[0 1]);

l23 = log(2)/log(3);
text(0.5,l23,'log_{2}3');

% ytick = get(gca,'ytick');
% near=abs(ytick-l23);
% [y,i] = min(near);
% ytick(i) = l23;
% yticklabel = get(gca,'yticklabel');
% yticklabel = strvcat(yticklabel(1:i-1,:),'log_{2}3',yticklabel(i+1:end,:));
% set(gca,'ytick',ytick);
% set(gca,'yticklabel',yticklabel);

cross = fminbnd('minfun',4,8);
plot([cross;cross],[0;1],':');
text(cross,0.95,sprintf('%4.1f', cross));

print('-deps', 'compare.eps');