We propose, analyze and test two novel fully discrete decoupled linearized algorithms for a nonlinearly coupled reaction-diffusion N -species competition model with harvesting or stocking effort. The time-stepping algorithms are first and second order accurate in time and optimally accurate in space. Stability and optimal convergence theorems of the decoupled schemes are proven rigorously. We verify the predicted convergence rates of our analysis and efficacy of the algorithms using numerical experiments and synthetic data for analytical test problems. We also study the effect of harvesting or stocking and diffusion parameters on the evolution of species population density numerically, and observe the co-existence scenario subject to optimal harvesting or stocking.