Operating criteria for power systems, such as the (N − 1)-criterion, are often based on evaluating whether the system is vulnerable to a specific set of contingencies. Therefore, a major part of power system security is concerned with establishing regions in parameter space where the system is vulnerable to specific contingencies. In this article we exploit the possibility of using Monte Carlo simulations to build an approximation of the region, in parameter space, where the power system will remain stable following a given contingency.