% Figure 4.3
eps=[1 2 4]; sig=[1 2 3];
e=linspace(0.5,1,101);
sigmodel=e'*eps;
sigr=ones(101,1)*sig;
diff=abs(sigr-sigmodel); maxdiff=max(diff');
sumsquares=diag(diff*diff'); rms=sqrt(sumsquares/3);
plot(e,maxdiff ); xlabel('E'); ylabel('cost function')
hold on; plot(e,rms,'r-'); legend('max error','RMS
error','Location','North')
%-----------------------------------------------------------------------------