% Diagramma la capacità di un BSC

clear 
close all
figure

p=[0:.001:1];
for i=2:length(p)-1,
    H(i)=-p(i)*log2(p(i))-(1-p(i))*log2(1-p(i));
end
H(1)=0;
H(length(p))=0;
C=ones(1,length(H))-H;

plot(p,C,'-r','linewidth',1.2);
axis([0 1 0 1.2])
xlabel('p')
ylabel('C_{BSC}    [bit/simbolo]')
grid