We study the probability of stability of a large complex system of size N within the framework of a generalized May model, which assumes a linear dynamics of each population size ni (with respect to its equilibrium value):The ai > 0's are the intrinsic decay rates, Jij is a real symmetric (N × N ) Gaussian random matrix and √ T measures the strength of pairwise interaction between different species. Unlike in May's original homogeneous model, each species has now an intrinsic damping ai that may differ from one another. As the interaction strength T increases, the system undergoes a phase transition from a stable phase to an unstable phase at a critical value T = Tc. We reinterpret the probability of stability in terms of the hitting time of the level b = 0 of an associated Dyson Brownian Motion (DBM), starting at the initial position ai and evolving in 'time' T . In the large N → ∞ limit, using this DBM picture, we are able to completely characterize Tc for arbitrary density µ(a) of the ai's. For a specific flat configuration ai = 1 + σ i−1 N , we obtain an explicit parametric solution for the limiting (as N → ∞) spectral density for arbitrary T and σ. For finite but large N , we also compute the large deviation properties of the probability of stability on the stable side T < Tc using a Coulomb gas representation.