% remake figure 5.17

figure
load run5;plot(mean(Wx'), '--'); hold on

mu=rp.mu; n=1:5000;
a=-0.99; lambda=abs(a); Ew1=-a+a*(1-mu*lambda).^n;
a=0.99;  lambda=abs(a); Ew2=-a+a*(1-mu*lambda).^n;
plot(Ew1); plot(Ew2)
load run5b;plot(mean(Wx'),'--');
line([1 5000], [0 0]);hold off
text(250,0.62,'a = - 0.99'); text(250,-0.62,'a = + 0.99')
title(' figure 5.17 \mu=0.001')
     legend('Experiment', 'Theory')
print -dpsc 5_17