Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina


IMPACT FACTOR 2018: 1.005

CiteScore 2018: 1.01

SCImago Journal Rank (SJR) 2018: 0.237
Source Normalized Impact per Paper (SNIP) 2018: 0.541

ICV 2017: 162.45

Open Access
Online
ISSN
2391-5471
See all formats and pricing
More options …
Volume 16, Issue 1

Issues

Volume 13 (2015)

The Greek parameters of a continuous arithmetic Asian option pricing model via Laplace Adomian decomposition method

Sunday O. Edeki
  • Corresponding author
  • Department of Mathematics, Covenant University, KM. 10 Idiroko Road, Canaan Land, Ota, Nigeria
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Tanki Motsepa
  • International Institute for Symmetry Analysis and Mathematical Modelling, Department of Mathematical Sciences, North-West University, Private Bag X2046, Mmabatho 2735, South Africa
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Chaudry Masood Khalique
  • International Institute for Symmetry Analysis and Mathematical Modelling, Department of Mathematical Sciences, North-West University, Private Bag X2046, Mmabatho 2735, South Africa
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Grace O. Akinlabi
  • Department of Mathematics, Covenant University, KM. 10 Idiroko Road, Canaan Land, Ota, Nigeria
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2018-12-14 | DOI: https://doi.org/10.1515/phys-2018-0097

Abstract

The Greek parameters in option pricing are derivatives used in hedging against option risks. In this paper, the Greeks of the continuous arithmetic Asian option pricing model are derived. The derivation is based on the analytical solution of the continuous arithmetic Asian option model obtained via a proposed semi-analytical method referred to as Laplace-Adomian decomposition method (LADM). The LADM gives the solution in explicit form with few iterations. The computational work involved is less. Nonetheless, high level of accuracy is not neglected. The obtained analytical solutions are in good agreement with those of Rogers & Shi (J. of Applied Probability 32: 1995, 1077-1088), and Elshegmani & Ahmad (ScienceAsia, 39S: 2013, 67–69). The proposed method is highly recommended for analytical solution of other forms of Asian option pricing models such as the geometric put and call options, even in their time-fractional forms. The basic Greeks obtained are the Theta, Delta, Speed, and Gamma which will be of great help to financial practitioners and traders in terms of hedging and strategy.

Keywords: Option pricing; Black-Scholes model; Adomian decomposition method; Laplace transform

PACS: 89.65.Gh; 05.10.Gg; 44.05.+e

1 Introduction

In financial mathematics, the Greeks also referred to as sensitivity parameters are partial derivatives of the option prices with respect to some fundamental parameters. These Greeks are of great interest for hedging and risk management [1, 2]. Different dimensions to the risk associated with an option position are measured accordingly by different and unique Greeks. The following basic Greeks: Delta, Speed, Theta, and Gamma are studied with respect to (w.r.t.) Asian option while their associated mathematical expressions follow in the later part of this paper. Asian option is a special form of option contract whose value is hinged on the average value of the associated underlying asset over the option life time. Asian options are path dependent in nature unlike other options such as the European, American, lookback options and so on [3, 4, 5, 6, 7, 8, 9, 10]. Basically, Asian options are of two kinds viz: geometric Asian option (GAO) and a rithmetic Asian option (AAO). The GAO is noted to have a closed form solution. However, the AAO is difficult to price in terms of closed form solution [11, 12]. Hence, many researchers have developed solution techniques to that effect [13, 14, 15, 16, 17, 18, 19, 20,]. Other numerical methods that are of interest are [21, 22, 23, 24, 25, 26, 27, 28, 29, 30]. Some vital research approaches involving neural networks in relation to stochastic differential equations (SDEs) and/or option pricing are captured in [31, 32, 33, 34, 35, 36,].

In this paper, we propose the Laplace-Adomian decomposition method (LADM) for the first time in literature as a semi-analytical method, for analytical solution of a continuous arithmetic Asian option pricing model. Thereafter, the basic Greeks of the AAO model are explicitly derived. The remaining parts of the paper are organized as follows: in section 2, a brief note on the Asian option pricing model is given. In section 3, the proposed solution method (LADM) is presented, section 4 contains the application and the Greek-terms while in section 5, concluding remark is made.

2 Asian option pricing model

The stock price S (t) at time t, is assumed to satisfy a geometric Brownian motion (GBM) governed by the stochastic dynamic:

dSt=Strdt+σdWt,tR+(1)

where σ is a volatility coefficient, r a drift term indicating average rate of growth, and W (t) , t ∈ [0, T] a standard Brownian motion. The payoff for an Asian option with arithmetic average strike is given as:

ΞT=maxST1T0TSςdς,0.(2)

Note, the option price at t ∈ [0, T]is a risk neutral pricing formula defined as [3]:

Ξt=EerTtΞTFt(3)

whereE (·) and F t denote mathematical expectation operator and filtration respectively.

The payoff, ΞTis path-dependent. Hence, the introduction of a stochastic process [4]:

It=0TSςdςdIt=Stdt,S0=S0,(SDEform)(4)

where I (t) is the running sum of the strike price. Therefore, the corresponding Asian call option price is characterized by the model:

Ξt+12σ2S22ΞS2+rSΞS+SΞIrΞ=0(5)

which is satisfied by ΞS,I,tfor continuous arithmetic average strike such that t ≥ 0, and S > 0. Equation (5) is similar to the classical time-fractional Black-Scholes model at α = 1 in [21, 22] except for the averaging term SΞI.

This may later call for modification while numerical or semi-analytical methods are adopted [23, 24]. It is obvious that (5) will eventually lead to a greater computational problem because of its three-dimensional form. Hence, the need for a reduction to a lower level dimensional form using the following transformation variables [4, 26]:

ΞS,I,t=Smt,ω,ωS=kIT.(6)

Hence, (5) becomes:

mt+12σ2ω22mω21T+rωmω=0,mT,ω=φω.(7)

Obviously, (7) is now reformed to be in two dimensional whose solution will be used to obtain the Asian option price via the link in (6).

3 The Analysis of the Laplace Adomian decomposition method

In this sequel, we consider w.r.t. LADM [37, 38, 39, 40, 41, 42, 43, 44,] the following differential equation (partial or ordinary) of the form:

gζx,t=hx,t(8)

where g signifies a first order differential operator in t, which may be nonlinear, thereby including linear and nonlinear terms. Hence, (8) is decomposed as:

Ltζx,t+Rζx,t+Nζx,tgζx,t=hx,t(9)

where Lt=t,Ris a linear differential operator, N represents the nonlinear differential operator equivalent to an analytical nonlinear term, while h (x, t) is the associated source term. Hence, (9) becomes:

Ltζx,t=hx,tRζx,t+Nζx,t.(10)

We proceed by introducing the Laplace operator to differentiate the solution technique from the classical ADM. This is done as follows via the definitions:

Definition 1: Let f (t) be defined on t ∈ [0,∞), then the Laplace transform of f (t) is F (s) defined as:

Fs=L^ft=0ftestdt.(11)

Definition 2: For a continuous function f (t) such that F (s) = L̂{f (t)}, f (t) called the inverse Laplace transform (ILT) is defined as:

L^1Fs=ft.(12)

Definition 3: For an nth order differential equation, the associated Laplace transform is:

L^fnt=snL^fti=0n1sn1ifi0,L^tnft=1nFns,(13)

where (n) denotes the nth derivative with respect to t and with respect to s associated with f (n) (t) and F (n) (s) respectively.

The Laplace Transform (LT) is incorporated in the ADM [37, 38, 39, 40, 41, 42, 43, 44,] by taking the LT of both sides of (10) as follow:

L^Ltζx,t=L^hx,tRζx,t+Nζx,t.(14)

By using the derivative properties as noted in (13), therefore (14) becomes:

sζx,sζx,0=L^hx,tRζx,t+Nζx,t.(15)

It thus implies that:

ζx,s=1sζx,0+L^hx,tRζx,tNζx,t.(16)

Hence, for non-homogeneous cases (NHC) and homogeneous cases (HC), we have:

ζx,s=1sζx,0+L^hx,t1sL^Rζx,t+Nζx,t,NHC,1sζx,01sL^Rζx,t+Nζx,t,HC.(17)

Applying the ILT −1 (·)to (17) gives:

ζx,t=ζx,0+L^11sL^hx,tL^11sL^Rζx,t+Nζx,t,ζx,0L^11sL^Rζx,t+Nζx,t,(18)

Next, the LADM proposes representing the solution as an infinite series given as:

ζx,t=n=0ζnx,t(19)

with ζn (x, t) to be computed recursively. Also, ^the nonlinear term (x, t) is defined as:

Nζx,t=n=0Anζ0,ζ1,ζ2,,ζn(20)

and An referred to as Adomian polynomials is given as:

An=1n!nλnNi=0nλiζiλ=0.(21)

Therefore, substituting (19) and (20) in (18) gives:

n=0ζnx,t=ζx,0+L^11sL^hx,tQ.(22)

where:

Q=L^11sL^Rn=0ζnx,t+n=0An.

From (22), the solution ξ (x, t) is therefore determined via the recursive relation:

ζ0=ζx,0+L^11sL^hx,tζn+1=L^11sL^Rζn+NAn,n0(23)

while ξ (x, t) is finalized as:

ζx,t=limjn=0jζnx,t.(24)

4 Illustrative examples and applications

Here, the consideration of the analytical solution is made based on the proposed LADM with two illustrative cases.

Case I: Consider (5) via (6-7) in an operator form as follows:

mt=1T+rωmω12σ2ω2mωωm0,ω=φω,(25)

where the subscripts denote partial derivatives w.r.t. the subscripted ^variables.

Hence, ^taking the Laplace transform of (25) gives:

mω,s=1smω,0+L^1T+rωmω12σ2ω2mωω.(26)

Thus, applying the ILT L^1to (26) gives:

mω,t=mω,0+L^11sL^1T+rωmω12σ2ω2mωω.(27)

Therefore, with the initial condition, and the infinite series solution form: mω,t=n=0mnω,tthe following recursive relation is obtained via the proposed LADM:

m0=mω,0mn+1=L^11sL^1T+rωmn12σ2ω2mn,n0.(28)

Note: the prime notations in (27) denote derivatives w.r.t. ω.

Therefore, the recursive relation in (28) yields:

m0=mω,0m1=L^11sL^1T+rωm012σ2ω2m0,m2=L^11sL^1T+rωm112σ2ω2m1,m3=L^11sL^1T+rωm212σ2ω2m2,m4=L^11sL^1T+rωm312σ2ω2m3,m5=L^11sL^1T+rωm412σ2ω2m4,mk=L^11sL^1T+rωmk112σ2ω2mk1.(29)

Thus, we obtain the following by subjecting (29) to the initial condition:

m0=mω,0=1rT1erTωerT.(30)m1=1T+rωterT,m2=12!1T+rωrt2erT,m3=13!1T+rωr2t3erT,m4=14!1T+rωr3t4erT,m5=15!1T+rωr4t5erT,,mk=1k!1T+rωrk1tkerT,k1

Hence,

mω,t=j=0mj=1rT1erTωerT1T+rωt+t2r2!+erT=1rT1erTωerT1r1T+rωrt+rt22!+erT=1rT1erTωerT1ri=1rtii!1T+rωerT=1rT1erTωerT1r1+i=0rtii!1T+rωerT=1rT1erTωerT1r1+ert1T+rωerT=1rT1erTtωerTt,ω0.(31)

But from (6), ΞS,I,t=Smt,ω,ωS=kIT.

Therefore,

ΞS,I,t=SrT1erTtkITerTt.(32)

Equation (32) is the analytical solution of (5) corresponding to the continuous arithmetic Asian option pricing model.

Case II: Suppose (25) via (29) is considered based on a different initial condition:

m0=mω,0=ωr1rTSerT.(33)

Then, by the same approach, we have:

m1=T1+rωrtSerT,m2=12!T1+rωrt2SerT,m3=13!T1+rωrt3SerT,m4=14!T1+rωrt4SerT,m5=15!T1+rωrt5SerT,m6=16!T1+rωrt6SerT,,mp=1p!T1+rωrtpSerT,p1.

So,

mω,t=j=0mj=ωr1rT+1+ert1T+rωSerT.(34)...ΞS,I,t=ωr1rT+1+ert1T+rωS2erT.(35)

4.1 The Greeks of the Asian Option Model

Here, the Greeks (G1-G4) of the continuous AOPM are briefly introduced as their mathematical expressions are given. This is considered for a certain option value, ΞS,I,t=Ξ.

G1: The Theta-Greek of an option measures the rate of change of the option price w.r.t. the passage of time. Mathematically, the Delta is obtained by differentiating once the option value w.r.t the time variable say: Ξt=Ξt.

G2: The Delta-Greek of an option measures the rate of change of the option price w.r.t. the underlying asset price. It is the sensitivity of the option to the price of the asset. Mathematically, the Delta is obtained by differentiating once the option value w.r.t the spatial variable say: ΞS=Ξs.

G3: The Gamma-Greek of an option defines the rate of change of the Delta-Greek w.r.t. the spatial variable. Mathematically, the Gamma is obtained by differentiating the ΞSS=2Ξs2. Delta-Greek w.r.t the spatial variable say:

G4: The Speed-Greek of an option defines the rate of change of the Gamma-Greek w.r.t. the spatial variable. Mathematically, the Speed is obtained by differentiating the Gamma-Greek w.r.t the spatial variable say:ΞSSS=3Ξs3.

5 Conclusions

This paper considered the Greek parameters: Theta, Delta, Speed, and Gamma associated with a continuous arithmetic Asian option pricing model. The derivation is based on the analytical solution of the continuous arithmetic Asian option model obtained via a proposed semi-analytical method: Laplace-Adomian decomposition method (LADM). To the best of the Authors’ knowledge, the LADM is applied, for the first time, to the continuous arithmetic Asian option model for analytical solution. The solutions are provided in explicit form with few iterations, and less computational work is involved with high level of accuracy being maintained. For conformity, references are made to the analytical solutions obtained by Rogers & Shi (J. of Applied Probability 32: 1995, 1077-1088) [3], and Elshegmani & Ahmad (ScienceAsia, 39S: 2013, 67–69) [4]. The proposed method is highly recommended for analytical solution of other forms of Asian option pricing models such as the geometric put and call options, even in their time-fractional forms. The basic Greeks parameters obtained will be of great help to financial practitioners and traders in terms of hedging and portfolio management. Future research can include the application of the modified ADM, and the restarted ADM for speed and accuracy comparison.

Acknowledgement

The authors: SOE and GOA express sincere thanks to Covenant University for the provision of good working environment. Also, the constructive comments of the anonymous referees are highly appreciated.

References

  • [1]

    Esekon J.E., The Black-Scholes formula and the Greek parameters for a nonlinear Black-Scholes equation, Int. J. Pure Appl. Math., 2012, 76(2), 167-171 Google Scholar

  • [2]

    Kiptum P.J., Esekon J.E., Oduor O.M., Greek parameters of nonlinear Black-Scholes equation, Int. J. Math. Soft Comput., 2015, 5(2), 69-74 Google Scholar

  • [3]

    Rogers L.C.G., Shi Z., The value of an Asian option, J. Appl. Probability, 1995, 32, 1077-1088 Google Scholar

  • [4]

    Elshegmania Z.A., Ahmad R.R., Solving an Asian option PDE via the Laplace transform, Science Asia, 2013, 39S, 67–69 Google Scholar

  • [5]

    Falloon W., Turner D., The evolution of a market, In: Managing energy price risk, 1999, Risk Books, London, UK Google Scholar

  • [6]

    Edeki S.O., Owoloko E.A., Ugbebor O.O., The modified Black- Scholes model via constant elasticity of variance for stock options valuation, 2015 Progress in Applied Mathematics in Science and Engineering (PIAMSE), AIP Conference Proceedings, 2016, 1705, 4940289 Google Scholar

  • [7]

    Shokrollahi F., The evaluation of geometric Asian power options under time changed mixed fractional Brownian motion, J. Comp. Appl. Math., 2018, 344, 716- 724 Google Scholar

  • [8]

    Kemna A.G.Z., Vorst A.C.F., A pricing method for options based on average asset values, J. Banking and Finance, 1990, 14, 113-129 Google Scholar

  • [9]

    Edeki S.O., Ugbebor O.O., Owoloko E.A., On a dividend-paying stock options pricing model (SOPM) using constant elasticity of variance stochastic dynamics, Int. J. Pure Appl. Math., 2016, 106 (4), 1029-1036 Google Scholar

  • [10]

    Corsaro S., Kyriakou I., Marazzina D., Marino, Z., A general framework for pricing Asian options under stochastic volatility on parallel architectures, Europ. J. Operat. Res., 2019, 272 (3), 1082-1095 CrossrefGoogle Scholar

  • [11]

    Barucci E., Polidoro S., Vespri V., Some results on partial differential equations and Asian options, Math. Models Methods Appl. Sci. 2001, 11 (3), 475-497 CrossrefGoogle Scholar

  • [12]

    Geman H., Yor M., Bessel process, Asian options and perpetuities, Math. Finance 1993, 3 (4), 349-75 Google Scholar

  • [13]

    Vecer J., A new PDE approach for pricing arithmetic average Asian option, J. Comput. Finance, 2001, 4, 105-113 Google Scholar

  • [14]

    Zhang J., Theory of continuously-sampled Asian option pricing, (working paper), 2000, City University of Hong Kong Google Scholar

  • [15]

    Chen K., Lyuu Y., Accurate pricing formula for Asian options, J. Appl. Math. Comp., 2007, 188, 1711-1724 Google Scholar

  • [16]

    Elshegmani Z.A., Ahmad R.R., Zakaria R.H., New pricing formula for arithmetic Asian options using PDE Approach, Appl. Math. Sci, 2011, 5, 77, 3801–3809 Google Scholar

  • [17]

    Kumar A., Waikos A., Chakrabarty S.P., Pricing of average strike Asian call option using numerical PDE methods, 2011, arXiv:1106.1999v1 [q-fin.CP] Google Scholar

  • [18]

    Zhang B., Yu Y.,Wang W., Numerical algorithm for delta of Asian option, Sci. World J., 2015, 692847 

  • [19]

    Fadugba S.E., ThevMellin transformvmethod as an alternative analytic solution for the valuation of geometric Asian option, Appl. Comput. Math., 2014, 3, (6-1), 1-7 Google Scholar

  • [20]

    Elshegmani Z.A., Ahmad R.R., Analytical Solution for an arithmetic Asian option using mellin transforms, Int. J. Math. Analysis, 2011, 5 (26), 1259-1265 Google Scholar

  • [21]

    Edeki S.O., Ugbebor O.O., Owoloko E.A., Analytical solution of the time-fractional order Black-Scholes model for stock option valuation on no dividend yield basis, IAENG Int. J. Appl. Math., 2017, 47 (4), 407-416 Google Scholar

  • [22]

    Edeki S.O., Ugbebor O.O., Owoloko E.A. , He’s polynomials for analytical solutions of the Black-Scholes Pricing Model for Stock Option Valuation, Lecture Notes Eng. Comp. Sci., 2016, 2224, 632-634 Google Scholar

  • [23]

    Oghonyon J.G., Omoregbe N.A., Bishop S.A., Implementing an order six implicit block multistep method for third order ODEs using variable step size approach, Global J. Pure Appl. Math., 2016, 12 (2), 1635-1646 Google Scholar

  • [24]

    Akinlabi G.O., Edeki S.O., The solution of initial-value wavelike models via perturbation iteration transform method, Lecture Notes Eng. Comp. Sci., 2017, 2228, 1015-1018 Google Scholar

  • [25]

    Biazar J., Goldoust F., The Adomian decomposition method for the Black-Scholes equation, 3rd Int. Conf. on Appl. Math. Pharmaceutical Sci., Singapore, 2013, 321-323 Google Scholar

  • [26]

    Elshegmani Z.A., Ahmed R.R., Analytical solution for an arithmetic Asian option using mellin transforms, Int. J. Math. Analysis, 2011, 5 (26), 1259-1265 Google Scholar

  • [27]

    Akinlabi G.O., Edeki S.O., On approximate and closed-form solution method for initial-value wave-like models, International, J. Pure Appl. Math., 2016, 107 (2), 449-456 Google Scholar

  • [28]

    Edeki S.O., Akinlabi G.O., Odo C.E., Fractional complex transform for the solution of time-fractional advection-diffusion model, Int. J. Circuits, Systems and Signal Processing, 2017, 11, 425-432 Google Scholar

  • [29]

    Agarana M.C., Ede A.N., Application of differential transform method to vibration analysis of damped railway bridge on pasternak foundation under moving train, Lecture Notes Eng. Comp. Sci., 2016, 2224, 1177-1179 Google Scholar

  • [30]

    Edeki S.O., Akinlabi G.O., Nyamoradi N., Local fractional operator for analytical solutions of the K(2, 2)-focusing branch equations of time-fractional order, Int. J. Appl. Comput. Math, 2018, 4, 66 Google Scholar

  • [31]

    Guo Y., Globally robust stability analysis for stochastic Cohen-Grossberg neural networks with impulse control and time- varying delays, Ukrainian Math. J., 2018, 69 (8), 1220-1233 Google Scholar

  • [32]

    Moussi A., Lidouh A., Nqi F.Z., Estimators of sensitivities of an Asian option: numerical analysis, Int. J. Math. Analysis, 2014, 8 (17), 813-827 Google Scholar

  • [33]

    Guo Y., Mean square exponential stability of stochastic delay cellular neural networks, Electron. J. Qual. Theory Differ. Equ., 2013, 34, 1-10 Google Scholar

  • [34]

    Edeki S.O., Adeosun M.E., Owoloko E.A., Akinlabi G.O., Adinya I., Parameter estimation of local volatility in currency option valuation, Int. Rev. Model. Simul., 2016, 9(2), 130-133 Google Scholar

  • [35]

    Guo Y., Mean square global asymptotic stability of stochastic recurrent neural networks with distributed delays, Appl. Math. Comp., 2009, 215, 791-795 Google Scholar

  • [36]

    Aimi A., Guardasoni C., Collocation boundary element method for the pricing of Geometric Asian options, Engineering Analysis with Boundary Elements, 2018, 92, 90-100 Google Scholar

  • [37]

    Khuri S.A., A Laplace decomposition algorithm applied to a class of nonlinear differential equations, J. Appl. Math, 2001, 1 (4), 141-155 Google Scholar

  • [38]

    Fadaei J., Application of Laplace-Adomian decomposition method on linear and nonlinear system of PDEs, Appl. Math. Sci., 2011, 5 (27), 1307-1315 Google Scholar

  • [39]

    Pue-on P., Laplace Adomian Decomposition Method for Solving Newell-Whitehead-Segel Equation, Appl. Math. Sci., 2013, 7 (132), 6593-6600 Google Scholar

  • [40]

    Yousef H.M.,Md Ismail A.I.B.M., Application of the Laplace Adomian decomposition method for solution system of delay differential equations with initial value problem, AIP Conf. Proc., 2018, 1974, 020038 CrossrefGoogle Scholar

  • [41]

    Haq F. Shah , Ur Rahman G., Shahzad M., Numerical solution of fractional order smoking model via Laplace Adomian decomposition method, Alexandria Eng. J., 2018, 57 (2), 1061-1069 CrossrefWeb of ScienceGoogle Scholar

  • [42]

    Mishra V., Rani D., Newton-Raphson based modified Laplace Adomian decomposition method for solving quadratic Riccati differential equations, MATEC Web of Conferences, 2016, 57, 05001 CrossrefGoogle Scholar

  • [43]

    Singh R., Saha J., Kumar J., Adomian decomposition method for solving fragmentation and aggregation population balance equations, J. Appl. Math. Comput., 2015, 48, (1–2), 265-292 Google Scholar

  • [44]

    Matinfar M., Ghanbari M., The application of the modified variational iteration method on the generalized Fisher’s equation, J. Appl. Math. Comput., 2009, 31, (1-2), 165-175 Google Scholar

About the article

Received: 2018-05-05

Accepted: 2018-09-19

Published Online: 2018-12-14


Conflict of InterestConflict of Interests: The authors declare no conflict of interest regarding this paper.


Citation Information: Open Physics, Volume 16, Issue 1, Pages 780–785, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0097.

Export Citation

© 2018 S. O. Edeki et al., published by De Gruyter.. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Citing Articles

Here you can find all Crossref-listed publications in which this article is cited. If you would like to receive automatic email messages as soon as this article is cited in other publications, simply activate the “Citation Alert” on the top of this page.

[1]
Rajarama Mohan Jena, Snehashish Chakraverty, and Dumitru Baleanu
Mathematics, 2019, Volume 7, Number 8, Page 722
[2]
Rajarama Mohan Jena, Snehashish Chakraverty, and Dumitru Baleanu
Mathematics, 2019, Volume 7, Number 8, Page 689
[3]
Rajarama Mohan Jena and S. Chakraverty
SN Applied Sciences, 2019, Volume 1, Number 3

Comments (0)

Please log in or register to comment.
Log in