%% Segunda derivada - Central O(h^4) mejorada
function d2 = d2_central4(f, x, h)
    d2 = (-f(x+2*h) + 16*f(x+h) - 30*f(x) + ...
           16*f(x-h) - f(x-2*h)) / (12*h^2);
end
