Mathematical Problems in Engineering

Volume 2016, Article ID 7920482, 7 pages

http://dx.doi.org/10.1155/2016/7920482

## A Novel Method to Determine the Local Stability of the n-Species Lotka-Volterra System with Multiple Delays

^{1}College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China^{2}Department of Mathematics, Taizhou University, Taizhou 225300, China

Received 18 May 2016; Revised 4 November 2016; Accepted 20 November 2016

Academic Editor: Leonid Shaikhet

Copyright © 2016 Xiao-Ping Chen and Hao Liu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The n-species Lotka-Volterra system with discrete delays is considered. The local asymptotic stability of positive equilibrium is investigated based on a contour integral method. The main purpose of this paper is to propose a new and general algorithm to study the local asymptotic stability of the positive equilibrium for the -dimensional Lotka-Volterra system. Some numerical experiments are carried out to show the effectiveness of the proposed method.

#### 1. Introduction

In recent years, delay differential equations have been extensively used as models in population dynamics. In this paper, we study the local stability of the positive equilibrium of the following n-species Lotka-Volterra competition system with multiple delays: where denotes the population size of the th species at time and , , and are positive constants.

As a rudimentary model on mathematical ecology, Lotka-Volterra system has received a lot of attention and has been researched extensively; see, for example, [1–15] and the references cited therein. A basic model is the two-species competitive system with multiple delays which takes the form of Zhen and Ma [8] studied the local and global stability of the positive equilibrium of (2) by means of Lyapunov functions. There are also some other researches regarding this area; see, for example, [9–11]. Besides, there are some literatures on three-species models. Qin [12] first investigated the stability of positive equilibrium of a three-species Lotka-Volterra model with discrete delays. Ma et al. [13] studied the local stability of positive equilibrium of a three-species Lotka-Volterra model with discrete delays through constructing a proper Lyapunov function. For the general n-species models, Hofbauer and So [14] gave necessary and sufficient conditions for the global asymptotic stability of a positive equilibrium with the restriction . Later on, Campbell [15] extended the results of [14] and studied the local and global stability of the trivial equilibrium of a functional differential equation used to model artificial neural networks with discrete delays. However, in these literatures, we emphasize that the usual approach to study the stability of the equilibrium for system (1) relies on the use of Lyapunov functionals. Generally speaking, it is not easy to construct a Lyapunov functional for an n-dimensional Lotka-Volterra model. Besides, a new Lyapunov functional is required for each model under consideration. As far as we know, there is no good all-purpose algorithm to study the stability of n-species Lotka-Volterra competition system with multiple delays at present. Therefore, it is necessary to find a general algorithm to study the local stability of the positive equilibrium of system (1) which owns more adaptability. In this paper, we propose a novel and general method to study the local asymptotic stability of positive equilibrium of system (1) by the contour integral method. To the best of our knowledge, this paper is the first to deal with the research of analyzing the local stability of Lotka-Volterra system by the contour integral method.

The paper is organized as follows. In Section 2, we recall an algorithm to find the rightmost eigenvalues of the nonlinear eigenvalue problem (NEP) in some disk and propose a novel algorithm for determining the local asymptotic stability of the positive equilibrium for the general n-species Lotka-Volterra system. Section 3 is devoted to some numerical experiments. Finally, some concluding remarks are made in Section 4.

For convenience, we use the following notations: and denote the trace and the determinant of a matrix , respectively; denotes -norm of a matrix ; denotes the identity matrix with suitable size; is the th column vector of ; denotes the transpose of a vector or a matrix; is the boundary of an open disk and its radius is ; denotes the imaginary unit of a complex number.

#### 2. Local Asymptotical Stability of the Positive Equilibrium

For convenience, in this paper, we impose the following hypothesis for the n-species Lotka-Volterra competition system (1): and , , where .

It is well known (see [16]) that if condition is satisfied, then system (1) has a unique positive equilibrium and we can compute the unique positive equilibrium by the following algebraic equation:

Let then system (1) reduces to the following system:

According to expression (5), the variational system of (1) with respect to the positive equilibrium is given by the linearized system of (5). That is,

Rewrite (6) in the matrix form as follows: where and .

It is easy to see that the characteristic matrix of system (7) is

To establish an effective way to determine whether or not the positive equilibrium for system (1) is asymptotically stable, it suffices to study the asymptotic stability of system (7). It is well known that the local stability of trivial solution of system (7) depends on the locations on the complex plane of the characteristic roots of . And these characteristic roots of system (7) can be computed [17] as the eigenvalues of the following nonlinear eigenvalue problem: where , .

As we know, if all the rightmost eigenvalues of NEP (9), which are the eigenvalues whose real parts are the biggest, have (strict) negative real parts, system (7) will be (asymptotically) stable; otherwise, it will be unstable; see, for example, [18].

Next, we briefly recall the algorithm to compute the rightmost characteristic roots of system (7); see [18] for details.

Assume that is a positive integer; let and the complex moments be By , two Hankel matrices and are defined as follows:

Then, we consider that the above integration is evaluated via a trapezoidal rule on the circle . Let be the number of sample points on the circle and . From the trapezoidal rule, the complex moments can be approximately calculated by the following formula: Based on the Trace-Theorem of Devidenko can be rewritten as

Using , two Hankel matrices and are formed as follows:

Then, the algorithm to compute the rightmost characteristic roots of system (7) is described as follows.

*Algorithm 1 (see [18]). *The following is an algorithm for determining the rightmost eigenvalues of NEP (9) in the disk :(1)Compute the radius of the disk whose center is .(2)Choose the number of sample points on the circle which is the boundary of the disk .(3)Compute .(4)Compute by using (13).(5)Form two Hankel matrices and by using (14) and (15), respectively.(6)Compute all the eigenvalues of the GEP .(7)Compute the smallest singular value of to determine whether is the eigenvalue of NEP (9).(8)Compute the rightmost eigenvalue of NEP (9) in the disk .

*Remark 2. *Algorithm 1 allows for using a larger positive integer .

*Remark 3. *In actual computation, is considered as the eigenvalue of NEP (9) if the smallest singular value of is less than .

Considering what was mentioned above and the stability theory, we can have the following stability criterion for system (1). That is, if all the rightmost eigenvalues of NEP (9) have strict negative real parts, then the positive equilibrium of system (1) will be asymptotically stable; otherwise, it will be unstable.

Therefore, we have the algorithm to determine the local asymptotic stability of the positive equilibrium for system (1) as follows.

*Algorithm 4. *The following is an algorithm for determining the local asymptotic stability of the positive equilibrium for system (1):(1)Compute the unique positive equilibrium of system (1) by solving the algebraic equation (3).(2)Run Algorithm 1 to find the rightmost eigenvalue of NEP (9).(3)**if** the rightmost eigenvalue of NEP (9) has (strict) negative real parts** then**(4)The positive equilibrium of system (1) is asymptotically stable.(5)**else**(6)The positive equilibrium of system (1) is not stable.(7)**end if**

#### 3. Numerical Examples

In this section, we report some numerical examples to study the effectiveness of Algorithm 4 for a set of systems (1). All computations are performed under MATLAB (version 2012b). denotes the rightmost eigenvalue obtained by Algorithm 1.

*Example 1 (see [10]). *We consider the following system: which has a unique positive equilibrium .

For the case when all are equal, we set . By the numerical simulation shown in Figure 1 in which , the positive equilibrium of system (16) is asymptotically stable. And from Figure 2 with , the positive equilibrium loses its stability.