On a Fractional Master Equation

A fractional order time-independent form of the wave equation or diffusion equation in two dimensions is obtained from the standard time-independent form of the wave equation or diffusion equation in two-dimensions by replacing the integer order partial derivatives by fractional Riesz-Feller derivative and Caputo derivative of order α, β, 1 < R α ≤ 2 and 1 < R β ≤ 2 respectively. In this paper, we derive an analytic solution for the fractional time-independent form of the wave equation or diffusion equation in two dimensions in terms of the Mittag-Leffler function. The solutions to the fractional Poisson and the Laplace equations of the same kind are obtained, again represented bymeans of theMittag-Leffler function. In all three cases, the solutions are represented also in terms of Fox’sH-function.


Introduction
The standard time-independent form of the wave equation or diffusion equation in twodimensions where k > 0 is the wave number, is mathematically considered as the master equations to different classes of partial differential equations, namely, Poisson equation and Laplace equation.It represents the time-independent form of the wave equation or diffusion equation obtained while applying the technique of separation of variables to reduce the complexities of the solution procedure of the original equations.This equation appears in physical phenomena and engineering applications such as heat conduction, acoustic radiation, water wave propagation, and even in biology.For estimating the geodesic sea floor properties, the proper prediction of acoustic propagation in shallow water as well as at low frequencies is In electromagnetics, the two-dimensional time-independent form of the wave equation or diffusion equation appears as the governing equation for waveguide problems.There is huge mathematical and engineering interest in electromagnetic wave scattering problems driven by many applications such as modeling radar, sonar, acoustic noise barriers, atmospheric particle scattering, and ultrasound since both the incident and scattered electric field satisfy the two-dimensional time-independent form of the wave equation or diffusion equation which is also known as the scalar Helmholtz equation see Budiarto and Takada 3 .This paper introduces a new fractional-model time-independent form of the wave equation or diffusion equation in two-dimensions, in which both the space variables x and y are allowed to take fractional order changes.Such models are described in Mainardi et al. 4 and are defined as where k is the wave number given by k 2π/λ where λ is the wavelength, E x, y is the field variable of interest, which could be acoustic pressure, wave elevation, or electromagnetic potential, among many other possibilities, and Φ x, y is a nonlinear function in the field.x D α θ is the Riesz-Feller space fractional derivative of order α and asymmetry parameter skewness θ, and y D β * is the Caputo fractional derivative of order β.These fractional derivatives are integrodifferential operators, and are defined in the section on mathematical preliminaries.
The Mittag-Leffler function is a special function having an essential role in the solutions of fractional order integral and differential equations.Recently, this function is frequently used in modeling phenomena of fractional order appearing in physics, biology, engineering and applied sciences.After being introduced and studied by Mittag-Leffler 5 , Wiman 6 and Agarwal 7 , the Mittag-Leffler function, in its two forms: has been studied in details by Dzherbashyan 8 .Both functions 1.2 -1.3 are entire functions of order ρ 1/α and type σ 1.In 1920, Hille and Tamarkin 9 have presented a solution of the Abel-Volterra type integral equation The objective of this paper is to develop a solution of the fractional time-independent form of the wave equation or diffusion equation 1.1 in terms of Mittag-Leffler function and then in Fox's H-function, using the Laplace and Fourier transforms and their inverse transforms.Mathematically, the Poisson and the Laplace equations are the two special cases of the two-dimensional time-independent form of the wave equation or diffusion equation.We apply this fact to fractional case also.
This paper is divided as follows.Section 2 is devoted to mathematical preliminaries used to solve the Cauchy problems.In Section 3, we derive the solution of the fractional Laplace equation in Mittag-Leffler function and then in Fox's H-function.Section 4 is devoted to the fractional Poisson equation.The solution of the fractional master equation is given in Section 5. Its proof and the convergence and the series representation of the H-function are given in the appendix.The Mellin-Barnes representation and the series representation of the special functions in the fundamental solutions are given in the appendix.

Mathematical Preliminaries
The space fractional Riesz-Feller derivative x D α θ of order α and skewness θ is defined as where 0 < α ≤ 2, |θ| ≤ min α, 2 − α and the Caputo derivative of order α with respect to t, t > 0, and an empty product is always interpreted as unity; m, n, p, q ∈ N 0 with 0 , where we employ the usual notations: The Wright's generalized hypergeometric function, Wright 21,22 ; is defined by the series representation C is the complex number field.The Mittag-Leffler function is a special case of this function, The Laplace transform of the function N x, t with respect to t is and its inverse transform with respect to s is given by γ being a fixed real number.The Fourier transform of a function N x, t with respect to x is defined as the inverse Fourier transform with respect to p: The space of functions for the above two transforms is where

2.11
The Laplace transform of the Caputo fractional derivative see, e.g., Podlubny 16

2.12
From Mainardi et al. 4 , the Fourier transform of the Riesz-Feller derivative is given by where ψ θ α f p |p| α e i sign p θπ/2 , p ∈ R, and f F f x ; p ∞ −∞ e ipx f x dx.We also need the property of H-function Mathai et al. 11 2.14 and the following result Podlubny, 16 : 2.15

The Exact Solution of the Fractional Laplace Equation
The exact solution of the fractional Laplace equation 1.2 in terms of three special functions namely, Mittag-Leffler function, Fox's H-function, and Mainardi function, respectively, are derived in the following subsections.

The Exact Solution in Terms of Mittag-Leffler Function
Theorem 3.1.Consider the fractional Laplace equation x, y 0, where x D α θ is the Riesz-Feller fractional derivative of order α and symmetry θ and y D β * is the Caputo fractional derivative of order β.Then the solution of 3.1 subject to the initial conditions is where indicates the Fourier transform with respect to the space variable x.
Proof.If we apply Laplace transform with respect to the space variable y and use 2.12 , 3.1 becomes Now applying the Fourier transform with respect to the space variable x to obtain

3.5
Taking the inverse Fourier transform on both sides, we obtain the solution g p E β,2 ψ θ α p y β e −ipx dp.
From the above Cauchy problem, by introducing the initial conditions f x δ x where δ x is the Dirac delta function and g x 0, the exact solution or the Green function in Mittag-Leffler function is derived with the help of the Fourier convolution: where the exact solution in terms of the Mittag-leffler function G x, y is

Exact Solution in Terms of Fox's H-Function
From the above theorem, exact solution in Fox's H-function for the fractional Laplace equation where , using 2.14 .

3.11
International Journal of Differential Equations

Fractional Poisson Equation
As a generalization of the Cauchy problem 3.1, the Mittag-Leffler solution of the fractional Poisson equation is given as.

Cauchy Problem
The solution to the fractional Poisson equation x D α θ is the Riesz-Feller fractional derivative of order α and symmetry θ and y D β * is the Caputo fractional derivative of order β.Then, the solution of 4.2 subject to the initial conditions is where indicates the Fourier transform with respect to the space variable x.
Solution.Applying the Laplace transform with respect to the space variable y and Fourier transform with respect to the space variable x and using initial conditions, we have

4.4
Taking the inverse Laplace transform and using Laplace convolution,

4.5
International Journal of Differential Equations 9 Taking the inverse Fourier transform on both sides and using Fourier convolution, we obtain the solution 4.6

The Fractional Master Equation
In this section, we solve the fractional master equation, namely, the fractional timeindependent form of the wave or diffusion equation fractional Helmholtz equation 1.2 to generate a solution in terms of the Mittag-Leffler function, using the Laplace and Fourier transforms and their inverses.The wave number k is selected such a way that k 2 > ψ θ α p |p| α e i sign p θπ/2 , p ∈ R. Thus, the model is made suitable to handle the scattering problems of electromagnetic waves of large wave number or short wavelength.

5.2
Where is the Fourier transform with respect to x.The proof of this theorem and the H-function solution of the fractional time-independent form of the wave equation or diffusion equation are discussed in the appendix.

Conclusion
The closed form solutions in terms of the Mittag-Leffler function and If we apply the Laplace transform with respect to the space variable y and use 2.12 , 5.1 becomes Applying Fourier transform with respect to x, Thus, Taking Laplace inverse transform with respect to s, A.4 Using inverse Fourier transform with respect to p, the solution is It is seen that the above solution is a master solution of the solutions to the fractional Laplace equation 3.6 and the fractional Poisson equation 4.6 .

B. The Fox's H-Function Solution of the Fractional Time-Independent form of the Wave Equation or Diffusion Equation.
The Fox's H-function solution of the fractional time-independent form of the wave equation or diffusion equation is where the exact solution in terms of the H-function G θ α,β x, y is  where z −2 α e iθ π/2 |x| −α y β is given as follows.We have where z 2 α e iθ π/2 |x| −α y β .Let us assume that the poles of the integrand are simple.The region of convergence is L L ∞ is a loop beginning and ending at ∞ and encircling all the poles of Γ 1/2 − α/2 s once in the negative direction but none of the poles of Γ βs .By calculating the residues at the poles of Γ 1/2 − α/2 s where the poles are given by 1/2 − α/2 s −ν, ν 0, 1, 2, . .., we will get the series representation as

2 . 2 where
• • • is the integer part.The main results on Mittag-Leffler functions of 1.3 , 1.4 are available in the handbook of Erdélyi et al. 12, 13 , the monographs by Dzherbashyan 8, 14 , some recent books by Mathai et al. 11 , Kiryakova 15 , Podlubny 16 , Kilbas et al. 17 , and Mainardi 18 .The H-function is defined by means of a Mellin-Barnes type integral in the following manner Mathai and Haubold 19 : , ∞ ; C being the complex number field.For the details about the contour L and existence conditions see;Mathai et al. 11 and Kilbas and Saigo 20 .

Theorem 5 . 1 .
The analytical solution of the fractional time-independent form of the wave equation or diffusion equationx D α θ E x, y y Dβ * E x, y k 2 E x, y Φ x, y , 5.1 Fox's H-function are obtained for the fractional Laplace and the fractional Poisson equations.It is seen that the solutions in terms of the Mittag-Leffler function as well as in the H-function to International Journal of Differential Equations the fractional time-independent form of the wave equation or diffusion equation or the fractional Helmholtz equation are the master solutions to the solutions of the fractional Laplace equation and the fractional Poisson equation.Appendices A. The Mittag-Leffler Solution of the Fractional Time-Independent form of the Wave Equation or Diffusion Equation

Γ
βs Γ 1/2 − α/2 s Γ α/2s −z −s ds, D.1 International Journal of Differential Equations very essential.It also provides the solution to such problems, refer Liu et al. 1 .It also solves the problems in pattern formation in animal coating, see Murray and Myerscough 2 .
5in terms of Mittag-Leffler functions.Fox 10 used the H-function to give the most generalized symmetrical Fourier kernel.This function is defined in terms of Mellin-Barnes integrals and is a generalization of the Meijer G-function.Most of the special functions are available as the special cases of this functionMathai et al. 11 .
From Mathai et al. 11 , and Prudnikov et al. 23 the cosine transform of the H-function is defined as where L R is the space of summable functions on R with norm f ∞ the space of summable functions on R with norm f ∞ −∞ |f t |dt < ∞.