%remake plot 5.18

figure
load run5
mu = rp.mu; a  = rp.a; sigu2 = 0.995;  sigv2 = (1-(a^2)) * sigu2; n = 1:5000; 
a2=a^2; p5mu=0.5*mu;
j = sigu2*(1-a2)*(1+p5mu*sigu2) + sigu2*(a2+p5mu*a2*sigu2-0.5*mu*sigu2)*(1-mu*sigu2).^(2*n);
semilogy(j,'--'); hold on
semilogy(mean((E').^2));semilogy(j,'-');
%legend('Theory', 'Experiment')
title('figure 5.18 \mu=0.001')
print -dps 5_18