%  This routine plots the order star accuracy domain of BI4/5_0.45
%  ---------------------------------------------------------------
%
%  Find order star accuracy BI4/5_0.45
%
algor = 1;
j = sqrt(-1);
%
%  Sweep through s
%
errvec = zeros(121,121);
s = [-3:0.05:3];
w = [-3:0.05:3];
for i=1:121,
  ss = s(i);
  %
  %  Sweep through w
  %
  for k=1:121,
    ww = w(k);
    lambd = ss + j*ww;
    shat = -damp(-lambd,algor);
    ohat = freq(lambd,algor);
    errvec(k,i) = abs(ss - shat) + abs(ww - ohat);
  end,
end
%
%  Plot the results
%
tol = 1.0e-4;
contour(s,w,errvec,[tol tol],'k-')
grid on
axis equal
title('ttt')
xlabel('xxx')
ylabel('yyy')
hold on
tol = 1.0e-3;
contour(s,w,errvec,[tol tol],'k-')
tol = 1.0e-2;
contour(s,w,errvec,[tol tol],'k-')
text(-2.6,0.7,'tt2')
text(-1.6,1.3,'tt3')
text(-0.7,-0.8,'tt4')
%
print -deps nsds_hw3_19c.eps
%
return
