%
%  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
%
%  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

p =

  21.8413          
  -0.5210 + 6.0863i
  -0.5210 - 6.0863i
   4.2886          
   3.2137 + 3.3423i
   3.2137 - 3.3423i

z

z =

   0.8496 + 6.9172i
   0.8496 - 6.9172i
  -4.5705          
  -3.3792 + 3.7377i
  -3.3792 - 3.7377i

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