function m = getSlope(f)

    eps = 1e-9;
    % finite difference
    m = @(x0) (f(x0 + eps) - f(x0)) / eps;

end