function [s]=pl2sp(c);
th=angle(c);r=abs(c);
x=r.*cos(th)./(r.^2+1);
y=r.*sin(th)./(r.^2+1);
z=((r.^2-1)./(r.^2+1)+1)/2;
s=[x y z];