
good = [
           10            0            0      0.23357        0.222      0.23443
           10          2.5            0        0.202      0.18229        0.205
           10            5            0      0.12586      0.10271        0.125
           10          7.5            0        0.078     0.038571     0.069857
           10           10            0     0.040286     0.010857     0.034429
];

hold off
subplot(2,2,1)
semilogy(good(:,2),good(:,5),'k-')
hold on
semilogy(good(:,2),good(:,6),'k--')
semilogy(good(:,2),good(:,4),'k-.')
% semilogy(good(:,2),bad(:,4),'k:')
xlabel 'SNR (dB)'
ylabel 'BER'
% title 'BER vs. SNR'
hold off

% a = axis;

x = 0.5; % 11; % a(1) + (a(2)-a(1))/20;
y = 0.3e-1; % 1.2e-1; % a(3) + (a(4)-a(3))/20;

text(x,y, ['.-.-.- good EM init     '; ...
           '----  true MAP          '; ...
           '\_\_\_\_\_ known channel']);

