% make figure 7.9abcd
base=3;
strb=['a'; 'b'; 'c'; 'd'];
for i=1:4,
   figure;
   eval(['load run' num2str(i)])
   z=sum(real(E).^2,2);
   [y ind]=sort(z);
   newE=E(ind,:);
   t=66;
   z_dct=sum((real(newE(1:400-t,:)).^2))/(400-t);
   semilogy(z_dct);
   hold on
   eval(['load ../Chap_5_EQ/run' num2str(i)])
   semilogy(mean(E'.^2),'r')
   eval(['load ../Chap_9/run' num2str(i)])
   semilogy(mean(A.^2),'m')
   hold off
   str=['title(''figure 7.9' strb(i) ''')'];
   eval(str)
   eval(['print -dpsc 7_9' strb(i) ]);
end