Set theoretic analysis plays a crucial part in understanding safety requirements for control systems. We apply the lattice theoretic insights of fixed point theorems to construct the maximal invariant set and verify safety constraint satisfaction. We present a Kleene iteration-based algorithm with a novel initialization—based on omega-limit sets—for exactly computing the maximal invariant set under a general safety constraint, and we analyze this algorithm for the special case of Markov chains under polyhedral safety constraints. We also establish sufficient conditions for finite termination of the for a general discrete-time system, which is applicable to Markov chains. In addition to proving that the new initialization of the Kleene iteration algorithm requires no more iterations than the previously considered initialization, we present an example where our method computes the exact maximal invariant set in finitely many iterations and the previous method does not.