English
- %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
