The system of differential equations is given below: dx_i/dt = omega - mu * x_i + \\epsilon * x_i * y_i dy_i/dt = r * y_i * (1-y_i/k) - gamma * x_i * y_i + sum_{