% Diagramma la capacità di un BEC

clear 
close all
figure

p=[0:.001:1];
for i=1:length(p),
    C(i)=1-p(i);
end

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