function pplot(x,y)
if abs(imag(y))<1e-8
   if real(y)<0
      str='b.';
   else
      str='r.';
   end
else
   if real(y)<=0;
      str='bo';
   else
      str='ro';
   end
end
plot(x,real(y),str)