%% Primera derivada - Central O(h^4) mejorada
function df = d1_central4(f, x, h)
    df = (-f(x+2*h) + 8*f(x+h) - 8*f(x-h) + f(x-2*h)) / (12*h);
end
