%  This routine plots the damping order star of BI4/5_0.45
%  -------------------------------------------------------
%
echo on
diary nsds_hw3_19a.dry
%
%  Find order star for BI4/5_0.45
%
algor = 1;
j = sqrt(-1);
%
%  Sweep through s
%
svec = zeros(161,201);
s = [-10:0.1:10];
w = [-8:0.1:8];
echo off
for i=1:201,
  ss = s(i);
  %
  %  Sweep through w
  %
  for k=1:161,
    ww = w(k);
    lambd = ss + j*ww;
    shat = -damp(-lambd,algor);
    svec(k,i) = ss - shat;
  end,
end
echo on
%
%  Plot the results
%
contour(s,w,svec,[0 0],'k-')
grid on
axis equal
title('ttt')
xlabel('xxx')
ylabel('yyy')
[p,z] = pz(1);
p
z
text(real(p(2)),imag(p(2)),'x')
text(real(p(3)),imag(p(3)),'x')
text(real(p(4)),imag(p(4)),'x')
text(real(p(5)),imag(p(5)),'x')
text(real(p(6)),imag(p(6)),'x')
text(real(z(1)),imag(z(1)),'o')
text(real(z(2)),imag(z(2)),'o')
text(real(z(3)),imag(z(3)),'o')
text(real(z(4)),imag(z(4)),'o')
text(real(z(5)),imag(z(5)),'o')
text(-7.5,5.2,'neg')
text(-7.5,-5.2,'neg')
text(6.8,5.2,'pos')
text(6.8,-5.2,'pos')
text(-3.0,3.0,'pos')
text(-3.0,0.0,'pos')
text(-3.0,-3.0,'pos')
text(1.9,2.7,'neg')
text(2.2,0.0,'neg')
text(1.9,-2.7,'neg')
%
print -deps nsds_hw3_19a.eps
%
%  Close diary and exit
%
diary off
echo off
return
