% [h0_var Gauss GM Known]
gm = [
0.0001 0.5018 0.1272 0.0248
0.001 0.5018 0.1256 0.0248
0.005 0.5018 0.1634 0.0248
0.01 0.5018 0.215 0.0248
0.05 0.5018 0.3356 0.0248
0.1 0.5018 0.3804 0.0248
];

last=size(gm,1);

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

xlabel 'Model h_0 variance'
ylabel 'BER'
title 'BER vs. Model h_0 variance'

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

a = axis;

x = 2e-4;
y = 0.6e-1;

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