X(1)=1/9;
for k=1:70
  X(k+1)=doubling(X(k));
end
plot(X,'.')