We consider a system of differential equations with nonlinear Steklov
boundary conditions, related to the fractional problem (−Δ)sui=fi(x,ui)−βuipj=i∑aijujp, where i=i,…,k,
s∈(0,1), p>0, aij>0 and β>0. When k=2 we develop a
quasi-optimal regularity theory in C0,α, uniformly w.r.t. β,
for every α<αopt=min(1,2s); moreover we show that the
traces of the limiting profiles as β→+∞ are Lipschitz continuous
and segregated. Such results are extended to the case of k≥3 densities,
with some restrictions on s, p and aij