load ml.01.txt
G=ml;
hold off
subplot(2,2,1)
semilogy(G(:,2),G(:,3),'k:')
hold on
semilogy(G(:,2),G(:,4),'k-.')
semilogy(G(:,2),G(:,5),'k-')
semilogy(G(:,2),G(:,6),'k--')
xlabel 'SNR (dB)'
ylabel 'BER'
title '(a) ml.01: reasonable init'

axis([0 20 1e-3 1e0]);

xlabel 'SNR (dB)'
ylabel 'BER'

x = 1;
y = 0.5e-2;

text(x,y, ['.......  deterministic h'; ...
           '-.-.-. Gaussian h       '; ...
           '---- direct ML          '; ...
           '\_\_\_\_ known channel  ']);

load ml.02.txt
G=ml;
hold off
subplot(2,2,2)
semilogy(G(:,2),G(:,3),'k:')
hold on
semilogy(G(:,2),G(:,4),'k-.')
semilogy(G(:,2),G(:,5),'k-')
semilogy(G(:,2),G(:,6),'k--')
xlabel 'SNR (dB)'
ylabel 'BER'
title '(b) ml.02: bad init'

axis([0 20 1e-3 1e0]);

xlabel 'SNR (dB)'
ylabel 'BER'

print -deps bad.eps
