%% Example problem

%% (a)
% Consider the function f(x) = sin(x)-x. Let x=10^-7 and print out y=f(x) with 15 digits.
%
% How many of these digits are correct???

x = 1e-7;
y = sin(x) - x;
fprintf('y = %.15g\n',y)    % print y with 15 digits

%% (b)
% Plot the function f(x) for x in [0,10^-7].
%
% Why does the plot look like this???

x = linspace(0,1e-7,1000);  % vector of 1000 equidistant points 0,...,1e-7
y = sin(x) - x;             % vector of function values
plot(x,y)                   % plot vector y vs. vector x
xlabel('x')
title('function sin(x)-x')

