trying to understand why I get 2.2204e-16 rather than 0

12 views (last 30 days)
Im trying to understand why when i make x-y (where x=1.6 and y=1.6) i do not get 0 as a result.

Answers (1)

DGM
DGM on 28 Aug 2021
Edited: DGM on 28 Aug 2021
Simply put, the assumption that both x and y are equal to the same value is incorrect. One or both have accumulated rounding error in some undescribed process.
% x and y are assigned the same value
x = 1.6;
y = 1.6;
x-y
ans = 0
% y has accumulated some rounding error
k = pi;
y = k*(y/k);
x-y
ans = 2.2204e-16
% the difference is nonzero, but they look the same
[x y]
ans = 1×2
1.6000 1.6000
% but that's just because the error is small relative to the nominal value
[x y]-1.6
ans = 1×2
1.0e+-15 * 0 -0.2220
Not all numbers can be represented exactly in base 2, just as not all numbers (e.g. 1/3) can be represented exactly in base 10.
  1 Comment
Walter Roberson
Walter Roberson on 28 Aug 2021
Minor correction: all rational numbers can be represented exactly in base 2 or base 10 -- but it might take infinite precision to do so!! Since there is only a finite amount of resources in the universe, what we can represent in base 2 or base 10 in a computer is finite length rather than infinite length, not all rational numbers can be represented in base 2 or base 10 in a computer .

Sign in to comment.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!