Advanced
search

Exact solutions of (1+2)-dimensional non-linear time-space fractional PDEs

Manoj Kumar (Department of Mathematics, National Defence Academy, Pune, India)

Arab Journal of Mathematical Sciences

ISSN: 1319-5166

Article publication date: 7 June 2022

Issue publication date: 23 January 2024

506
pdf (154 KB) Article view Figure view Cited (10) cite article

Abstract

Purpose

In this paper, the author presents a hybrid method along with its error analysis to solve (1+2)-dimensional non-linear time-space fractional partial differential equations (FPDEs).

Design/methodology/approach

The proposed method is a combination of Sumudu transform and a semi-analytc technique Daftardar-Gejji and Jafari method (DGJM).

Findings

The author solves various non-trivial examples using the proposed method. Moreover, the author obtained the solutions either in exact form or in a series that converges to a closed-form solution. The proposed method is a very good tool to solve this type of equations.

Originality/value

The present work is original. To the best of the author's knowledge, this work is not done by anyone in the literature.

Keywords

Citation

Kumar, M. (2024), "Exact solutions of (1+2)-dimensional non-linear time-space fractional PDEs", Arab Journal of Mathematical Sciences, Vol. 30 No. 1, pp. 30-42. https://doi.org/10.1108/AJMS-11-2021-0282

Publisher

:

Emerald Publishing Limited

Copyright © 2022, Manoj Kumar

License

Published in the Arab Journal of Mathematical Sciences. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) license. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this license may be seen at http://creativecommons.org/licences/by/4.0/legalcode

1. Introduction

In the past few decades, Fractional Calculus has drawn the attention of many researchers due to its wide applicability in all disciplines. In contrast to ordinary derivatives, fractional derivatives are non-local in nature and carry the past information [3–6]. Differential equations of fractional orders have become an essential tool to understand real-life problems. It has been established that fractional order partial differential equations (FPDEs) provide an appropriate framework for the description of anomalous and non-Brownian diffusion. They are more effective while developing processes having memory effects [7,8]. Recently, various models such as fractal foam drainage model [9], Klein-Gordon model [10], a fractal model for the soliton motion [11] and so on in micro-gravity space have been discussed. The fractal Schrödinger system using the fractal derivatives has been studied in references [12,13].

Several analytic approximate and numerical methods have been developed to solve fractional differential equations and FPDEs in the literature. For solving linear-differential equations transform methods such as Laplace transform [14], Fourier transform [15], Mellin transform [16,17], fractional Fourier transform [18], natural transform [19], Sumudu transform [1], Elzaki transform [20], Differential transform [21], Jafari transform [22] and so on are useful. Further, several decomposition/iterative methods such as Adomian decomposition [23], homotopy perturbation [24], Daftardar-Gejji and Jafari [2], variational iteration [25] and so on are developed to solve linear/non-linear FPDEs that give solutions in terms of convergent series. Moreover, these methods have become popular as they do not involve discretization and are free from rounding off errors. Furthermore, various hybrid methods i.e. combinations of integral transform and decomposition/iterative methods such as iterative Laplace transform [26], fractional Laplace homotopy perturbation transform [27], homotopy perturbation Sumudu transform [28], Sumudu decomposition [29], Laplace decomposition [30], Laplace homotopy analysis [31], homotopy perturbation transform [32], etc. have proven to be quite effective. In 2016, Wang and Liu [33] introduced a hybrid method known as Sumudu transform iterative method (STIM). STIM is a combination of Sumudu transform and the Daftardar-Gejji and Jafari method (DGJM). Besides, it is observed that DGJM with Sumudu transform requires less computational time to solve non-linear fractional models as compared to the traditional methods, and gives more accuracy. In this paper, we extend STIM along with its error analysis for solving general (1+2)-dimensional time-space FPDEs. Moreover, the utility and efficiency of STIM is demonstrated by solving various non-trivial examples.

The organization of this paper is as follows: In section 2, the author gives some basic definitions, properties of fractional calculus and Sumudu transform. In section 3, the author develops STIM for (1+2)-dimensional time-space FPDEs, whereas its error-analysis is presented in section 4. In section 5, the author solves various non-linear (1+2)-dimensional time-space FPDEs. Finally the author draws conclusions in section 6.

2. Preliminaries

Useful definitions and properties of fractional calculus and Sumudu transform are presented here.

Definition 2.1.

[34] Mittag-Leffler function with one parameter μ is defined as

(1) Eμ(z)=∑k=0∞zkΓ(μk+1),Re(μ)>0,z∈C.
Whereas the Mittag-Leffler function with two parameters μ and ν is defined as [35]
(2) Eμ,ν(z)=∑k=0∞zkΓ(μk+ν),Re(μ)>0,z,ν∈C.

Definition 2.2.

[36] Caputo partial fractional derivative of order μ > 0 with respect to the variable t, of w(t, x, y) is defined as

∂μw(t,x,y)∂tμ=1Γ(n−μ)∫0t(t−τ)n−μ−1∂nw(τ,x,y)∂τndτ,n−1<μ<n,∂nw(t,x,y)∂tn,μ=n∈N.

Note:

  1. If μ = m, where m is a positive integer then

dμtqdtμ=0,ifq∈{0,1,2,…,m−1},Γ(q+1)Γ(q−μ+1)tq−μ,ifq∈Nandq≥m,orq∉Nandq>m−1.

  1. The Caputo derivative of one parameter Mittag-Leffler function is

    dμEμ(ktμ)dtμ=kEμ(ktμ),μ>0,k∈R.

Note that we consider the fractional partial derivatives ∂mνw(t,x,y)∂xmν and ∂pγw(t,x,y)∂ypγ, where m,p∈N as the sequential fractional derivatives [37,38] i.e.

∂mνw(t,x,y)∂xmν=∂ν∂xν∂ν∂xν…∂νw(t,x,y)∂xν︸m−times,∂pγw(t,x,y)∂ypγ=∂μ∂yμ∂μ∂yμ…∂μw(t,x,y)∂yμ︸p−times.
Definition 2.3.

[39] The Sumudu transform (ST) over the set of functions D={ϕ(t)|∃L,θk>0,k=1,2 such that|ϕ(t)|<Le|t|/θk if t∈(−1)k×[0,∞)} is defined as

S[ϕ(t)](ρ)=F(ρ)=∫0∞1ρe−tρϕ(t)dt=∫0∞e−tϕ(ρt)dt,ρ∈(−θ1,θ2).

The Sumudu transform of the power function tμ is

(3)S[tμ](ρ)=Γ(μ+1)ρμ,μ>−1.
Definition 2.4.

[39] The inverse Sumudu transform (IST) of F(ρ) is denoted by ϕ(t), and is defined by the following integral:

ϕ(t)=S−1[F(ρ)]=12πi∫z−i∞z+i∞1ρetρF(ρ)dρ,whereRe(1/ρ)>z,z∈C.

Inverse Sumudu transform of ρμ is

(4)S−1[ρμ]=tμΓ(μ+1),μ>−1.
Definition 2.5.

[40] Sumudu transform of a time-fractional Caputo derivative of order μ > 0 of a real valued function w(t, x, y), where m−1<μ≤m,m∈N is defined as

(5) S∂μw(t,x,y)∂tμ=ρ−μS[w(t,x,y)]−∑k=0m−1ρ−μ+k∂kw(0,x,y)∂tk.

3. Sumudu transform iterative method (STIM) for (1+2)-dimensional time and space FPDEs

In this section, the author extends STIM [33] to solve (1+2)-dimensional time-space FPDEs.

Consider the following general non-linear time-space FPDE in (1+2)-dimensional (l,m,n,p,q∈N):

(6)∂μw(t,x,y)∂tμ=Fx,y,w,∂νw∂xν,…,∂mνw∂xmν,∂γw∂yγ,…,∂pγw∂ypγ,∂kw(0,x,y)∂tk=hk(x,y),k=0,1,2,…,l−1,
where μ ∈ (l − 1, l], ν ∈ (n − 1, n], γ ∈ (q − 1, q] and w is a function of three variables t, x, y and Fx,y,w,∂νw∂xν,…,∂mνw∂xmν,∂γw∂yγ,…,∂pγw∂ypγ a known non-linear operator. Taking ST on both sides of eqn (6) and simplifying, we get
(7)S[w(t,x,y)]=∑k=0l−1ρk∂kw(0,x,y)∂tk+ρμSFx,y,w,∂νw∂xν,…,∂mνw∂xmν,∂γw∂yγ,…,∂pγw∂ypγ.

Taking IST of eqn (7), we get

(8)w(t,x,y)=S−1∑k=0l−1ρkhk(x,y)+S−1ρμSFx,y,w,∂νw∂xν,…,∂mνw∂xmν,∂γw∂yγ,…,∂pγw∂ypγ.

Eqn (8) is of the following form

(9)wt,x,y=gx,y+Nx,y,w,∂νw∂xν,…,∂mνw∂xmν,∂γw∂yγ,…,∂pγw∂ypγ,
where
(10)g(x,y)=S−1∑k=0l−1ρkhk(x,y),N=S−1ρμSFx,y,w,∂νw∂xν,…,∂mνw∂xmν,∂γw∂yγ,…,∂pγw∂ypγ,
in which g(x, y) is a known function and N a known non-linear operator. The author applies DGJM to solve eqn (9), in which the solution is expressed in terms of the following infinite series
(11)w(t,x,y)=∑i=0∞wi,
where wi are calculated recursively. According to DGJM the non-linear operator N is decomposed as
(12)N=Nx,y,w0,∂νw0∂xν,…,∂mνw0∂xmν,∂γw0∂yγ,…,∂pγw0∂ypγ+∑j=1∞Nx,y,∑i=0jwi,∂ν∑i=0jwi∂xν,…,∂mν∑i=0jwi∂xmν,∂γ∑i=0jwi∂yγ,…,∂pγ∑i=0jwi∂ypγ−∑j=1∞Nx,y,∑i=0j−1wi,∂ν∑i=0j−1wi∂xν,…,∂mν∑i=0j−1wi∂xmν,∂γ∑i=0j−1wi∂yγ,…,∂pγ∑i=0j−1wi∂ypγ.
(13)S−1ρμSFx,y,∑i=0∞wi,∂ν∑i=0∞wi∂xν,…,∂mν∑i=0∞wi∂xmν,∂γ∑i=0∞wi∂yγ,…,∂pγ∑i=0∞wi∂ypγ=S−1ρμSFx,y,w0,∂νw0∂xν,…,∂mνw0∂xmν,∂γw0∂yγ,…,∂pγw0∂ypγ+∑j=1∞S−1ρμSFx,y,∑i=0jwi,∂ν∑i=0jwi∂xν,…,∂mν∑i=0jwi∂xmν,∂γ∑i=0jwi∂yγ,…,∂pγ∑i=0jwi∂ypγ−∑j=1∞S−1ρμSFx,y,∑i=0j−1wi,∂ν∑i=0j−1wi∂xν,…,∂mν∑i=0j−1wi∂xmν,∂γ∑i=0j−1wi∂yγ,…,∂pγ∑i=0j−1wi∂ypγ.

Using eqns (11, 13) in eqn (9), we get

∑i=0∞wi=S−1∑k=0l−1ρkhk(x,y)+S−1ρμSFx,y,w0,∂νw0∂xν,…,∂mνw0∂xmν,∂γw0∂yγ,…,∂pγw0∂ypγ+∑j=1∞S−1ρμSFx,y,∑i=0jwi,∂ν∑i=0jwi∂xν,…,∂mν∑i=0jwi∂xmν,∑i=0jwi,∂γ∑i=0jwi∂yγ,…,∂pγ∑i=0jwi∂ypγ−S−1ρμSFx,y,∑i=0j−1wi,∂ν∑i=0j−1wi∂xν,…,∂mν∑i=0j−1wi∂xmν,∂γ∑i=0j−1wi∂yγ,…,∂pγ∑i=0j−1wi∂ypγ.

Thus, the author defines the recursive relation to calculate wi′s,i∈{0}∪N as follows:

(14)w0=S−1∑k=0l−1ρkhkx,y,w1=S−1ρμSFx,y,w0,∂νw0∂xν,…,∂mνw0∂xmν,∂γw0∂yγ,…,∂pγw0∂ypγ,wr+1=S−1ρμSFx,y,∑i=0rwi,∂ν∑i=0rwi∂xν,…,∂mν∑i=0rwi∂xmν,∂γ∑i=0rwi∂yγ,…,∂pγ∑i=0rwi∂ypγ−S−1ρμSFx,y,∑i=0r−1wi,∂ν∑i=0r−1wi∂xν,…,∂mν∑i=0r−1wi∂xmν,∂γ∑i=0r−1wi∂yγ,…,∂pγ∑i=0r−1wi∂ypγ, r≥1.
Note that the k − term STIM approximate solution of eqn (6) is: w(t, x, y) ≈ w0 + w1 + ⋯ + wk−1.

4. Error analysis

The author presents the error analysis of the proposed method by proving the following theorem.

Theorem 4.1.

Let N be a non-linear operator form a Hilbert space H→H and w(t, x, y) be the solution of eqn (6). Then the series ∑i=0∞wi converges to w(t, x, y) if ‖wm(t, x, y)‖ ≤ μ‖wm−1 (t, x, y)‖, ∀m∈N and 0 < μ < 1. Further, let sm=∑h=0m−1wh(t,x,y) be the m − term approximate solution of eqn (6). Then we have

w(t,x,y)−sm≤μm(1−μ)‖w0(t,x,y)‖.

Proof. It is clear that

wt,x,y−sm=wt,x,y−∑h=0m−1wht,x,y=∑h=m∞wht,x,y≤∑h=m∞‖wht,x,y‖≤∑h=m∞μh‖w0t,x,y‖≤μm+μm+1+⋯‖w0t,x,y‖≤μm1−μ‖w0t,x,y‖,
which is the required result. □

5. Illustrative examples

5.1 Time-space fractional non-linear Boussinesq equation

Consider the following non-linear time-space fractional (1+2)-dimensional Boussinesq equation [41]:

(15)∂μw∂tμ=∂ν∂xν(k1w+k2)∂ν(k1w+k2)∂xν+∂γ∂yγ(k1w+k2)∂γ(k1w+k2)∂yγ,
along with the initial condition
(16)w(0,x,y)=a+byγ,t>0,μ,ν,γ∈(0,1),
where k1, k2, a, b are arbitrary constants.

Taking ST on both sides of eqn (15) and using the property (5), we get

(17)S[w(t,x,y)]=w(0,x,y)+ρμS∂ν∂xν(k1w+k2)∂ν(k1w+k2)∂xν+∂γ∂yγ(k1w+k2)∂γ(k1w+k2)∂yγ.

Taking the IST on both sides of eqn (17), we obtain

(18)w(t,x,y)=S−1[w(0,x,y)]+S−1ρμS∂ν∂xν(k1w+k2)∂ν(k1w+k2)∂xν+∂γ∂yγ(k1w+k2)∂γ(k1w+k2)∂yγ.

Using the recurrence relation (14), we get

w0=w(0,x,y)=a+byγ,w1=S−1ρμS∂ν∂xν(k1w0+k2)∂ν(k1w0+k2)∂xν+∂γ∂yγ(k1w0+k2)∂γ(k1w0+k2)∂yγ,=bk12Γ(γ+1)2Γ(μ+1)tμ,wj=0,∀j≥2.

Hence, we obtain the following exact solution of (15-16)

w(t,x,y)=a+byγ+bk12Γ(γ+1)2Γ(μ+1)tμ.

Note that for a=95,b=e2 the STIM solution is same as obtained by invariant subspace method in reference [41].

5.2 Time-space fractional non-linear diffusion equation

Consider the following non-linear time-space fractional (1+2)-dimensional diffusion like FPDE:

(19)∂μw(t,x,y)∂tμ=12y2γ∂2νw(t,x,y)∂x2ν+x2ν∂2γw(t,x,y)∂y2γ,
along with the initial condition
(20)w(0,x,y)=y2γ,t>0,μ,ν,γ∈(0,1).

Taking ST and then IST on both sides of eqn (19), we obtain

(21)w(t,x,y)=S−1[w(0,x,y)]+S−1ρμS12y2γ∂2νw∂x2ν+x2ν∂2γw∂y2γ.

In view of the recurrence relation (14), we obtain

w0=S−1[w(0,x,y)]=y2γ,w1=S−1ρμS12y2γ∂2νw0∂x2ν+x2ν∂2γw0∂y2γ,=Γ(2γ+1)tμx2ν2Γ(μ+1),w2=S−1ρμS12y2γ∂2ν(w0+w1)∂x2ν+x2ν∂2γ(w0+w1)∂y2γ−S−1ρμS12y2γ∂2νw0∂x2ν+x2ν∂2γw0∂y2γ,=Γ(2ν+1)Γ(2γ+1)t2μy2γ4Γ(2μ+1),w3=Γ(2ν+1)Γ(2γ+1)2t3μx2ν8Γ(3μ+1),w4=Γ(2ν+1)2Γ(2γ+1)2t4μy2γ16Γ(4μ+1),w5=Γ(2ν+1)2Γ(2γ+1)3t5μx2ν32Γ(5μ+1),⋮

Hence, the series solution of (19-20) is

w(t,x,y)=y2γ+Γ(2γ+1)tμx2ν2Γ(μ+1)+Γ(2ν+1)Γ(2γ+1)t2μy2γ4Γ(2μ+1)+Γ(2ν+1)Γ(2γ+1)2t3μx2ν8Γ(3μ+1)+Γ(2ν+1)2Γ(2γ+1)2t4μy2γ16Γ(4μ+1)+Γ(2ν+1)2Γ(2γ+1)3t5μx2ν32Γ(5μ+1)+⋯,
which converges to
w(t,x,y)=λ1tμE2μ,μ+1(λt2μ)x2ν+1+λt2μE2μ,2μ+1(λt2μ)y2γ,
where λ1=Γ(2γ+1)2, λ2=Γ(2ν+1)2 and λ = λ1λ2.

5.3 Time-space fractional equation: population model-I

Consider the following time-space fractional generalized biological population model equation [42]:

(22)∂μw(t,x,y)∂tμ=∂2νw∂x2ν+∂2γw∂y2γ+hw(1−rw),t>0,μ,ν,γ∈(0,1),
along with the initial condition
(23)w(0,x,y)=Eνr8xνEγr8yγ.

Taking ST and IST of eqn (22), we get

wt,x,y=S−1w0,x,y+S−1ρμS∂2νw∂x2ν+∂2γw∂y2γ+hw1−rw.

Using the recurrence relation (14), we obtain

w0=S−1[w(0,x,y)]=Eνr8xνEγr8yγ,w1=S−1ρμS∂2νw0∂x2ν+∂2γw0∂x2γ+hw0(1−rw0)=Eνr8xνEγr8yγhtαΓ(α+1),w2=Eνr8xνEγr8yγh2t2αΓ(2α+1),w3=Eνr8xνEγr8yγh3t3αΓ(3α+1),⋮

Hence, the series solution of eqns (22-23) is

(24)w(t,x,y)=Eνr8xνEγr8yγ1+htαΓ(α+1)+h2t2αΓ(2α+1)+⋯,
which converges to the following closed form solution:
(25)w(t,x,y)=Eνr8xνEγr8yγEμ(htμ).

Note that for μ = ν = γ = 1, we have w(t,x,y)=er8(x+y+ht), which is the standard solution of the biological model given in reference [43].

5.4 Time-space fractional equation: population model-II

Consider the following non-linear time-space fractional order biological population model equation [44, 45]:

(26)∂μw(t,x,y)∂tμ=∂2νw2∂x2ν+∂2γw2∂x2γ+hw,t>0,μ,ν,γ∈(0,1),
along with the initial condition
(27)w(0,x,y)=k1+k2(xνyγ)12,h,k1,k2∈R.

Taking ST and IST of (26), we get

(28)w(t,x,y)=S−1[w(0,x,y)]+S−1ρμS∂2νw2∂x2ν+∂2γw2∂x2γ+hw.

In view of the recurrence relation (14), we get

w0=S−1[w(0,x,y)]=k1+k2(xνyγ)12,w1=S−1ρμS∂2νw02∂x2ν+∂2γw02∂x2γ+hw0=k1+k2(xνyγ)12htμΓ(μ+1),w2=k1+k2(xνyγ)12h2t2μΓ(2μ+1),w3=k1+k2(xνyγ)12h3t3μΓ(3μ+1),⋮

Hence, the series solution of (26-27) is

w(t,x,y)=k1+k2(xνyγ)121+htμΓ(μ+1)+h2t2μΓ(2μ+1)+h3t3μΓ(3μ+1)+⋯,,
which converges to the following closed form
w(t,x,y)=k1+k2(xνyγ)12Eμ(htμ).

6. Conclusions

The author presented a hybrid method along with its error-analysis for solving non-linear (1+2)-dimensional time-space FPDEs, where the derivatives are considered in Caputo sense. Moreover, the author solved various non-trivial examples of (1+2)-dimensional FPDEs to demonstrate the utility and efficiency of the proposed method. It has been observed that the proposed method is accurate and provides the solutions either in a convergent series or in exact form. Thus, the proposed method is a useful tool to solve non-linear (1+2)-dimensional time-space FPDEs.

Data Availability: The data that supports the findings of this study are available within the article.

References

[1]Watugala G. Sumudu transform: a new integral transform to solve differential equations and control engineering problems. Integrated Educ. 1993; 24(1): 35-43.

[2]Daftardar-Gejji V. and Jafari H. An iterative method for solving nonlinear functional equations. J Math Anal Appl. 2006; 316(2): 753-63.

[3]Atangana A., Secer A. A note on fractional order derivatives and table of fractional derivatives of some special functions. Abstract and applied analysis. Hindawi; 2013. 2013. 1-8.

[4]De Oliveira E.C., Tenreiro Machado J.A. A review of definitions for fractional derivatives and integral. Math Probl Eng. 2014; 2014: 1-6.

[5]Tarasov V.E. and Tarasova S.S. Fractional derivatives and integrals: what are they needed for?. Mathematics. 2020; 8(2): 164.

[6]Li C. and Zeng F. Numerical methods for fractional calculus, 24. CRC Press; 2015.

[7]Podlubny I. Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. Elsevier; 1998.

[8]Bhalekar S., Daftardar-Gejji V., Baleanu D., Magin R. Fractional Bloch equation with delay. Comput Mathematics Appl. 2011; 61(5): 1355-65.

[9]Wang K.-L. A study of the fractal foam drainage model in a microgravity space. Math Methods Appl Sci. 2021; 44: 10530-10540.

[10]Wang K. and Wang K. A new analysis for Klein-Gordon model with local fractional derivative. Alexandria Eng J. 2020; 59(5): 3309-13.

[11]Wang K. A new fractal model for the soliton motion in a microgravity space. Int J Numer Methods Heat Fluid Flow. 2020; 31: 442-451.

[12]Wang K. New variational theory for coupled nonlinear fractal Schrödinger system. Int J Numer Methods Heat Fluid Flow. 2021; 32: 589-597.

[13]Wang K.-L. A novel perspective for the fractal schrödinger equation. Fractals. 2021; 29(4): 2150093.

[14]Podlubny I. Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. Elsevier; 1998.

[15]Kemple, S. and Beyer, H. Global and causal solutions of fractional differential equations in: transform method and special functions. In: Varna96, Proceding of the 2nd International Workshop (SCTP). Singapore. 1997.

[16]Butera S. and Di Paola M. Fractional differential equations solved by using Mellin transform. Commun Nonlinear Sci Numer Simulation. 2014; 19(7): 2220-7.

[17]Klimek M. and Dziembowski D. On Mellin transform application to solution of fractional differential equations. Scientific Res Inst Mathematics Computer Sci. 2008; 7(2): 31-42.

[18]Almeida L.B. The fractional Fourier transform and time-frequency representations. IEEE Trans Signal Process. 1994; 42(11): 3084-91.

[19]Shah K., Khalil H., Khan R.A. Analytical solutions of fractional order diffusion equations by natural transform method. Iranian J Sci Technol Trans A: Sci. 2018; 42(3): 1479-90.

[20]Khalid M., Sultana M., Zaidi F., Arshad U. Application of Elzaki transform method on some fractional differential equations. Math Theor Model. 2015; 5(1): 89-96.

[21]Wang K.-L. and Wang K.-J. A modification of the reduced differential transform method for fractional calculus. Therm Sci. 2018; 22(4): 1871-5.

[22]Jafari H. A new general integral transform for solving integral equations. J Adv Res. 2021; 32: 133-8.

[23]Adomian G. Solving frontier problems of physics: the decomposition method. Kluwer Boston. 1994.

[24]He J.-H. Homotopy perturbation technique. Computer Methods Appl Mech Eng. 1999; 178(3): 257-62.

[25]He J.-H. Variational iteration method–a kind of non-linear analytical technique: some examples. Int J non-linear Mech. 1999; 34(4): 699-708.

[26]Jafari H., Nazari M., Baleanu D., Khalique C. A new approach for solving a system of fractional partial differential equations. Comput Mathematics Appl. 2013; 66(5): 838-43.

[27]Khan Y., Faraz N., Kumar S., Yildirim A. A coupling method of homotopy perturbation and Laplace transformation for fractional models. University “Politehnica” Bucharest Scientific Bull Ser A: Appl Mathematics Phys. 2012; 74(1): 57-68.

[28]Singh J. and Devendra S. Homotopy perturbation Sumudu transform method for nonlinear equations. Adv Theor Appl Mech. 2011; 4(4): 165-75.

[29]Kumar D., Singh J., Rathore S. Sumudu decomposition method for nonlinear equations. Int Math Forum. 2012; 7: 515-21.

[30]Jafari H., Yousefi S., Firoozjaee M., Momani S., Khalique C. M. Application of legendre wavelets for solving fractional differential equations. Comput Mathematics Appl. 2011; 62(3): 1038-45.

[31]Zurigat M. Solving fractional oscillators using Laplace homotopy analysis method. Ann Univ Craiova-Mathematics Computer Sci Ser. 2011; 38(4): 1-11.

[32]Singh B.K. and Kumar P. Homotopy perturbation transform method for solving fractional partial differential equations with proportional delay. SeMA J. 2018; 75(1): 111-25.

[33]Wang K. and Liu S. A new Sumudu transform iterative method for time-fractional Cauchy reaction–diffusion equation. SpringerPlus. 2016; 5(1): 865.

[34]Mittag-Leffler G. Sur la nouvelle fonction Eα (x). CR Acad Sci Paris. 1903; 137: 554-8.

[35]Wiman A. Uber den fundamentalsatz in der teorie der funktionen Eα(x). Acta Mathematica. 1905; 29(1): 191-201.

[36]Purohit S. Solutions of fractional partial differential equations of quantum mechanics. Adv Appl Math Mech. 2013; 5(5): 639-51.

[37]Miller K.S., Ross B. An introduction to the fractional calculus and fractional differential equations. New York: John-Wily and Sons; 1993.

[38]Chikrii A., Matychyn I.. Riemann–Liouville, Caputo, and sequential fractional derivatives in differential games. Advances in dynamic games. Birkhäuser Boston; 2011. 61-81.

[39]Belgacem F.B.M., Karaballi A.A.. Sumudu transform fundamental properties investigations and applications. Int J Stochastic Anal. 2006; May 2016: 1-23.

[40]Amer Y., Mahdy A., Youssef E. Solving systems of fractional nonlinear equations of Emden Fowler type by using Sumudu transform method. Glob J Pure Appl Mathematics. 2018; 14(1): 91-113.

[41]Choudhary S. and Daftardar-Gejji V. Solving systems of multi-term fractional PDEs: invariant subspace approach. Int J Model Simulation, Scientific Comput. 2019; 10(1): 1941010.

[42]Shakeri F. and Dehghan M. Numerical solution of a biological population model using he's variational iteration method. Comput Mathematics Appl. 2007; 54(7-8): 1197-209.

[43]Liu Y., Li Z., Zhang Y. Homotopy perturbation method to fractional biological population equation. Fractional Differential Calculus. 2011; 1(1): 117-24.

[44]El-Sayed A., Rida S., Arafa A. Exact solutions of fractional-order biological population model. Commun Theor Phys. 2009; 52(6): 992.

[45]Gurtin M.E. and MacCamy R.C. On the diffusion of biological populations. Math Biosciences. 1977; 33(1-2): 35-49.

Corresponding author

Manoj Kumar can be contacted at: mkmath14@gmail.com

Related articles

Manage cookies

We’re listening — tell us what you think

Join us on our journey