gm = [
0 0.4958 0.4958 0.029
0.25 0.4958 0.3414 0.029
0.5 0.4958 0.2118 0.029
0.75 0.4958 0.1604 0.029
0.9 0.4958 0.082 0.029
0.95 0.4958 0.0402 0.029
0.96 0.4958 0.0332 0.029
0.97 0.4958 0.0294 0.029
0.975 0.4958 0.0302 0.029
0.98 0.4958 0.0314 0.029
0.99 0.4958 0.0326 0.029
];

mid=6;
last=11;

hold off
subplot(2,2,1)
semilogy(gm(1:last,1),gm(1:last,2),'k--')
hold on
semilogy(gm(1:last,1),gm(1:last,3),'k-.')
semilogy(gm(1:last,1),gm(1:last,4),'k-')

xlabel 'Model \alpha'
ylabel 'BER'
title '(a) BER vs. Model \alpha'

% axis([0 1 0.5e-1 0.5]);

a = axis;

x = a(1) + 1*(a(2)-a(1))/20;
y = 0.7e-1;

text(x,y, ['---  Gaussian Model     '; ...
           '.-.-. Gauss-Markov Model'; ...
           '\_\_\_\_ Known Channel  ' ]);

hold off
subplot(2,2,3)
semilogy(gm(mid:last,1),gm(mid:last,3),'k-.')
hold on
semilogy(gm(mid:last,1),gm(mid:last,4),'k-')

xlabel 'Model \alpha'
ylabel 'BER'
title '(b) BER vs. Model \alpha'

b = axis;

axis([gm(mid,1) 1 2.5e-2 5e-2])

hold off
