hold off
load gm.32.txt
G=gm
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--')
semilogy(G(:,2),G(:,7),'k+--')
% axis tight
% axis([0 20 1e-3 1e0]);
axis([0 20 2e-2 2e-1]);
xlabel 'SNR (dB)'
ylabel 'BER'
title 'gm.32: \alpha = 0.08, var\_h = 0.05'
hold off

a = axis;

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

x = 22;
y = 0.5e-1;

text(x,y, ['.......  deterministic EM'; ...
           '-.-.-. Gauss-Markov EM   '; ...
           '---- Gauss-Markov ML     '; ...
           '-+-+ Gaussian ML         '; ...
           '\_\_\_\_ known channel   ']);

print -deps gm.32.eps

%	    1		 2	      3	           4		5	     6         7
%
%           13            5      0.15367      0.15867     0.092667        0.154     0.159
%           13           10         0.09        0.094    0.0086667        0.073     0.094
%           13           15     0.035667        0.041   0.00033333        0.024     0.041
%
%           N, 		SNR, 	    det,           g,       known,       ML GM,    ML G
