% Sauer 5.1.5
% Kod till övningen 2017-09-28 i SF1514/18/19

clear

f = @(x) 1./x;
fpp = @(x) 2./x.^3; % exakt andraderivata

x = 1;
h = [0.1 0.01 0.001]';

format long
D = (f(x+h) - 2*f(x) + f(x-h))./h.^2 % approximation
E = abs(D-fpp(x)) % fel

loglog(h, E, 'o-')
% Noggrannhetsordningen är lutningen av linjen, dvs 2
