% Diagramma del piano di Shannon

clear 
close all
figure
E1=[-10:40];
p=[0:.005:8];
p1=ones(length(E1),1);
for i=2:length(p),
   E(i)=((2^p(i)-1)/p(i));
end
E(1)=log(2);
plot(E,p,'-r','linewidth',1.3);
hold
plot(E1,p1,'-b','linewidth',1.3);
axis([-1 32 0 10])
xlabel('E_b/N_0')
ylabel('R_b/W')
grid