% make figure 7.8
Z3=zeros(4,491);
t=66;
for i=1:4,
   eval(['load run' num2str(i)]);
   z=sum(real(E).^2,2);
   [y ind]=sort(z);
   newE=E(ind,:);
   Z3(i,:)=sum((real(newE(1:400-t,:)).^2))/(400-t); 
z1(i,:)=sum(real(E).^2,1);
end
figure(3); semilogy(Z3')
hold ontitle('figure 7.8')
print -dpsc 7_8