Research Article | Open Access

# Numerical Solution for Complex Systems of Fractional Order

**Academic Editor:**Turgut Öziş

#### Abstract

By using a complex transform, we impose a system of fractional order in the sense of Riemann-Liouville fractional operators. The analytic solution for this system is discussed. Here, we introduce a method of homotopy perturbation to obtain the approximate solutions. Moreover, applications are illustrated.

#### 1. Introduction

Fractional models have been studied by many researchers to sufficiently describe the operation of variety of computational, physical, and biological processes and systems. Accordingly, considerable attention has been paid to the solution of fractional differential equations, integral equations, and fractional partial differential equations of physical phenomena. Most of these fractional differential equations have analytic solutions, approximation, and numerical techniques [1–3]. Numerical and analytical methods have included finite difference methods such as Adomian decomposition method, variational iteration method, homotopy perturbation method, and homotopy analysis method [4–7].

The idea of the fractional calculus (i.e., calculus of integrals and derivatives of any arbitrary real or complex order) was planted over 300 years ago. Abel in 1823 investigated the generalized tautochrone problem and for the first time applied fractional calculus techniques in a physical problem. Later Liouville applied fractional calculus to problems in potential theory. Since that time the fractional calculus has drawn the attention of many researchers in all areas of sciences (see [8–10]).

One of the most frequently used tools in the theory of fractional calculus is furnished by the Riemann-Liouville operators. It possesses advantages of fast convergence, higher stability, and higher accuracy to derive different types of numerical algorithms. In this paper, we will deal with scalar linear time-space fractional differential equations. The time is taken in sense of the Riemann-Liouville fractional operators. Also, This type of differential equation arises in many interesting applications. For example, the Fokker-Planck partial differential equation, bond pricing equations, and the Black-Scholes equations are in this class of differential equations (partial and fractional).

In [11], the author used complex transform to obtain a system of fractional order (nonhomogeneous) keeping the equivalency properties. By employing the homotopy perturbation method, the analytic solution is presented for coupled system of fractional order. Furthermore, applications are imposed such as wave equations of fractional order.

#### 2. Fractional Calculus

This section concerns with some preliminaries and notations regarding the fractional calculus.

*Definition 2.1. *The fractional (arbitrary) order integral of the function of order is defined by
When , we write , where denoted the convolution product (see [12]), and and as where is the delta function.

*Definition 2.2. *The fractional (arbitrary) order derivative of the function of order is defined by

*Remark 2.3. *From Definitions 2.1 and 2.2, , we have
The Leibniz rule is
where

*Definition 2.4. *The Caputo fractional derivative of order is defined, for a smooth function by
where , (the notation stands for the largest integer not greater than ).

Note that there is a relationship between Riemann-Liouville differential operator and the Caputo operator
and they are equivalent in a physical problem (i.e., a problem which specifies the initial conditions).

In this paper, we consider the following fractional differential equation: where are complex valued functions, analytic in the domain and .

The above equation involves well-known time fractional diffusion equations.

#### 3. Complex Transforms

In this section, we will transform the fractional differential equation (2.8) into a coupled nonlinear system of fractional order has similar form. It was shown in [11] that the complex transform where is a complex valued function of complex variable , reduces (2.8) into the system where Also, it was shown that the complex transform reduces the nonhomogenous equation into the system where

#### 4. Numerical Solution

Let us put where and are arbitrary functions; are the linear parts of and , respectively. While and are the nonlinear parts of and , respectively. Moreover, let us set the homotopy system where Hence we obtain the following system: where Consequently, we have the approximate solution Thus, we impose a nonlinear integral equation in the following formula: Now we can sake the main result of this section.

Theorem 4.1. *Consider the fractional differential system (3.6) subject to the initial conditions
**The homotopy perturbation technique implies that the initial value problem ((3.6)–(4.9)) can be expressed as a nonlinear integral equation of the form (4.8).*

We proceed to prove the analytical convergence of our solution.

Theorem 4.2. *Suppose the sequence of the homotopy series and is defined for . Assume the initial approximation inside the domain of the solution . If for all , where , then the solution is absolutely convergent when .*

* Proof. *Let be the sequence of partial sum of the homotopy series. Our aim is to show that is a Cauchy sequence. Consider
For , we have
Hence
therefore, is a Cauchy sequence in the complex Banach space and consequently yields that the series solution is convergent. This completes the proof.

Recently the homotopy methods are used to obtain approximate analytic solutions of the time-fractional nonlinear equation and time-space-fractional nonlinear equation (see [12–17]).

#### 5. Applications

In this section, we will consider the pump wave equations along the fiber (Schrödinger equations). These types of equations are the fundamental equations for describing non-relativistic quantum mechanical behavior taking the form

Under the transform such that either or , we have the uncoupled system where . Subject to the initial conditions Operating (5.2) by , we have By the same computation as in Section 5, we receive Thus the solution is given by Moreover, under the same transform, (5.1) reduces to coupled system Operating (5.7) by , we have Therefore,

#### 6. Conclusion

We suggested two types of complex transforms for systems of fractional differential equations. We concluded that the complex fractional differential equations can be transformed into coupled and uncoupled system of homogeneous and nonhomogeneous types. Moreover, we employed the homotopy perturbation scheme for solving the nonlinear complex fractional differential systems. The convergence of the method is discussed in a domain that contains the initial solution. The Schrödinger equation is illustrated as an application. This type of equation is used in the quantum mechanics, which describes how the quantum state of a physical system changes with time. In the standard quantum mechanics, the wave function is the most complete explanation that can be specified to a physical system. Solutions of the Schrdinger’s equation describe not only molecular, atomic, and subatomic systems, but also macroscopic systems (see Figure 1).

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

#### References

- R. W. Ibrahim, “Existence and uniqueness of holomorphic solutions for fractional Cauchy problem,”
*Journal of Mathematical Analysis and Applications*, vol. 380, no. 1, pp. 232–240, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH - R. W. Ibrahim, “On holomorphic solution for space- and time-fractional telegraph equations in complex domain,”
*Journal of Function Spaces and Applications*, vol. 2012, Article ID 703681, 10 pages, 2012. View at: Google Scholar | Zentralblatt MATH - R. W. Ibrahim, “Approximate solutions for fractional differential equation in the unit disk,”
*Electronic Journal of Qualitative Theory of Differential Equations*, vol. 64, 11 pages, 2011. View at: Google Scholar - S. Momani and Z. Odibat, “Comparison between the homotopy perturbation method and the variational iteration method for linear fractional partial differential equations,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 910–919, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH - R. Y. Molliq, M. S. M. Noorani, and I. Hashim, “Variational iteration method for fractional heat- and wave-like equations,”
*Nonlinear Analysis. Real World Applications*, vol. 10, no. 3, pp. 1854–1869, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH - K. Sayevand, A. Golbabai, and A. Yildirim, “Analysis of differential equations of fractional order,”
*Applied Mathematical Modelling*, vol. 36, no. 9, pp. 4356–4364, 2012. View at: Publisher Site | Google Scholar - H. Jafari and M. A. Firoozjaee, “Homotopy analysis method for solving KdV equations,”
*Surveys in Mathematics and Its Applications*, vol. 5, pp. 89–98, 2010. View at: Google Scholar - I. Podlubny,
*Fractional Differential Equations*, Academic Press, London, UK, 1999. - A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo,
*Theory and applications of fractional differential equations*, North-Holland Mathematics Studies, Elsevier, 2006. - H. M. Srivastava and S. Owa,
*Univalent Functions, Fractional Calculus, and Their Applications*, Halsted Press, New York, NY, USA, 1989. - R. W. Ibrahim, “Complex transforms for systems of fractional differential equations,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 814759, 11 pages, 2012. View at: Publisher Site | Google Scholar - E. Hesameddini and H. Latifizadeh, “Homotopy analysis method to obtain numerical solutions of the Painlevé equations,”
*Mathematical Methods in the Applied Sciences*, vol. 35, no. 12, pp. 1423–1433, 2012. View at: Publisher Site | Google Scholar - R. W. Ibrahim and H. Jalab, “Analytic solution for fractional differential equation in the unit disk,”
*Wulfenia Journal*, vol. 19, no. 8, pp. 105–114, 2012. View at: Google Scholar - M. G. Sakar, F. Erdogan, and A. Yıldırım, “Variational iteration method for the time-fractional Fornberg-Whitham equation,”
*Computers & Mathematics with Applications*, vol. 63, no. 9, pp. 1382–1388, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH - A. Salah, M. Khan, and M. Gondal, “A novel solution procedure for fuzzy fractional heat equationsby homotopy analysis transform method,”
*Neural Computing & Applications*, vol. 2012, 3 pages, 2012. View at: Publisher Site | Google Scholar - A. El-Ajou, O. Abu Arqub, and S. Momani, “Homotopy analysis method for second-order boundary value problems of integrodifferential equations,”
*Discrete Dynamics in Nature and Society*, vol. 2012, Article ID 365792, 18 pages, 2012. View at: Google Scholar - K. Vishal and S. Das, “Solution of the nonlinear fractional diffusion equation with absorbent term and external force using optimal homotopy-analysis method,”
*Zeitschrift fr NAturforschung A*, vol. 67, pp. 203–209, 2012. View at: Google Scholar

#### Copyright

Copyright © 2012 Rabha W. Ibrahim. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.