% script sin1.m
% Loss of significance example x - sin(x)
% To see the loss of significance show the
% results in binary form (hexadecimal)

f = @(x) x - sin(x);

for x = [0.0001234567890123456, 0.001, 0.1, 1, 1.9, 2.0]
   fprintf('%25.15e %25.15e %s\n', x, f(x), num2hex(f(x)));
end