A numerical scheme for solutions of the Chen system

In this paper, we will develop the Bessel collocation method to find approximate solutions of the Chen system, which is a three‐dimensional system of ODEs with quadratic nonlinearities. This scheme consists of reducing the problem to a nonlinear algebraic equation system by expanding the approximate solutions by means of the Bessel polynomials with unknown coefficients. By help of the collocation points and the matrix operations of derivatives, the unknown coefficients of the Bessel polynomials are calculated. The accuracy and efficiency of the proposed approach are demonstrated by two numerical examples and performed with the aid of a computer code written in MAPLE. In addition, comparisons between our method and the homotopy perturbation method numerical solutions are made with the accuracy of solutions. Copyright © 2012 John Wiley & Sons, Ltd.