u=1;
while (u+1 ~= 1) u=u*0.5; end
u=u*2
dig=log(u)/log(2)
