p=0.9*exp(j*pi/10);
a=poly([p conj(p)]);
b=[1 0 0];
theta=0:pi/100:2*pi;
c=1*exp(j*theta);
plot(real(c),imag(c));
set(gca, 'PlotBoxAspectRatio',[1 1 1]);
zvals=roots(b);
pvals=roots(a);
hold on
plot (real(pvals), imag(pvals),'x');
plot (real(zvals), imag(zvals),'o');