This is a multi-species version of the classic Lotka-Volterra competition model. The dynamics of each species is given by:
dBi/dt = ri*(1-sumj = 1,nsp(ci,j*Bj))
where:
Bi is the biomas of species i;
ri is the maximum specific growth rate for species i;