% make figure 7.7

Z2=zeros(4,491);
load run2
tail=[0 5 8 66]; 
z=sum(real(E).^2,2);
[y ind]=sort(z);
newE=E(ind,:);
for i=1:4,
   Z2(i,:)=sum((real(newE(1:400-tail(i),:)).^2))/(400-tail(i));
end
figure(1);semilogy(Z2','k')
title('figure 7.7')
print -dps 7_7