
% Plot reference lines
x=[1:100]/100;
H1=log2((x).^(-x))+log2(((1-x)).^(-(1-x)));
plot(x,H1,'LineStyle',':',...
    'Color',[0.501960813999176 0.501960813999176 0.501960813999176]);
H2=log2((x).^(-x))+2*log2((((1-x)/2)).^(-(1-x)/2));
hold on
plot(x,H2,'LineStyle',':',...
    'Color',[0.501960813999176 0.501960813999176 0.501960813999176]);
hold off
H3=log2((x).^(-x))+3*log2((((1-x)/3)).^(-(1-x)/3));
hold on
plot(x,H3,'LineStyle',':',...
    'Color',[0.501960813999176 0.501960813999176 0.501960813999176]);
hold off

