function R=cantor(I)
%Cantor set.  -JJL
x1=(1/3)*I;
x2=x1+(2/3);
x=[x1,x2];
plot(x,zeros(size(x)),'.')
if nargout>0,R=[x1,x2];end