function y=napaka(x,h)

% Numericno racunamo prvi odvod za 1/(1+x) po simetricni 
% pettockovni formuli, vse vrednosti funkcije zaokrozimo
% na 6 decimalk (abs. napaka je pod 5.10^-7).

y0=1/(1+x-2*h);
y0=round(1e6*y0)/1e6;
y1=1/(1+x-h);
y1=round(1e6*y1)/1e6;
y3=1/(1+x+h);
y3=round(1e6*y3)/1e6;
y4=1/(1+x+2*h);
y4=round(1e6*y4)/1e6;

apr=(y0-8*y1+8*y3-y4)/(12*h); % aproksimacija
toc=-1/(1+x)^2;               % tocna vrednost odvoda 
y=abs(apr-toc);               % absolutna napaka