%% Segunda derivada - Central O(h^2)
function d2 = d2_central(f, x, h)
    d2 = (f(x+h) - 2*f(x) + f(x-h)) / h^2;
end
