#### Abstract

This paper numerically simulates three-dimensional Kolmogorov model with chaotic dynamic behavior by barycentric Lagrange interpolation collocation method. Some numerical examples are studied for finding some new chaotic behaviors and demonstrating some existing chaotic dynamic behaviors of the Kolmogorov model. Results obtained by the present method indicate that the method has merits of small operations and good numerical stability.

#### 1. Introduction

J.P. Berrut [1, 2] introduced barycentric Lagrange interpolation and studied its numerical stability and convergence. Barycentric Lagrange interpolation is unconditionally stable at the Chebyshev points. S.P. Li, Z.Q. Wang [3, 4] gave some algorithms of barycentric Lagrange interpolation collocation method (BLICM). Some authors [3, 5–8] solved all sorts of equations and showed the BLICM has merits of small operations and high precision (see [3, 4, 9]). This paper numerically simulates some three-dimensional Kolmogorov models with chaotic dynamic behavior. The purpose of this paper is to find some new chaotic behaviors and verify the existing chaotic dynamic behaviors by the BLICM.

Three-dimensional Kolmogorov models comprise a significant class of ecological models that are used widely in ecology to represent the dynamic behavior of prey and predators, which are expressed in the following form:where represents the population density of the species and represents the per capita growth rate of the species.

In ecology, the most frequently used model is the Lotka-Volterra system; that is, each per capita growth function is affine and chosen as the logistic growth. In this circumstance, model (1) reads as

Model (2) is a totally competitive system if all parameters are positive.

There are several famous functional responses in the Kolmogorov model which are referred to as Holling type I, type II, type III, type IV, Monod-Haldane type, Hassel-Verley type, Beddington-DeAngelis type functional response, etc.

In this paper, we consider the following three-dimensional Kolmogorov model with functional response:with the following initial condition:where represent the population density of the species and is the known functional response. are known constants. are unknown functions of time .

Many researchers [10–20] studied the dynamics of three dimensional Kolmogorov model with different type of functional response in theory. They found some chaotic dynamics [13, 18–23] of the three-dimensional Kolmogorov model. Chaos and hyperchaos exist in many natural processes and are one of the main contents of nonlinear science research. Although many kinds of numerical methods of the Kolmogorov model have been announced, simple and efficient methods have always been the direction that scholars strive to pursue. This paper suggests the BLICM to solve the three-dimensional Kolmogorov model. Model (3) is adopted as an example to elucidate the solution process.

#### 2. Barycentric Lagrange Interpolation Collocation Method

First of all, we discrete computational interval by Chebyshev points into and construct following linear iterative format of model (3):

The format (5) is convergent; then , , .

Next, we transform format (5) into the following linear algebraic equations.where is order matrix. is order unit matrix and is a symbol of diagonal matrix composed of vectors. is, respectively, barycentric interpolation primary function and is center of gravity Lagrange interpolation weight. The vector

The vector

The last and the first line of equations (6) are replaced separately by the equation of the initial condition (9) in turn.

So, we can get a numerical solution of (3) and (4).

#### 3. Numerical Experiments

In this section, some numerical examples are studied to demonstrate the accuracy of the present method. The examples are computed using MatlabR2017a. In numerical experiments, the number of nodes . The accuracy of iteration control and the initial iteration value .

*Experiment 1. *We consider the following three-species food chain model [18]:where is real parameters, which satisfy the initial condition Results of Experiment 1 are given in Figures 1-2.

Figure 1 is obtained by using the current method with . Among them, is the time series plot; is the phase diagram of ; is the phase diagram of ; is the graph projected on -plane; is the three-dimensional space graph. Figure 2 is obtained by using the current method with . is the time series plot; is the phase diagram of ; is the three-dimensional space graph; is the graph projected on -plane; is the the graph projected on -plane; is the graph projected on -plane.

Figure 1 gets some new chaotic behaviors. Figure 2 verifies the existing chaotic dynamic behaviors [18]. Our study suggests that model (10) will go chaotic when the rate of the self-reproduction of the prey is large. Our numerical results are in good agreement with the theory [18].

*Experiment 2. *We consider the following tritrophic food chain model [19]:In comparison with [19], we get better numerical results; results of Experiment 2 are given in Figures 3-4. Our numerical results are in good agreement with the theory [19].

Figures 3 and 4 are obtained by using the current method with the initial conditions : and the initial conditions : , respectively. and are the time series plot of ; and are the phase diagram of ; and are the three-dimensional space graph; and are the graph projected on -plane; and are the graph projected on -plane; and are the graph projected on -plane. It can be seen that the change of initial conditions leads to significant changes in the time series diagrams of , and .

*Experiment 3. *We consider the following tritrophic food chain model [20]:where is real parameters, which satisfy the initial condition Our numerical simulations and experimental results of [20] are confirmed. Results of Experiment 3 are given in Figure 5.

Figure 5 is obtained by using the current method. is time series plot and is the three-dimensional space graph with ; is time series plot and is the three-dimensional space graph with ; is time series plot and is the three-dimensional space graph with . We can see that, from the figure, the fluctuation range of at is very large when ; the fluctuation range of at is obviously reduced when ; at almost all tend to be stationary when .

*Experiment 4. *We consider the following model with Beddington-DeAngelis-type functional response [21]:where , and are real parameters, which satisfy the initial condition We choose the parameters . Results of Experiment 4 are given in Figures 6-7.

Figure 6 is time series plots and the three-dimensional space graphs of a predator-prey model with Beddington-DeAngelis type functional response for Experiment 4 with different parameters obtained by using the current method. is time series plot and is the three-dimensional space graph with ; is time series plot and is the three-dimensional space graph with ; is time series plot and is the three-dimensional space graph with . Figure 7 is obtained by using the current method with and , respectively. and are the graph projected on -plane; and are the graph projected on -plane; and are the graph projected on -plane.

*Experiment 5. *We consider the following model with Monod-Haldane type response function [22]:where , and are real parameters, which satisfy the initial condition We choose the parameter .

Our numerical simulation results are in good agreement with the theory [22]. We give some new chaotic behaviors. Results of Experiment 5 are given in Figures 8-9.

Figure 8 is time series plots and the three-dimensional space graphs of a predator-prey system with Monod-Haldane type response function for Experiment 5 with different parameters obtained by using the current method. is time series plot and is the three-dimensional space graph with ; is time series plot and is the three-dimensional space graph with ; is time series plot and is the three-dimensional space graph with . Figure 9 is obtained by using the current method with and , respectively. and are the graph projected on -plane; and are the graph projected on -plane; and are the graph projected on -plane.

*Experiment 6. *We consider the following model with a Watt-type functional response [23]:Our numerical simulation results are in good agreement with the theory [23]. Compared with [23], we get better numerical results and give some new chaotic behaviors. Results of Experiment 6 are given in Figures 10–15. We choose the parameter; see Tables 1-2.

Figures 10–12 are obtained by using the current method. , , and are the time series plot of ; , , and are the phase diagram of ; , , and are the phase diagram of ; , , and are the graph projected on -plane; , , and are the three-dimensional space graph. Figures 13–15 are obtained by using the current method. , , and are the time series plot of ; , , and are the phase diagram of ; , , and are the phase diagram of ; , , and are the graph projected on -plane; , , and are the three-dimensional space graph.

#### 4. Conclusions and Remarks

In this paper, the three dimensional Kolmogorov model was solved by using the barycentric Lagrange interpolation collocation method. These numerical experiments illustrate that numerical results of the present method are the same as experimental results.

All computations are performed by the MatlabR2007b software packages.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

#### Authors’ Contributions

Mingjing Du carried out the main part of this manuscript. Wei Zhang obtained data and analyze numerical simulation results, which greatly improved the quality of the paper. Junmei Li carried out numerical simulation and gave some numerical results. Yulan Wang made the program of the barycentric interpolation collocation method in the revised version and corrected the derived results. All authors read and approved the final manuscript. When we resubmitted the revised version, we forgot to change the author list.

#### Acknowledgments

The authors thank the reviewers for their valuable suggestions, which greatly improved the quality of the paper. This paper is supported by the Natural Science Foundation of Inner Mongolia , Inner Mongolia Maker Collaborative Innovation Center of Jining Normal University, and the National Natural Science Foundation of China .