#### Abstract

We present a new version of the reproducing kernel Hilbert space method (RKHSM) for the solution of systems of fractional integrodifferential equations. In this approach, the solution is obtained as a convergent series with easily computable components. Several illustrative examples are given to demonstrate the effectiveness of the present method. The method described in this paper is expected to be further employed to solve similar nonlinear problems in fractional calculus.

#### 1. Introduction

In this paper, we consider the following system of fractional integrodifferential equations: where , , , and is derivative of order in the sense of Caputo and , subject to the initial conditions:

In the last two decades, fractional calculus has found diverse applications in various scientific and technological fields [1, 2], such as thermal engineering, acoustics, electromagnetism, control, robotics, viscoelasticity, diffusion, edge detection, turbulence, signal processing, and many other physical and biological processes. Fractional differential equations have also been applied in modeling many physical and engineering problems. Most systems of fractional integrodifferential equations do not have exact solutions, so numerical techniques are used to solve such systems. The homotopy perturbation method, the Adomian decomposition method, and other methods are used to give an approximate solution to linear and nonlinear problems; see [3–13] and the references therein.

In our previous work [14], we proposed a reproducing kernel Hilbert space method for solving integrodifferential equations of fractional order based on the reproducing kernel theory [14, 15]. In this paper, we will generalize the idea of the RKHSM to provide a numerical solution for systems of fractional integrodifferential equations (1). To demonstrate the effectiveness of the RKHSM algorithm, several numerical experiments of linear and nonlinear systems of fractional equations (1) will be presented.

This paper is organized as follows. An introduction of the algorithm for solving systems of fractional integrodifferential equations is given in Section 2. In Section 3, we introduce several examples to show the efficiency of the method. Finally, a conclusion is given in Section 4.

#### 2. The Algorithm

After homogenizing the initial conditions (2), we apply the operator , the Riemann-Liouville fractional integral of order [2, 16–20], to both sides of (1) to have where

It is clear that (3) is equivalent to (1), so every solution of the integral equation (3) is also a solution of our original problem (1) and vice versa.

To solve (3) by means of the reproducing kernel Hilbert space method, first, we need to construct a reproducing kernel of certain spaces is absolutely continuous, , and in which every function satisfies the homogenous initial conditions of (1).(i)The inner product of the space is absolutely continuous real value function, is given by and norm . In [21], Li and Cui proved that is a reproducing kernel Hilbert space and its reproducing kernel is given by (ii)The inner product of the space are absolutely continuous real value functions, is given by and norm . is a reproducing kernel Hilbert space and its reproducing kernel is given by (iii)The inner product of the space are absolutely continuous real value functions, is given by and norm . is a reproducing kernel Hilbert space and its reproducing kernel is given by where .

The method of obtaining the reproducing kernel can be found in [15].

Let such that . Then , , are bounded linear operators.

Let be a countable dense set in . Let and , where is the adjoint operator of .

By Gram-Schmidt process we can construct an orthonormal system of , where

Theorem 1. *Let be a dense set in . Then is a complete system of .*

For the proof, see [14].

Theorem 2. *Let be a dense set in and the solution of (3) is unique on . Then the solution of (3) is given by , where .*

For the proof, see [14].

One can get an approximate solution by taking finitely many terms in the series representation of and .

Since is a Hilbert space, then .

Theorem 3. *The approximate solution and its derivatives are uniformly convergent to , , .*

*Proof. *By the reproducing kernel property of and Schwarz inequality, we can obtain
where is a constant.

By the representation of we can obtain
Since , , is uniformly bounded about and , we have
and so
Thus and its derivatives are uniformly convergent to .

#### 3. Numerical Results

In this paper, three numerical examples are given to show the accuracy of this method. The computations are performed by Mathematica 8.0. We compare the results by this method with the exact solution of each example.

*Example 1. *Consider the following linear system of fractional integrodifferential equations:
The exact solution for is , .

After homogenizing the initial conditions and using this method, taking , , and , the graphs of the approximate solutions for different values of and are plotted in Figure 1. From Figure 1, it is clear that the approximate solutions are in good agreement with the exact solutions when , and the solution continuously depends on the fractional derivative.

**(a)**

**(b)**

**(c)**

**(d)**

*Example 2. *Consider the following nonlinear system of fractional integrodifferential equations:
The exact solution for is , .

After homogenizing the initial conditions and using this method, taking , , and , the graphs of the approximate solutions for different values of and are plotted in Figure 2.

**(a)**

**(b)**

**(c)**

**(d)**

*Example 3. *Consider the following nonlinear system of fractional integrodifferential equations:
The exact solution for is , .

After homogenizing the initial conditions and using this method, taking , , and , the graphs of the approximate solutions for different values of and are plotted in Figure 3.

**(a)**

**(b)**

**(c)**

**(d)**

#### 4. Conclusion

In this paper, we introduce a new algorithm for solving systems of fractional integrodifferential equations. The approximate solution obtained by this method and its derivative are both uniformly convergent. The obtained results demonstrate the reliability of the algorithm and its wider applicability to linear and nonlinear systems of fractional differential equations.

#### Conflict of Interests

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