Numerical Simulation of a Class of Three-Dimensional Kolmogorov Model with Chaotic Dynamic Behavior by Using Barycentric Interpolation Collocation Method

This paper numerically simulates three-dimensional Kolmogorov model with chaotic dynamic behavior by barycentric Lagrange interpolation collocationmethod. Some numerical examples are studied for finding some newchaotic behaviors anddemonstrating 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][6][7][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)  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 threedimensional 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][11][12][13][14][15][16][17][18][19][20] studied the dynamics of three dimensional Kolmogorov model with different type of functional response in theory.They found some chaotic dynamics [13,[18][19][20][21][22][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.
So, we can get a numerical solution of (3) and (4).

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  = 40.The accuracy of iteration control  = 10 −10 and the initial iteration value  0 =  0 =  0 = 0;  1 =  1 =  1 = .Experiment .We consider the following three-species food chain model [18]: where  is real parameters, which satisfy the initial condition (0) = 1, (0) = 2, (0) = 1.Results of Experiment 1 are given in Figures 1-2. Figure 1 is obtained by using the current method with  = 1.Among them, ( 1 ) is the time series plot; ( 1 ) is the phase diagram of ; ( 1 ) is the phase diagram of ; ( 1 ) is the graph projected on (, )-plane; ( 1 ) is the three-dimensional space graph.Figure 2 is obtained by using the current method with  = 4.2.( 2 ) is the time series plot; ( 2 ) is the phase diagram of ; ( 2 ) is the three-dimensional space graph; ( 2 ) is the graph projected on (, )-plane; ( 2 ) 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 .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 (1): (0) = 1, (0) = 1, (0) = 1 and the initial conditions ( 2 ) are the three-dimensional space graph; ( 1 ) and ( 2 ) are the graph projected on (, )-plane; ( 1 ) and ( 2 ) are the graph projected on (, )-plane; ( 1 ) and ( 2 ) 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 .
Figure 5 is obtained by using the current method.() is time series plot and () is the three-dimensional space graph with  = 7; () is time series plot and () is the threedimensional space graph with  = 9; () is time series plot and () is the three-dimensional space graph with  = 12.We can see that, from the figure, the fluctuation range of , ,  at  = 900 is very large when  = 7; the fluctuation range of , ,  at  = 900 is obviously reduced when  = 9; , ,  at  = 900 almost all tend to be stationary when  = 12.
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.( 1 ) is time series plot and ( 1 ) is the three-dimensional space graph with  = 1.4; ( 2 ) is time series plot and ( 2 ) is the threedimensional space graph with  = 1.575; ( 3 ) is time series plot and ( 3 ) is the three-dimensional space graph with  = 1.62.Experiment .We consider the following model with a Watttype 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.

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.

Table 1 :
Parameters used in Experiment 6.

Table 2 :
Parameters used in Experiment 6.