pos=0;
while pos<1
    step=(1-pos)/2;
    pos=pos+step;
    disp(pos)
end
