This paper centres on the application of the new piecewise successive linearization method (PSLM) in solving the chaotic and nonchaotic Chen system. Numerical simulations are presented graphically and comparison is made between the PSLM and Runge-Kutta-based methods. The work shows that the proposed method provides good accuracy and can be easily extended to other dynamical systems including those that are chaotic in nature.

1. Introduction

In this paper, we consider the Chen dynamical system [1], which is a three-dimensional system of ordinary differential equations with quadratic nonlinearities, defined as where , , are dependent variables and , , are positive constants. The Chen system (1.1) does not admit a closed form solution and it is known to exhibit both chaotic and nonchaotic behaviors for distinct parameter values. Bifurcation studies [24] have shown that with the parameters and , system (1.1) exhibits nonchaotic behavior and chaotic behavior when and , respectively. Finding accurate and efficient methods for solving chaotic systems such as the Chen system has been an active research area for many researchers who use such systems as a benchmark for testing the accuracy and performance of methods of solutions under their development. Other variants of the Chen system such as fractional order and four-dimensional Chen systems have been the subject of study by many scholars.

The solution of nonchaotic, chaotic and hyperchaotic Chen system equations has been obtained by different analytical and numerical methods such as Runge-Kutta-based numerical schemes. In recent years, several analytical and seminumerical methods have been utilized in solving the Chen system equations [5, 6]. Because the solutions for the Chen systems rapidly oscillate within small intervals, conventional analytical and numerical algorithms have been implemented in a sequence of small multiple intervals in order to ensure accuracy convergence of the solution techniques. Recent multistage methods that have been reported for solving Chen systems include the multistage differential transform method [7, 8], the multistage homotopy analysis method [9], the multistep homotopy perturbation method [10] and the multistage variational iteration method [11].

The aim of this paper is to obtain the solution of the Chen system using the successive linearization method (SLM) that is implemented in a sequence of intervals. This new algorithm is called the piecewise successive linearization method (PSLM). The standard SLM was developed primarily for solving nonlinear boundary value problems (see, e.g., [1217]). To the best of our knowledge, the method has not been extended to initial value problem (IVPs) and in particular chaotic IVPs of the type (1.1). We remark that the standard SLM is not suitable for solving dynamical systems of the type (1.1), hence the need to modify the method and implement it on multiple intervals. The paper focuses on the use of the proposed PSLM in solving the Chen system (1.1) using selected parameters that yield both chaotic and nonchaotic results. Numerical results are presented graphically and are compared with Runge-Kutta-based results and excellent agreement is observed.

The remainder of the paper is organized as follows. In Section 2, the basic description of the standard SLM is briefly presented. In Section 3, the PSLM is introduced. The application of the proposed method is numerically investigated in Section 4 to generate results which are discussed therein. Finally, conclusions in Section 5 close the paper.

2. Successive Linearization Method (SLM)

In this section, we give a brief description of the basic idea behind the successive linearization method (SLM) [1315]. This will be followed by a description of the multistage extension of the SLM algorithm which is applied to the Chen system (1.1).

The SLM approach assumes that the solution of (1.1) can be expressed as Substituting (2.1) in the governing equation (1.1) and neglecting nonlinear terms in we obtain subject to the initial conditions where Starting from a suitable initial approximation , , , the solutions for , , (for ) can be obtained by iteratively solving the linear system (2.2)–(2.4). A suitable initial guess to start off the SLM algorithm (2.2)–(2.4) is obtained by solving the linear part of (1.1) subject to the problem's initial condition where are given initial values. For this system of first-order equations the initial approximation can easily be solved analytically in any given interval . However, the solution of the SLM algorithm (2.2)–(2.4) cannot be solved exactly. Spectral collocation methods (or any other numerical method) can be used to solve (2.2)–(2.4) in any given interval . After obtaining solutions for (2.3), the approximate solution , , and is determined as the series solution Similar expressions are obtained for and . A SLM solution is said to be of order if the above series is truncated at , that is, if

3. Piecewise Successive Linearization Method

It was observed through numerical experimentation that, in its standard form, the SLM will not give accurate solutions for initial value problems of the Chen system type. Thus, it can be remarked that the SLM is ideally suited for boundary value problems whose solutions do not rapidly change in behaviour or oscillate over small regions of the domain of the governing problem. The SLM solution can thus be considered to be local in nature and may not be suitable for initial value problems at very large values of the independent variable . A simple way of ensuring the validity of the approximations for all values of is to determine the solution in a sequence of equal intervals, which are subject to continuity conditions at the end points of each interval. To extend this solution over the interval , we divide the interval into sub-intervals , where . We solve (2.2)–(2.4) in each subinterval . When applied to the multiple subintervals the SLM will be referred to as the piecewise SLM (PSLM).

Let , , be the solution of (2.2)–(2.4) in the first subinterval and , , be the solutions in the subintervals for . The initial conditions used in obtaining the solutions in the subinterval are obtained from the initial conditions of the subinterval . Thus, we solve subject to the initial conditions The initial approximations for solving (3.1) are obtained as solutions of the system subject to the initial condition An exact solution of (3.3) can easily be found. To solve (3.1) we use the Chebyshev spectral collocation method on each interval . Before applying the spectral method the region is transformed to the interval on which the spectral collocation method is defined. This can be achieved by using the linear transformation in each interval (for ), After the transformation, the interval is discretized using the Chebyshev-Gauss-Lobatto collocation points [18, 19] which are the extrema of the th-order Chebyshev polynomial The basic idea behind the pseudo-spectral collocation method is the introduction of a differentiation matrix which is used to approximate the derivatives of the unknown variables , , at the collocation points as the matrix vector product where and is the vector function at the collocation points . Similar expressions can be obtained for and . The Chebyshev derivative matrix is of size and its entries are defined [18, 19] as with Applying the Chebyshev spectral collocation method in (3.1) gives where and is an identity matrix, is an matrix of zeros. The denotes an diagonal matrix corresponding to that is evaluated at the collocation points and is an vector corresponding to the quantities that have been evaluated at the collocation points. Thus, starting from the initial approximation obtained as a solution of (3.3)-(3.4), the recurrence formula can be used to obtain the solutions , , in the interval . The solution approximating , for example, in the entire interval is given by Similar expressions can be obtained for and . We observe that when , the proposed piecewise successive linearization method (PSLM) becomes equivalent to the original SLM algorithm.

4. Results and Discussion

In this section, numerical simulations are given to validate the proposed piecewise successive linearization method. The PSLM results were obtained using collocation points and ten iterations (that is ) in each interval. Unless otherwise specified, the width of each interval was taken to be . We also fix the values of the parameters and with for the nonchaotic case and for the chaotic case. The initial conditions are , , and . The results for all the simulations done in this paper are computed for the time span . In order to assess the accuracy and performance of the proposed PSLM approach, the present results are compared with those obtained with the MATLAB built-in solver . The solver integrates a system of ordinary differential equations using explicit 4th and 5th Runge-Kutta (4,5) formula, the Dormand-Prince pair [20].

The results of the PSLM simulation of the Chen system for the nonchaotic case are shown in Figures 1 and 2. In Figure 1 the results of the comparison between the PSLM and Runge-Kutta results for the time series evolution of the state variables are given. It can be seen that there is good agreement between the two results. Figure 2 depicts the phase portraits for the nonchaotic case. We observe that the system converges to a periodic orbit and results in a limit cycle as shown in Figure 2 when the parameters , , and are used.

The results of the PSLM simulation of the Chen system for the chaotic case are shown in Figures 3 and 4. In Figure 3, the PSLM results are compared with the Runge-Kutta results when parameters , , and are used. For this set of parameters, the system becomes chaotic and the chaotic attractors are shown in Figure 4. We observe that there is good agreement between the PSLM and Runge-Kutta results in the time series illustration depicted in Figure 3. This shows that the proposed method of solution has great potential of being a useful tool for solving complex dynamical systems with chaotic and nonchaotic behaviour.

5. Conclusion

In this paper, we carefully develop a novel piecewise successive linearization method (PSLM) which seeks to overcome the difficulties arising when the successive linearization method is extended to initial value problems with complex behaviour. The proposed method of solution is implemented on the Chen system and results are generated using known parameters that give both chaotic and nonchaotic outcomes. The validity of PSLM is verified by comparing its results with the results of built in Matlab-based numerical Runge-Kutta method, . Close agreement of the two sets of results is observed thus demonstrating the accuracy of the PSLM approach for the particular problem considered. In conclusion, the PSLM may be considered as a robust refinement and extension of the existing SLM. The extended method may also be valid for other complex nonlinear initial value problems with chaotic behaviour. This paper can be used as a standard example for other applications in engineering and applied sciences.