- figure
- for i=1:4,
- eval(['load run1', num2str(i)])
- E_priori=E./Gamma(:,6:501);
- semilogy(mean(E_priori.^2))
- hold on
- end
- axis([0 200 1e-3 1e1])
- xlabel('Number of iterations, n')
- ylabel('Mean- squared error')
- gtext('W = 2.9')
- gtext('W = 3.1')
- gtext('W = 3.3')
- gtext('W = 3.5')
- print -deps 12_1