%machepstd.m
x=1;
y=1;
z=x+y;
n=0;
data=[0 1 2];
while z>x
   n=n+1;
   y=y/2;
   z=x+y;
   data=[data;n y z];
end
