We analyse the dynamics of a prototype model for competing species with diffusion coefficients (d1d2) in a heterogeneous environment Ω. When diffusion is switched off, at each pointx∊ Ω we have a pair of ODE's: thekinetic. If for somex∊ Ω kinetic has a unique stable coexistence state, we show that there existsuch that for everythe RD-model ispersistent, in the sense that it has a compact global attractor within the interior of the positive cone and has a stable coexistence state. The same result is true if there existxu,xv∊ Ω such that the semitrivial coexistence states (u, 0) and (0,v) of thekineticare globally asymptotically stable atx=xuandx=xv, respectively. More generally, our main result shows that, for most kinetic patterns, stable coexistence of xspopulations can be found for some range of the diffusion coefficients.Singular perturbation techniques, monotone schemes, fixed point index, global analysis ofpersistence curves, global continuation and singularity theory are some of the technical tools employed to get the previous results, among others. These techniques give us necessary and/or sufficient conditions for the existence and uniqueness of coexistence states, conditions which can be explicitly evaluated by estimating some principal eigenvalues of certain elliptic operators whose coefficients are solutions of semilinear boundary value problems.We also discuss counterexamples to the necessity of the sufficient conditions through the analysis of the local bifurcations from the semitrivial coexistence states at the principal eigenvalues. An easy consequence of our analysis is the existence of models having exactly two coexistence states, one of them stable and the other one unstable. We find that there are also cases for which the model hasthree or morecoexistence states.