function small = epsilon()

small = 1

while ( 1+ small/2> 1 )
small = small/2;

end
