s=poly(0,"s");
p2a=(2*%i)^(1/2), p2a2 = exp(%i*%pi)*p2a
//alternately
p=s^2-2*%i;
roots(p)


p2b=(1-%i*sqrt(3))^(1/2), p2b2 = exp(%i*%pi)*p2b
// alternately
p=s^2-1+%i*sqrt(3);
roots(p)

p3a = (-1)^(1/3), p3a2=exp(%i*2*%pi/3)*p3a, p3a3=exp(%i*2*%pi/3)*p3a2
// alternately
p=s^3+1;
roots(p)

p3b = 8^(1/6), for j=1:5, exp(%i*%pi*j/3)*p3b, end
// alternately
p=s^6-8;
roots(p)

s=poly(0,"s");
p=4+s^4
roots(p)
