An algorithm for estimating the stability domain of linear discrete-time systems with saturations