An effective method is presented to compute the loop gain margins exactly for two-input–two-output processes. The stable regions in parameter space are first obtained by determining the stability boundaries and the loop gain margins found within the stable regions. Examples are provided for illustration and comparison with other methods.