So far, a lot of works published on fractional order linear/non-linear differential equations but there are still works have to be done. In this work, we aim to extend the Hermite Collocation method (HCM) for obtaining solution to a system of fractional order differential equations with variable coefficients and specified initial conditions. The technique constructs an analytical solution of the form of a truncated Hermite series with unknown coefficients. The orthogonal Hermite polynomials have the importance in the theory of light fluctuations and quantum states and, in particular, some problems of coastal hydrodynamics and meteorology  . This method is the adaptation of Taylor collocation method with Hermite polynomials and first has been used to solve higher-order linear Fredholm integro differential equation in  and the development of the method can be found in  .
This paper is organized as follows. Section 2 involves some basic definitions and properties of fractional calculus. In Section 3, the theory and definitions of Hermite collocation method and the construction of this method for fractional order systems are presented. In Section 4, the matrix relations for initial conditions are defined and the Section 5 deals with the error estimate for the method. Section 6 involves some illustrative examples. Finally, the last section concludes with some remarks based on the reported research.
2. Preliminary and Notations
We first recall the following known definitions and preliminary facts of fractional derivatives and integrals which are used throughout this paper.
Definition 2.1. (   ) Let. The Riemann Liouville fractional integral of a function f of order is defined by
where is the gamma function and. For consistency, we take, which is identity operator and holds.
Definition 2.2. (   ) The Riemann Liouville fractional derivative of order of a function f is defined by
where and is defined as the integer part of. Again, for consistency, , then, it follows where.
Alternatively, we recall the following definition of Caputo (  ) for fractional derivatives and Caputo’s definition is much more suitable for identifiable physical states, i.e. initial or boundary conditions. Therefore, all derivatives will be in Caputo sense throughout the paper.
Definition 2.3. The Caputo fractional derivative of order of a function f on an interval is defined by
Some properties of Caputo derivative are given as follows:
1) (  ) Let and or then
2) (  ) Let and. If or, then
3) (  ) For every
4) Let and for some then,
3. Establishing Hermite-Collocation Method for Fractional Order Systems with Variable Coefficients
In this section, we will consider the following system of fractional order differential equations (FDEs) with variable coefficients, ,
where and, are continuous functions on. The initial conditions are defined as
In Equation (5), are some given constants and we denote for simplicity. Here, we assume that the approximate solution of the problem is given in terms of truncated Hermite series,
where defines the unknown Hermite coefficients of the solution and N is a positive integer which is chosen sufficiently small for avoiding the laborious work such that. Therefore, the fundamental matrix relation of Equation (4) can be written as
where and are defined as follows:
Now, we need to define the Caputo derivatives of . By using Equation (6), therefore, we write (see  ),
where and are defined as and respectively. Now, we will describe the matrix representation of the truncated Hermite series in terms of rational power of the indepandant variable x, by using the following generalized formula:
where and Now, in terms of N being odd or even, we denote the truncated series in matrix notation such as follows (see  ):
If N is an odd number:
if N is even then,
Hence, we have or equivalently,
and then, substitution of Equation (11) into Equation (6) yields,
Now, the nath order Caputo derivative of Equation (12) is written as
Here, the matrix B is defined as follows ():
Hence, if we substitute Equation (14) into Equation (13) we have:
Therefore, the matrices in Equation (15), for, are clearly shown by
where the each submatrix, consisting of rows and columns. Consequently, the above matrix equation can be written as,
where appears as consisting of k rows and columns. Hence, inserting the collocation points, , , into Equation (7) then, we have
where and G are of the form:
Apart from this, arranging Equation (16) for each collocation points then, we can write explicitly as,
Therefore, the matrix form is equivalent to
and each submatrix in is denoted by,
Consequently, now we denote Equation (7) of the form:
Then, by writing Equation (18) in Equation (19), the matrix form of the system of FDEs is written by
hence, the system of FDEs is simply shown by
Now, Equation (21) constructs an algebraic system. To obtain the solution of the above system, the augmented matrix is written as follows:
Solving the above system, as a result, we obtain the desired Hermite coefficients in the truncated Hermite series. Hence, writing in Equation (12) we evaluate the unknowns of the system of FDEs (Equation (4)).
4. Matrix Relations for Initial Conditions
In generally, we look for the solution of the system of FDEs under specified conditions. However, preceding calculations do not involve these conditions. Therefore, we need to incorporate these conditions into the work. Then, we have to establish the new form of Equation (22) which involves initial conditions, Equation (5). Now, we start by writing Equation (5) explicitly for each same as below:
Hence, by using the above relations, we obtain t-conditions for each unknown,. For example, for we obtain t conditions such as follows:
Therefore, the conditions in matrix notation fulfils,
and for, we define
Now writing Equation (16) into Equation (23) for, we obtain
Now, calling U as,
then, the Hermite polynomial coefficients matrix which corresponds to the given initial conditions (Equation (5)), can be written as
In Equation (25), U involves kt rows and columns. Consequently, deleting rows in Equation (21) and then replacing these rows by Equation (25), we obtain the whole augmented matrix of the system, , as follows:
Hence, the system of algebraic equations of which unknowns are the hermite polynomial coefficients are shown by
Theorem 1. If, (i.e.) then,
By the above theorem, the matrix of Hermite coefficients, A is uniquely determined by Equation (27). Finally, substitution of these coefficients into the truncated Hermite series gives the desired solution of the form:
5. Error Estimate for the Solution
The truncated Hermite series, Equation (29), is the approximate solution of Equation (4) with the given initial conditions, (Equation (5)). Since this solution should approximately satisfy the Equation (4) hence, the residuals
give the error at the particular points ,. Let us now denote the residuals by as an error function. The error should be approximately zero or where is any positive constant. If the is prescribed before then, the truncation limit for N is increased until becomes smaller than (see   ).
6. Numerical Applications
The technique which we have developed to solve fractional order systems is quite feasible and accurate. To show the accuracy of the method the following system of FDEs with variable coefficients are solved. All the numerical calculations have been performed by using MatlabR2007b.
Example 6.1. We first consider the problem, which is mentioned in  (),
with given initial conditions,
Now, we will look for a solution to the system of FDEs in terms of Hermite polynomials of the form;
Here, we will take into consideration:. Since and, then it requires that. As it should be, therefore, we can select for convenience. Now, Equation (4) can be rewritten as:
By using Equations ((6) and (20)) then, the matrix form of the system, Equation (30), is written by
where, , , , , , , and the rest are zero. From here, we evaluate that. Hence, the matrix form of the Equation (30) is deduced as,
For, the collocation points are and. Then the matrices in Equation (34) become,
Hence, we have
Therefore, evaluating Equation (34), we obtain W as,
then, the augmented matrix for the system, or, is obtained as,
where the matrix consisting of k rows and columns and similarly, the matrix form of the initial conditions, Equation (31), is obtained from Equation (24) such as follows:
Hence, we have,
Then, by substituting the related matrices into Equation (37), the augmented matrix is obtained as
Moreover, deleting the last two rows of Equation (36) and replacing the matrix in Equation (38), we have
Since then, the solution of the resulting linear system, gives coefficient matrix, which is equivalent to
In conclusion, writing these coefficients into:
we obtain the solutions of the system of FDEs as follows
Example 6.2. In  , the authors have modeled the pollutant problem in a lake which connected by channels by the following fractional order system (see Figure 3),
where they considered:, and the initial conditions were de- fined as In  , the authors solved the following
Figure 1. Approximate solution of the system in Example 1 by HCM.
Figure 2. (a)Approximate solution of the system in Example 1 by Differential Transform Method, (b) Adomian Decomposition Method.
Figure 3. Pollutant problem scheme of three lakes which connected by the channels.
ordinary system by Bessel Polynomial Collocation method (BCM) with the assumptions: and the initial conditions;
Now, we will solve the fractional form of Equation (44), which is defined as in Equation (44) by HCM method. We consider here the case: and. Since then, it requires that. As a result, we can choose. Therefore, the solution will be of the form:
The fundamental matrix form of the system of Equation (44) is obtained from Equations ((6) and (20)) such as follows,
which is equivalent to:
Then, by performing the calculations, we obtain the following matrices:
then, the agumented matrix is obtained at collocation points as follows:
Since, then, the coefficient matrix, is obtained. When these coefficients substituted into
Figure 4. Solution of the system in Example 2 by HCM and comparison with BCM method for, and (a) Solution of, (b) Solution of, (c) Solution of.
the solution of the system, Equation (44), is obtained as follows;
Figure 4 shows the plots of and, which are the solutions of Example 6.2 respectively. In these plots, the results have been compared by BCM method and our method (HCM) for. Furthermore, each plots also shows the results for, which exists first time in the literature and there is a clear difference between the solution of the fractional order sytem and ordinary differential equation system although, there is a small change between.
The basic goal of this work is to employ HCM method to obtain solution for a system of fractional order differential equations. These types of systems with variable coefficients are usually difficult to solve analytically. However, the presented method provides considerable simplifications in the solution. The coefficients of truncated Hermite series can be evaluated easily by the help of any symbolic computer packages. The obtained results demonstrate the reliability of the algorithm and give us a wider applicability to fractional higher order systems.
 Odibat, Z.M. (2009) Computational Algorithms for Computing the Fractional Derivatives of Functions. Mathematics and Computers in Simulation, 79, 2013-2020.
 Ford, N.J. and Joseph Connolly, A. (2009) Systems-Based Decomposition Schemes for the Approximate Solution of Multi-Term Fractional Differential Equations. Journal of Computational and Applied Mathematics, 229, 382-391.
 Ford, N.J. and Charles Simpson, A. (2001) The Numerical Solution of Fractional Differential Equations: Speed versus Accuracy. Numerical Algorithms, 26, 333-346.
 Sweilam, N.H., Khader, M.M. and Al-Bar, R.F. (2007) Numerical Studies for a Multi Order Fractional Differential Equations. Physics Letters A, 371, 26-33.
 Momani, S. and Odibat, Z. (2006) Analytical Solution of a Time-Fractional Navier-Stokes Equation by Adomian Decomposition Method. Applied Mathematics and Computation, 177, 488-494.
 Odibat, Z. and Momani, S. (2007) Numerical Approach to Differential Equations of Fractional Order. Journal of Computational and Applied Mathematics, 207, 96-110.
 Odibat, Z. and Momani, S. (2008) Numerical Methods for Nonlinear Partial Differential Equations of Fractional Order. Applied Mathematical Modelling, 32, 28-39.
 Odibat, Z. and Momani, S. (2006) Application of Variational Iteration Method to Equation of Fractional Order. International Journal of Nonlinear Sciences and Numerical Simulation, 7, 271-279.
 Wu, J.L. (2009) A Wavelet Operational Method for Solving Fractional Partial Differential Equations Numerically. Applied Mathematics and Computation, 214, 31-40.
 Ertürk, V.S. and Momani, S. (2008) Solving Systems of Fractional Differential Equations Using Differential Transform Method. Journal of Computational and Applied Mathematics, 215, 142-151.
 Yalcnbas, S., Konuralp, A., Demir, D.D. and Sorkun, H.H. (2010) The Solution of the Fractional Differential Equation with the Generalized Taylor Collocation Method. International Journal of Research & Reviews in Applied Sciences, 4, 296-303.
 Kumar, P. and Agrawal, O.P. (2006) Numerical Scheme for the Solution of Fractional Differential Equations of Order Greater than One. Journal of Computational and Nonlinear Dynamics, 1, 178-185.
 Dattoli, G. (2004) Laguerra and Generalized Hermite Polynomials: The Point of View of the Operation Method. Integral Transforms and Special Functions, 15, 93-99.
 Akgonüllü, N., Sahin, N. and Sezer, M. (2010) A Hermite Collocation Method for the Approximate Solutions of Higher-Order Linear Fredholm Integro-Differential Equations, Numerical Methods for Partial Differential Equations, 27, 1708-1721.
 Caputo, M. (1967) Linear Models of Dissipation Whose Q Is Almost Frequency Independent II. Geophysical Journal of the Royal Astronomical Society, 13, 529-539.
 Biazar, J., Shahbala, M. and Ebrahimi, H. (2010) VIM for Solving the Pollution Problem of a System of Lakes. Journal of Control Science and Engineering, 2010, Article ID: 829152.
 Yüzbas, S., Sahin, N. and Sezer, M. (2012) A Collocation Approach to Solving the Model of Pollution for a System of Lakes. Mathematical and Computer Modelling, 55, 330-341.