Your original number has too many digits to store in a double precision number. E.g.,
num2strexact(73167176531330624919225119674426574742355349194934)
ans =
7.3167176531330625585928818171611107149744948903936e49
So what you are seeing as output is exactly what is being stored (i.e., the exact decimal conversion of the binary bit pattern being stored), and there aren't enough mantissa bits in double precision to store more precision.
In fact, the nearest numbers in the double precision set to your original number are:
>> x = 73167176531330624919225119674426574742355349194934
>> num2strexact(x)
ans =
7.3167176531330625585928818171611107149744948903936e49
>> num2strexact(x+eps(x))
ans =
7.3167176531330635970522535241266364210737607344128e49
>> num2strexact(x-eps(x))
ans =
7.3167176531330615201335101101955850088752290463744e49
Best Answer