% temp good, replace with real one later

good = [
           10            0            0        0.287      0.23657      0.24557
           10            2            0      0.23557      0.18243        0.199
           10            4            0      0.19843        0.132        0.151
           10            6            0      0.17014        0.086      0.10214
           10            8            0        0.139        0.036        0.058
           10           10            0        0.112     0.010143     0.025571
];

bad = [
           10            0            0      0.67029      0.23657            0
           10            2            0      0.68957      0.18243            0
           10            4            0        0.701        0.132            0
           10            6            0        0.701        0.086            0
           10            8            0      0.70057        0.036            0
           10           10            0      0.70286     0.010143            0
];

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, ['.......... EM init #1   '; ...
           '.-.-.- EM init #2       '; ...
           '----  true MAP          '; ...
           '\_\_\_\_\_ known channel']);

