# On a Fractional Master Equation

Authors: Anitha Thomas
Journal: International Journal of Differential Equations
 (2011)

---

## Abstract

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<ℜ(α)≤2 and 1<ℜ(β)≤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 by means of the Mittag-Leffler function. In all three cases, the solutions are represented also in terms of Fox's H-function.


---

## Body

## 1. Introduction

The standard time-independent form of the wave equation or diffusion equation in two-dimensions(1.1)∇2Ψ(x,y)+k2Ψ(x,y)=0,
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 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].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
(1.2)xDθαE(x,y)+0DyβE(x,y)+k2E(x,y)=Φ(x,y)k>0, x∈ℜ, y∈ℜ+, 1<ℜ(α)≤2, 1<ℜ(β)≤2, 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. xDθα is the Riesz-Feller space fractional derivative of order α and asymmetry parameter (skewness) θ, and yD*β 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:
(1.3)Eα(z)=∑n=0∞znΓ(nα+1),α∈C,R(α)>0,(1.4)Eα,β(z)=∑n=0∞znΓ(nα+β),α,β∈C,R(α)>0,R(β)>0,
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
(1.5)ϕ(x)-λΓ(α)∫0xϕ(t)(x-t)1-αdt=f(x),0<x<1,
in 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 function (Mathai et al. [11]).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. Section2 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.


## 2. Mathematical Preliminaries

The space fractional Riesz-Feller derivativexDθα of order α and skewness θ is defined as
(2.1)Dxθαf(x)=Γ(1+α)π[sin⁡[(α+θ)π2]∫0∞f(x+ξ)-f(x)ξ1+αdξ]+Γ(1+α)π[sin⁡[(α-θ)π2]∫0∞f(x-ξ)-f(x)ξ1+αdξ],
where 0<α≤2,|θ|≤min(α,2-α) and the Caputo derivative of order α with respect to t,t>0,(2.2)tD*αN(x,t)={1Γ(m-α)∫0tN(m)(x,u)(t-u)α-m+1du,m-1<R(α)<m,dmdtmN(x,t),α=m,m∈N,
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]):
(2.3)Hp,qm,n(z)[z|(bq,Bq)(ap,Ap)]=12πi∫L{∏j=1mΓ(bj+Bjs)}{∏j=1nΓ(1-bj-Bjs)}{∏j=m+1qΓ(1-bj-Bjs)}{∏j=n+1pΓ(aj+Ajs)}z-sds,
and an empty product is always interpreted as unity; m,n,p,q∈N0 with 0≤n≤p, 1≤m≤q, Aj,Bj∈ℜ+,ai,bj∈ℂ, i=1,…,p;  j=1,…,q such that Ai(bj+k)≠Bj(ai-l-1), k,l∈N0; i=1,…,n; j=1,…,m,  where we employ the usual notations: N0=(0,1,…,); ℜ=(-∞,∞),ℜ+=(0,∞); ℂ 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].The Wright's generalized hypergeometric function, Wright ([21, 22]); is defined by the series representation
(2.4)pψq(z)=pψq[z|(bq,Bq)(ap,Ap)]=∑r=0∞[∏j=1pΓ(aj+Ajs)]zr[∏j=1qΓ(bj+Bjs)]r!,
where z∈ℂ,aj,bj∈ℂ,Aj,Bj∈ℜ+;i=1,…,p;j=1,…,q;∑j=1qbj-∑j=1pAj>-1;ℂ is the complex number field. The Mittag-Leffler function is a special case of this function,
(2.5)Eα,β(z)=1ψ1[z|(β,α)(1,1)]=H1,21,1[-z|(0,1)(1-β,α)(0,1)].
The Laplace transform of the function N(x,t) with respect to t is
(2.6)L[N(x,t)]=∫0∞e-stN(x,t)dt=N*(x,s),x∈R,R(s)>0,
and its inverse transform with respect to s is given by
(2.7)L-1[N*(x,s)]=12πi∫γ-i∞γ+i∞estN*(x,s)ds=N(x,t),γ being a fixed real number. The Fourier transform of a function N(x,t) with respect to x is defined as
(2.8)F[N(x,t)]=∫-∞∞eipxN(x,t)dx=N̂(p,t),x∈R,
the inverse Fourier transform with respect to p:
(2.9)F-1[N̂(p,t)]=12π∫-∞∞e-ipxN̂(p,t)dp=N(x,t).
The space of functions for the above two transforms is ℒℱ=ℒ(ℜ+)×ℱ(ℜ), where ℒ(ℜ+) is the space of summable functions on ℜ+ with norm ∥f∥=∫0∞|f(t)|dt<∞ and ℱ(ℜ) is the space of summable functions on ℜ with norm ∥f∥=∫-∞∞|f(t)|dt<∞. From Mathai et al. [11], and Prudnikov et al. [23] the cosine transform of the H-function is defined as
(2.10)∫0∞tρ-1cos⁡(kt)Hp,qm,n[atμ|(bq,Bq)(ap,Ap)]dt=2ρ-1πkρHp+2,qm,n+1[a(2k)μ|(bq,Bq)(2-ρ/2,μ/2),(ap,Ap),((1-ρ)/2,μ/2)],
where  (2.11)R(ρ)+μmin⁡1≤j≤m⁡R[bjBj]>0,R(ρ)+μmax⁡1≤j≤n⁡[(aj-1)Aj]<1,|arga|<12πθ,θ=∑j=1nAj-∑j=n+1pAj+∑j=1mBj-∑j=m+1qBj>0.
The Laplace transform of the Caputo fractional derivative (see, e.g., Podlubny [16])
(2.12)L[D0tαN(x,t)]=sαN*(x,s)-∑r=1nsr-10Dtα-rN(x,t)∣t=0,n-1<R(α)≤n,n∈N.
From Mainardi et al. [4], the Fourier transform of the Riesz-Feller derivative is given by
(2.13)F{Dxθαf(x);p}=-ψαθf̂(p),
where ψαθf̂(p)=|p|αei(signp)θπ/2,p∈ℜ,  and f̂=ℱ[f(x);p]=∫-∞+∞eipxf(x)dx.We also need the property ofH-function (Mathai et al. [11])
(2.14)Hp,qm,n[z|(b1,B1),…(bq-1,Bq-1),(a1,A1)(a1,A1),…(ap,Ap)]=Hp-1,q-1m,n-1[z|(b1,B1)…(bq-1,Bq-1)(a2,A2),…(ap,Ap)],m≥1,q>m,
and the following result (Podlubny, [16]):
(2.15)L-1(sβ-1sα-a;t)=tα-βEα,α-β+1(atα),R(s)>0,R(α-β)>-1,|asα|<1.


## 3. 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.


### 3.1. The Exact Solution in Terms of Mittag-Leffler Function

Theorem 3.1.
Consider the fractional Laplace equation(3.1)xDθαE(x,y)+yD*βE(x,y)=0,x∈R,y⩾0,1<ℜ(α)≤2, |θ|≤min(α,2-α), 1<ℜ(β)≤2 with initial conditions E(x,0)=f(x), ∂/∂yE(x,y)∣y=0=g(x), x∈ℜ, limx→±∞E(x,y)=0, where xDθα is the Riesz-Feller fractional derivative of order α and symmetry θ and yD*β is the Caputo fractional derivative of order β. Then the solution of (3.1) subject to the initial conditions is
(3.2)E(x,y)=12π∫-∞∞f̂(p)Eβ,1[(ψαθ(p))yβ]e-ipxdp+y2π∫-∞∞ĝ(p)Eβ,2[(ψαθ(p))yβ]e-ipxdp,
where ̂ indicates the Fourier transform with respect to the space variable x.Proof.
If we apply Laplace transform with respect to the space variabley and use (2.12), (3.1) becomes
(3.3)xDθαE*(x,s)+sβE*(x,s)-sβ-1f(x)-sβ-2g(x)=0.
Now applying the Fourier transform with respect to the space variable x to obtain
(3.4)[sβ-ψαθ(p)]Ê*(p,s)=sβ-1f̂(p)+sβ-2ĝ(p),Ê*(p,s)=sβ-1f̂(p)sβ-ψαθ(p)+sβ-2ĝ(p)sβ-ψαθ(p).
Using the result (2.15), we get
(3.5)Ê(p,y)=f̂(p)Eβ,1[ψαθ(p)yβ]+ĝ(p)yEβ,2[ψαθ(p)yβ].
Taking the inverse Fourier transform on both sides, we obtain the solution
(3.6)E(x,y)=12π∫-∞∞f̂(p)Eβ,1[(ψαθ(p))yβ]e-ipxdp+y2π∫-∞∞ĝ(p)Eβ,2[(ψαθ(p))yβ]e-ipxdp.
Hence Theorem 3.1 follows.
From the above Cauchy problem, by introducing the initial conditionsf(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:
(3.7)E(x,y)=∫-∞∞Gα,βθ(x-τ,y)f(τ)dτ,
where the exact solution in terms of the Mittag-leffler function G(x,y) is
(3.8)Gα,βθ(x,y)=12π∫-∞∞Eβ,1[(ψαθ(p))yβ]e-ipxdp.


### 3.2. Exact Solution in Terms of Fox'sH-Function

From the above theorem, exact solution in Fox'sH-function for the fractional Laplace equation
(3.9)xDθαE(x,y)+0DyβE(x,y)=0,x∈R,y⩾0,1<ℜ(α)≤2,  |θ|≤min(α,2-α),  1<ℜ(β)≤2 with initial conditions E(x,0)=f(x)=δ(x),∂/∂yE(x,y)∣y=0=0, x∈ℜ, limx→±∞E(x,y)=0, is given by
(3.10)E(x,y)=∫-∞∞Gα,βθ(x-τ,y)f(τ)dτ,
where (3.11)Gα,βθ(x,y)=12π∫-∞∞Eβ,1[(ψαθ(p))yβ]e-ipxdp,=1π∫0∞cos⁡(px)H1,21,1[(-ψαθ(p))yβ|(0,1),(0,β)(0,1)]dp,=1π∫0∞cos⁡(px)H1,21,1[-(|p|αyβeiθ(π/2))|(0,1),(0,β)(0,1)]dp,=1π|x|H3,21,2[-2αeiθ(π/2)|x|-αyβ|(0,1),(0,β)(1/2,α/2),(0,1),(0,α/2)],R(β)<2,=1π|x|H2,11,1[-2αeiθ(π/2)|x|-αyβ|(0,β)(1/2,α/2),(0,α/2)],using(2.14).


## 3.1. The Exact Solution in Terms of Mittag-Leffler Function

Theorem 3.1.
Consider the fractional Laplace equation(3.1)xDθαE(x,y)+yD*βE(x,y)=0,x∈R,y⩾0,1<ℜ(α)≤2, |θ|≤min(α,2-α), 1<ℜ(β)≤2 with initial conditions E(x,0)=f(x), ∂/∂yE(x,y)∣y=0=g(x), x∈ℜ, limx→±∞E(x,y)=0, where xDθα is the Riesz-Feller fractional derivative of order α and symmetry θ and yD*β is the Caputo fractional derivative of order β. Then the solution of (3.1) subject to the initial conditions is
(3.2)E(x,y)=12π∫-∞∞f̂(p)Eβ,1[(ψαθ(p))yβ]e-ipxdp+y2π∫-∞∞ĝ(p)Eβ,2[(ψαθ(p))yβ]e-ipxdp,
where ̂ indicates the Fourier transform with respect to the space variable x.Proof.
If we apply Laplace transform with respect to the space variabley and use (2.12), (3.1) becomes
(3.3)xDθαE*(x,s)+sβE*(x,s)-sβ-1f(x)-sβ-2g(x)=0.
Now applying the Fourier transform with respect to the space variable x to obtain
(3.4)[sβ-ψαθ(p)]Ê*(p,s)=sβ-1f̂(p)+sβ-2ĝ(p),Ê*(p,s)=sβ-1f̂(p)sβ-ψαθ(p)+sβ-2ĝ(p)sβ-ψαθ(p).
Using the result (2.15), we get
(3.5)Ê(p,y)=f̂(p)Eβ,1[ψαθ(p)yβ]+ĝ(p)yEβ,2[ψαθ(p)yβ].
Taking the inverse Fourier transform on both sides, we obtain the solution
(3.6)E(x,y)=12π∫-∞∞f̂(p)Eβ,1[(ψαθ(p))yβ]e-ipxdp+y2π∫-∞∞ĝ(p)Eβ,2[(ψαθ(p))yβ]e-ipxdp.
Hence Theorem 3.1 follows.
From the above Cauchy problem, by introducing the initial conditionsf(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:
(3.7)E(x,y)=∫-∞∞Gα,βθ(x-τ,y)f(τ)dτ,
where the exact solution in terms of the Mittag-leffler function G(x,y) is
(3.8)Gα,βθ(x,y)=12π∫-∞∞Eβ,1[(ψαθ(p))yβ]e-ipxdp.


## 3.2. Exact Solution in Terms of Fox'sH-Function

From the above theorem, exact solution in Fox'sH-function for the fractional Laplace equation
(3.9)xDθαE(x,y)+0DyβE(x,y)=0,x∈R,y⩾0,1<ℜ(α)≤2,  |θ|≤min(α,2-α),  1<ℜ(β)≤2 with initial conditions E(x,0)=f(x)=δ(x),∂/∂yE(x,y)∣y=0=0, x∈ℜ, limx→±∞E(x,y)=0, is given by
(3.10)E(x,y)=∫-∞∞Gα,βθ(x-τ,y)f(τ)dτ,
where (3.11)Gα,βθ(x,y)=12π∫-∞∞Eβ,1[(ψαθ(p))yβ]e-ipxdp,=1π∫0∞cos⁡(px)H1,21,1[(-ψαθ(p))yβ|(0,1),(0,β)(0,1)]dp,=1π∫0∞cos⁡(px)H1,21,1[-(|p|αyβeiθ(π/2))|(0,1),(0,β)(0,1)]dp,=1π|x|H3,21,2[-2αeiθ(π/2)|x|-αyβ|(0,1),(0,β)(1/2,α/2),(0,1),(0,α/2)],R(β)<2,=1π|x|H2,11,1[-2αeiθ(π/2)|x|-αyβ|(0,β)(1/2,α/2),(0,α/2)],using(2.14).


## 4. Fractional Poisson Equation

As a generalization of the Cauchy problem 3.1, the Mittag-Leffler solution of the fractional Poisson equation(4.1)xDθαE(x,y)+0DyβE(x,y)=Φ(x,y)
is given as.


### 4.1. Cauchy Problem

The solution to the fractional Poisson equation(4.2)xDθαE(x,y)+0DyβE(x,y)=Φ(x,y),y⩾0,x∈R,1<R(α)≤2,
where Φ(x,y) is a nonlinear function, 1<ℜ(α)≤2, |θ|≤min(α,2-α), 1<ℜ(β)≤2 with initial conditions E(x,0)=f(x), ∂/∂yE(x,y)∣y=0=g(x), x∈ℜ, limx→±∞E(x,y)=0 where xDθα is the Riesz-Feller fractional derivative of order α and symmetry θ and yD*β is the Caputo fractional derivative of order β. Then, the solution of (4.2) subject to the initial conditions is
(4.3)E(x,y)=12π∫-∞∞f̂(p)Eβ,1[(ψαθ(p))yβ]e-ipxdp+y2π∫-∞∞ĝ(p)Eβ,2[(ψαθ(p))yβ]e-ipxdp+12π∫0yξβ-1∫-∞∞Φ̂(p,y-ξ)Eβ,β[(ψαθ(p))ξβ]e-ipxdpdξ,
where ̂ indicates the Fourier transform with respect to the space variable x.Solution 4.1.
Applying the Laplace transform with respect to the space variabley and Fourier transform with respect to the space variable x and using initial conditions, we have
(4.4)[sβ-ψαθ(p)]Ê*(p,s)=sβ-1f̂(p)+sβ-2ĝ(p)+Φ̂*(p,s),Ê*(p,s)=sβ-1f̂(p)sβ-ψαθ(p)+sβ-2ĝ(p)sβ-ψαθ(p)+Φ̂*(p,s)sβ-ψαθ(p).
Taking the inverse Laplace transform and using Laplace convolution,
(4.5)Ê(p,y)=f̂(p)Eβ,1[(ψαθ(p))yβ]+ĝ(p)yEβ,2[ψαθ(p)ξβ]+∫0∞Φ̂(p,y-ξ)Eβ,β[(ψαθ(p))yβ]dξ.
Taking the inverse Fourier transform on both sides and using Fourier convolution, we obtain the solution
(4.6)E(x,y)=12π∫-∞∞f̂(p)Eβ,1[(ψαθ(p))yβ]e-ipxdp+y2π∫-∞∞ĝ(p)Eβ,2[(ψαθ(p))yβ]e-ipxdp+12π∫0yξβ-1∫-∞∞Φ̂(p,y-ξ)Eβ,β[(ψαθ(p))ξβ]e-ipxdpdξ.


## 4.1. Cauchy Problem

The solution to the fractional Poisson equation(4.2)xDθαE(x,y)+0DyβE(x,y)=Φ(x,y),y⩾0,x∈R,1<R(α)≤2,
where Φ(x,y) is a nonlinear function, 1<ℜ(α)≤2, |θ|≤min(α,2-α), 1<ℜ(β)≤2 with initial conditions E(x,0)=f(x), ∂/∂yE(x,y)∣y=0=g(x), x∈ℜ, limx→±∞E(x,y)=0 where xDθα is the Riesz-Feller fractional derivative of order α and symmetry θ and yD*β is the Caputo fractional derivative of order β. Then, the solution of (4.2) subject to the initial conditions is
(4.3)E(x,y)=12π∫-∞∞f̂(p)Eβ,1[(ψαθ(p))yβ]e-ipxdp+y2π∫-∞∞ĝ(p)Eβ,2[(ψαθ(p))yβ]e-ipxdp+12π∫0yξβ-1∫-∞∞Φ̂(p,y-ξ)Eβ,β[(ψαθ(p))ξβ]e-ipxdpdξ,
where ̂ indicates the Fourier transform with respect to the space variable x.Solution 4.1.
Applying the Laplace transform with respect to the space variabley and Fourier transform with respect to the space variable x and using initial conditions, we have
(4.4)[sβ-ψαθ(p)]Ê*(p,s)=sβ-1f̂(p)+sβ-2ĝ(p)+Φ̂*(p,s),Ê*(p,s)=sβ-1f̂(p)sβ-ψαθ(p)+sβ-2ĝ(p)sβ-ψαθ(p)+Φ̂*(p,s)sβ-ψαθ(p).
Taking the inverse Laplace transform and using Laplace convolution,
(4.5)Ê(p,y)=f̂(p)Eβ,1[(ψαθ(p))yβ]+ĝ(p)yEβ,2[ψαθ(p)ξβ]+∫0∞Φ̂(p,y-ξ)Eβ,β[(ψαθ(p))yβ]dξ.
Taking the inverse Fourier transform on both sides and using Fourier convolution, we obtain the solution
(4.6)E(x,y)=12π∫-∞∞f̂(p)Eβ,1[(ψαθ(p))yβ]e-ipxdp+y2π∫-∞∞ĝ(p)Eβ,2[(ψαθ(p))yβ]e-ipxdp+12π∫0yξβ-1∫-∞∞Φ̂(p,y-ξ)Eβ,β[(ψαθ(p))ξβ]e-ipxdpdξ.


## 5. The Fractional Master Equation

In this section, we solve the fractional master equation, namely, the fractional time-independent 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 k2>ψαθ(p)=|p|αei(signp)θπ/2,p∈ℜ. Thus, the model is made suitable to handle the scattering problems of electromagnetic waves of large wave number or short wavelength.Theorem 5.1.
The analytical solution of the fractional time-independent form of the wave equation or diffusion equation(5.1)xDθαE(x,y)+yD*βE(x,y)+k2E(x,y)=Φ(x,y),k>0, x∈ℜ, y∈ℜ+, with the initial conditions E(x,0)=f(x), ∂/∂yE(x,y)∣y=0=g(x), limx→±∞E(x,y)=0, is
(5.2)E(x,y)=12π∫-∞∞f̂(p)Eβ,1[-(k2-ψαθ(p))yβ]exp⁡(-ipx)dp+y2π∫-∞∞ĝ(p)Eβ,2[-(k2-ψαθ(p))yβ]exp⁡(-ipx)dp+12π∫0∞ξβ-1∫-∞∞Φ̂(p,y-ξ)Eβ,β[-(k2-ψαθ(p))ξβ]exp⁡(-ipx)dpdξ,
Where ̂ is the Fourier transform with respect to x.
The proof of this theorem and theH-function solution of the fractional time-independent form of the wave equation or diffusion equation are discussed in the appendix.


## 6. Conclusion

The closed form solutions in terms of the Mittag-Leffler function and Fox'sH-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 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.



---

*Source: 346298-2011-10-30.xml*