www.pudn.com > aft.rar > plot_9_7.m

% remake figure 9.6 
 
figure 
 
load run5 
semilogy(mean(A.^2)'); 
 
eval(['load ../Chap_5_EQ/run15']) 
Z=sum((E.^2),2)/200; 
hold on 
semilogy(Z,'r') 
hold off 
 
set(gca,'xlim',[0 500]) 
set(gca,'ylim',[1e-1 1e1]) 
title('figure 9.7') 
print -dpsc 9_7