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.