So just thinking out loud here but like in the 1d case... The set will be an interval [a,b] Let g(x) = f(x) - x
g(a) >= 0 and g(b) <= 0. Now we can use the intermediate value value theorem and binary search to find (get arbitrary close) a point with g(x) = 0 (i.e. fix point of f).
Not possible to do such searches in higher dimensions?