figure(1)

for i=1:4,
	eval(['load run' num2str(i)])
	P(1:488,i)=mean(A.^2)';
end

semilogy(P)
set(gca,'xlim',[1 200])
title('figure 13.6')
print -dpsc 13_6

figure(2)

load run5
semilogy(mean(A.^2)');

eval(['load ../9_7_lms/run15'])
Z=sum((E.^2),2)/200;
hold on
semilogy(Z)
hold off

set(gca,'xlim',[0 500])
set(gca,'ylim',[1e-1 1e1])
title('figure 13.7')
print -dpsc 13_7