x=linspace(-10,10,100);
y=roots(x.^3-1);
r=real(y);
i=imag(y);
plot(r,i,'x');