Abstract

The paper presents a new method, called the Polynomial Least Squares Method (PLSM). PLSM allows us to compute approximate analytical solutions for the Brusselator system, which is a fractional-order system of nonlinear differential equations.

1. Introduction

In recent years, in many practical applications in various fields such as physics, mechanics, chemistry, and biology (see, e.g., [16]), the problems being studied are modeled using fractional nonlinear equations.

For most of such fractional nonlinear equations, the exact solutions cannot be found and, as a consequence, a numerical solution or, if possible, an analytical approximate solution of these equations is sought.

Due to the complexity of this type of problems, a general approximation algorithm does not exist and, thus, various approximation methods, each with its strong and weak points, were proposed, including, among others:(i)the Adomian Decomposition Method [79],(ii)the Homotopy Analysis Method [1012],(iii)the Homotopy Perturbation Method [13, 14],(iv)the Laplace Transform Method [15, 16],(v)the Fourier Transform Method [17],(vi)the Variational Iteration Method [1820].

The objective of our paper is to present the Polynomial Least Squares Method (PLSM), which allows us to compute approximate analytical solutions for the Brusselator system. The fractional order Brusselator system was recently studied by several authors [2123] and can be expressed as follows.

We consider the following Brusselator system:together with the initial conditions:where , , , , are real constants, and denotes Caputo’s fractional derivative [15]:

In the next section we will introduce PLSM for the Brusselator system and in the third section we will compare the approximate solutions obtained by using PLSM with the approximate solutions from [20]. The computations show that the approximations computed by using our method present an error smaller than the error of the corresponding solutions from [20].

2. The Polynomial Least Squares Method

For problem ((1), (2)) we consider the remainder operators:

We will find approximate polynomial solutions of problem ((1), (2)) on the interval, solutions which satisfy the following conditions:

Definition 1. One calls an -approximate polynomial solution of problem ((1), (2)) an approximate polynomial solution satisfying relations ((5), (6)).

Definition 2. One calls a weak  -approximate polynomial solution of problem ((1), (2)) an approximate polynomial solution satisfying the relations: together with initial conditions (6).

Definition 3. One considers the sequence of polynomials: , , , , satisfying the conditions:One calls the sequence of polynomials convergent to the solution of problem ((1), (2)) if .

We will find weak -polynomial solutions of the type:where the constants , are calculated using the steps outlined as follows.(i)We attach to problem ((1), (2)) the following real functional:where are computed as functions of by using initial conditions (6).(ii)We computes the values of as the values which give the minimum of functional (10) and the values of again as functions of by using the initial conditions.(iii)Using the constants thus determined, we consider the polynomials:

The following convergence theorem holds.

Theorem 4. The necessary condition for problem ((1), (2)) to admit sequences of polynomials convergent to the solution of this problem is Moreover, , , such that, , , it follows that , , are weak -approximate polynomial solutions of problem ((1), (2)).

Proof. Based on the way the coefficients of the polynomials are computed and taking into account relations ((9)–(11)), the following inequality holds: It follows that We obtain From this limit we obtain that, such that, , , it follows that are weak -approximate polynomial solutions of problem ((1), (2)), .

Remark 5. Any -approximate polynomial solutions of problem ((1), (2)) are also weak -approximate polynomial solutions, but the opposite is not always true. It follows that the set of weak approximate solutions of problem ((1), (2)) also contains the approximate solutions of the system.

Taking into account the above remark, in order to find -approximate polynomial solutions of problem ((1), (2)) by PLSM we will first determine weak approximate polynomial solutions, . If , , then are also -approximate polynomial solutions of the system.

3. Application: the Fractional-Order Brusselator System

We consider the following fractional-order Brusselator system [20]:together with the initial conditions:

In [20] approximate solutions of (17) are computed using the Variational Iteration Method (VIM) for the case . Also, a comparison with numerical solutions is presented for the particular case , illustrating the applicability of the method.

3.1. The Case

For the case , using PLSM with , we obtain the following approximate polynomial solutions:

In Figures 1 and 2 we compare these approximations with the corresponding approximations of the same order computed by VIM (relations in [20]), obtaining a good agreement. In Figures 3 and 4 we compare the expressions of remainders (4) obtained by replacing the approximate solutions back in the equations. It is easy to observe that the errors obtained by using PLSM are smaller than the ones obtained by using VIM.

3.2. The Case

For the case , using PLSM with we obtain the following approximate polynomial solutions:

In this case both approximations (VIM and PLSM) consist of third-order polynomials.

We omitted the figures which compare our approximations with the ones given by VIM since they look almost the same as the corresponding ones from the case .

However, in the case it is possible to compute the absolute error corresponding to an approximate solution as the difference in absolute value between the approximate solution and the numerical solution (in this case computed by using the Wolfram Mathematica software).

Figures 5 and 6 present the comparison between the absolute errors corresponding to the approximate solutions from [20] obtained by VIM and the absolute errors corresponding to our approximate solutions.

Again, it is easy to observe that the errors obtained by using PLSM are smaller than the ones obtained by using VIM.

4. Conclusions

In this paper we present the Polynomial Least Squares Method, which is a relatively straightforward and efficient method to compute approximate solutions for the fractional-order Brusselator system.

The comparison with previous results illustrates the accuracy of the method, since we were able to compute more precise approximations than the previously computed ones.

In closing we mention the fact that, due to the nature of the method, it is relatively easy to extend PLSM for the general case of fractional systems of nonlinear differential equations.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.