% Plotting
BERpersub1 = squeeze(sum(BER1,2)./sum(N1,2));
BERpersub2 = squeeze(sum(BER2,2)./sum(N2,2));
BERpersub3 = squeeze(sum(BER3,2)./sum(N3,2));
for s = 1 : 3
  figure(1)
  plot(C(1,:,s),C(2,:,s),'r-x',C2(1,:,s),C2(2,:,s),'go--',C3(1,:,s),C3(2,:,s),'m^-.')
  hold on
  grid on
  title('Capacity Region for 2 Users')
  xlabel('Capacity for user 1 (bps/Hz)')
  ylabel('Capacity for user 2 (bps/Hz)')
end
legend('Perfect CSI','Imperfect CSI - Optimal','Imperfect CSI - Suboptimal')

BERpersub1 = squeeze(sum(BER1,2)./sum(N1,2));
BERpersub2 = squeeze(sum(BER2,2)./sum(N2,2));
BERpersub3 = squeeze(sum(BER3,2)./sum(N3,2));