function [] = testder_x3(x)
% si "ripulisce" il workspace e lo schermo
clear;
clc;

%x = -0.5;

NP = 5;

ek = 0:NP;

nk = 10.^ek;

for ik = 1:length(nk)
    h = nk(ik);
    Df(ik)=(fracx3(x+h)-fracx3(x-h))./(2*h);
end

Df0 = derfracx3(x);

Ek = abs(Df - (Df0));

Erk = Ek/Df0;


barra = '----------------------------------------------------------------';
titolo = 'Tabella degli errori nel calcolo approssimato della derivata';

disp(titolo);
disp(barra);
fprintf('  ek  \t     n  \t  Df0 \t\t |E_k| \t\t E_r \n');
disp(barra);
fprintf('  %d \t %7d \t %8.7g \t %4.3e  \t %4.3g     \n',...
    [ek; nk; Df; Ek; Erk]);

%%%%%%%%%%%%%%%%%%%%%%%%
function y=fracx3(x)

y = 1 ./ (1 + x.^3);

%%%%%%%%%%%%%%%%%%%%%%%%
function yd = derfracx3(x)

yd = - (3*x.^2) ./ (1+x.^3).^2;
