%% Primera derivada - Central O(h^2)
function df = d1_central(f, x, h)
    df = (f(x+h) - f(x-h)) / (2*h);
end
