# Numerical Solution via Operational Matrix for Solving Prabhakar Fractional Differential Equations

Authors: Farah Suraya Md Nasrudin; Chang Phang
Journal: Journal of Mathematics
 (2022)

---

## Abstract

In this work, we apply the operational matrix based on shifted Legendre polynomials for solving Prabhakar fractional differential equations. The Prabhakar derivative is defined in three-parameter Mittag-Leffler function. We achieve this by first deriving the analytical expression for Prabhakar derivative ofxp where p is positive integer, via integration. Hence, for the first time, the operational matrix method for Prabhakar derivative is derived by using the properties of shifted Legendre polynomials. Hence, we transform the Prabhakar fractional differential equations into a system of algebraic equations. By solving the system of algebraic equations, we were able to obtain the numerical solution of fractional differential equations defined in Prabhakar derivative. Only a few terms of shifted Legendre polynomials are needed for achieving the accurate solution.


---

## Body

## 1. Introduction

The operational matrix method is one of the powerful tools for solving fractional differential equations. This method uses the concept of replacing a symbol with another symbol, i.e., replacing symbol fractional derivative,Dα, with another symbol, which is an operational matrix, Pα. In [1], the authors had derived shifted Legendre operational matrix for solving fractional differential equations, defined in Caputo sense. Then, researchers started to apply the various types of polynomials to derive the operational matrix for solving various types of fractional calculus problems, including Genocchi operational matrix for fractional partial differential equations [2], Laguerre polynomials operational matrix for solving fractional differential equations with non-singular kernel [3], and Müntz–Legendre polynomial operational matrix for solving distributed order fractional differential equations [4].Recently, apart from the fractional differential equation defined in Caputo sense, this kind of operational matrix method had been extended to tackle another type of fractional derivative or operator, which includes the Caputo–Fabrizio operator [5] and Atangana–Baleanu derivative [6, 7]. In this research direction, the operational matrix method is either an operational matrix of derivative or operational matrix of integration based on certain polynomials. The operational matrix method is possible to apply to another type of fractional derivatives if there is an analytical expression for xp (where p is integer positive) in the sense of certain fractional derivatives or operators. Hence, we extend this operational matrix to tackle operator defined by one parameter Mittag–Leffler function, i.e. Antagana–Baleunu derivative [6] to the operator that defined by using three-parameter Mittag–Leffler function, so-called Prabhakar fractional integrals or derivative. In short, we aim to solve the following fractional differential equation defined in Prabhakar sense:(1)Dα,β,ωγyx=gx,yx,subject to the initial condition y0=a.On top of that, Prabhakar introduced a type of convolution-type integral operator, called Prabhakar integral in [8]. Recently, Prabhakar fractional integrals or derivative had received more and more attention by the researchers [9–11]. This kind of integral had been applied in anomalous dielectrics [12], viscoelasticity [13], kinetic equation [14], and diffusion equation [15]. Besides that, some new concepts or theories were derived to suit this Prabhakar operator, for example, in [16], Hyers–Ulam stability of fractional differential equations with Prabhakar derivatives was investigated, and stability analysis of fractional differential equations with Prabhakar derivative was studied in [17]. Furthermore, since this Prabhakar operator involves more parameters, existing numerical methods may not be applicable for solving the fractional differential equation defined in the Prabhakar sense. Hence, some works had been done including numerical approximation to Prabhakar fractional Sturm–Liouville problem [18], and the Prabhakar derivative was approximated by using series representations in [19]. However, the work for solving numerically the Prabhakar fractional differential equations is still less. Hence, we intend to derive a new operational matrix based on shifted Legendre polynomials to approximate Prabhakar derivative, hence solving the Prabhakar fractional differential equations by using a collocation scheme.This paper is organized as follows. We will briefly explain some preliminary concepts including the Prabhakar fractional integral and derivative in Section2. Section 2.3 presents analytical expression for Prabhakar integral and derivative for xp. Section 3 discusses the derivation of a new operational matrix based on shifted Legendre polynomials for Prabhakar fractional derivative. In Section 4, we explain the new scheme and the error analysis. Some examples for solving fractional differential equations defined in Prabhakar derivative using our proposed method via the new operational matrix will be presented in Section 5. Conclusion and some recommendations are highlighted in Section 6.


## 2. Preliminaries

### 2.1. Prabhakar Fractional Integral and Derivative

In this section, we will present some basic concepts related to Prabhakar fractional integral and derivative.Definition 1.
The one-parameter Mittag-Leffler function is defined as follows:(2)Eαz=∑k=0∞zkΓαk+1,Reα>0.Definition 2.
The three-parameter generalization of Mittag-Leffler function is given by(3)Eα,βγz=∑k=0∞Γγ+kΓγΓαk+βzkk!,Reα>0.Definition 3.
Forf∈L1a,b, Prabhakar fractional integral is defined by(4)Iα,β,ω,a+γfx=∫axx−τβ−1Eα,βγωx−ταfτdτ,where Reα>0, Reβ>0, and γ,α,β,ω∈ℂ.
Normally, the definition in (4) can be written as(5)Iα,β,ω,a+γfx=∫axeα,βγx−τ;ωfτdτ,where(6)eα,βγx;ω=xβ−1Eα,βγωxα.
Using (6), we have the following Prabhakar fractional derivative and regularized Prabhakar derivative.Definition 4.
For0<β<1, fx∈L1a,b, and Prabhakar fractional derivative (in Riemann–Liouville sense) is defined by(7)Dα,β,ω,a+γRfx=dmdxm∫axeα,m−β−γx−τ;ωfτdτ,where Reα>0, Reβ>0, and γ,α,β,ω∈ℂ.Definition 5.
Forfx∈ACma,b, 0≤a<x<b≤∞, the regularized Prabhakar derivative (in Caputo sense) is defined by(8)Dα,β,ω,a+γcfx=∫a+xeα,m−β−γx−τ;ωdmfτdτmdτ,where Reα>0, Reβ>0, and γ,α,β,ω∈ℂ.


### 2.2. Shifted Legendre Polynomials

The analytical form of the shifted Legendre polynomialsLjx of degree j is given by(9)Ljx=∑h=0j−1j+hj+h!j−h!h!2xh.A functionyx, which is square integrable in 0,1, can be expressed in terms of shifted Legendre polynomials as(10)yx=∑j=0∞cjLjx,where the coefficient cj is given by(11)cj=2j+1∫01yxLjxdx,j=0,1,2,….(11) can be written like that due to the orthogonality condition of shifted Legendre polynomials(12)∫01LixLjxdx=12i+1fori=j,0fori≠j.In this work, only the firstN+1 term of shifted Legendre polynomials is considered. So, we have(13)yx=∑j=0NcjPjx=CTLx,where(14)CT=c0,c1,c2,…,cN,Lx=L0x,L1x,L2x,…,LNxT.


### 2.3. Analytical Expression for Prabhakar Integral and Derivative forxp

In this section, we find the analytical expression for Prabhakar integration as well as Prabhakar derivative forfx=xp where p is positive integer.Theorem 1.
The Prabhakar integration forfx=xp where p is positive integer can be defined as follows:(15)Iα,β,ωγxp=∑k=0∞Γγ+kωkp!xαk+β+pΓγΓαk+β+p+1Γk+1.Proof.
For sake of simplicity, we leta=0 (that in equation (4)); by using Definition 3 and equations (5) and (6),(16)Iα,β,ωγxp=∫0xx−τβ−1∑k=0∞Γγ+kωkx−ταkΓγΓαk+βk!τpdτ=∑k=0∞Γγ+kωk∫0xx−τβ−1x−ταkτpdτΓγΓαk+βk!.
By using integration by parts, we obtain(17)∫0xx−ταk+β−1τpdτ=Γαk+βp!xαk+β+pΓαk+β+p+1.
Substituting (17) into (16) and after some algebra manipulation, we obtain(18)Iα,β,ωγxp=∑k=0∞Γγ+kωkΓαk+βp!xαk+β+p/Γαk+β+p+1ΓγΓαk+βk!=∑k=0∞Γγ+kωkp!xαk+β+pΓγΓαk+β+p+1Γk+1.
The expression in Theorem1 is equivalent to that in Lemma 4 [20]. If γ=α=1,β=2,ω=3,p=2, Theorem 1, Definition 3, or Lemma 4 in [20] gives the same result which is 2/81e3x−1/9x3−1/9x2−2/27x−2/81.Theorem 2.
The Prabhakar derivative of order0<γ,α,β<1 for fx=xp where p is positive integer can be defined as follows:(19)Dα,β,ωγxp=∑k=0∞Γ−γ+kωkΓp+1xαk−β+pΓ−γΓαk−β+p+1Γk+1,where γ is not equal to integer positive and α≠β.Proof.
Takingm=1 in Definition 4, for sake of simplicity, we let a=0 (that in equation (7)), and using equation (6) and applying the similar approach as in the proof of Theorem 1, we obtain the following results:(20)Dα,β,ωγxp=ddx∫0xeα,1−β−γx−τ;ωτpdτ=∫0x−x−τ−β−1β∑k=0∞Γ−γ+kωkx−ταkΓ−γΓαk−β+1k!τpdτ+∫0xx−τ−β∑k=0∞Γ−γ+kωkx−ταkkαx−τΓ−γΓαk−β+1k!τpdτ=−β∑k=0∞Γ−γ+kωk∫0xx−ταk−β−1τpdτΓ−γΓαk−β+1k!+α∑k=0∞Γ−γ+kωk∫0xx−ταk−β−1xpdτΓ−γΓαk−β+1Γk.
By using integration by parts, we obtain(21)∫0xx−ταk−β−1τpdτ=Γαk−βp!xαk−β+pΓαk−β+p+1.
Substituting (21) into (20) and after some algebra manipulation, we obtain(22)Dα,β,ωγxp=−β∑k=0∞Γ−γ+kωkΓp+1xαk−β+pΓ−γΓαk−β+p+1αk−βΓkk+α∑k=0∞Γ−γ+kωkΓp+1xαk−β+pΓ−γΓαk−β+p+1Γkαk−β,where α≠β. (22) can be further reduced as follows:(23)Dα,β,ωγxp=∑k=0∞Γ−γ+kωkΓp+1xαk−β+pΓ−γΓαk−β+p+1Γk+1.
In (22), since the denominator consists of αk−β, we must have α≠β. For m=1, Theorem 2 for xp is applicable for both Prabhakar derivative and regularized Prabhakar derivative. For regularized Prabhakar derivative [10], we have,(24)Dα,β,ωγCxp=0,p=0,1,…,m−1,m=β.
The result obtained here can be verified via the series representation of fractional calculus operators involving generalized Mittag-Leffler functions (see Theorem 2.1 in19).


## 2.1. Prabhakar Fractional Integral and Derivative

In this section, we will present some basic concepts related to Prabhakar fractional integral and derivative.Definition 1.
The one-parameter Mittag-Leffler function is defined as follows:(2)Eαz=∑k=0∞zkΓαk+1,Reα>0.Definition 2.
The three-parameter generalization of Mittag-Leffler function is given by(3)Eα,βγz=∑k=0∞Γγ+kΓγΓαk+βzkk!,Reα>0.Definition 3.
Forf∈L1a,b, Prabhakar fractional integral is defined by(4)Iα,β,ω,a+γfx=∫axx−τβ−1Eα,βγωx−ταfτdτ,where Reα>0, Reβ>0, and γ,α,β,ω∈ℂ.
Normally, the definition in (4) can be written as(5)Iα,β,ω,a+γfx=∫axeα,βγx−τ;ωfτdτ,where(6)eα,βγx;ω=xβ−1Eα,βγωxα.
Using (6), we have the following Prabhakar fractional derivative and regularized Prabhakar derivative.Definition 4.
For0<β<1, fx∈L1a,b, and Prabhakar fractional derivative (in Riemann–Liouville sense) is defined by(7)Dα,β,ω,a+γRfx=dmdxm∫axeα,m−β−γx−τ;ωfτdτ,where Reα>0, Reβ>0, and γ,α,β,ω∈ℂ.Definition 5.
Forfx∈ACma,b, 0≤a<x<b≤∞, the regularized Prabhakar derivative (in Caputo sense) is defined by(8)Dα,β,ω,a+γcfx=∫a+xeα,m−β−γx−τ;ωdmfτdτmdτ,where Reα>0, Reβ>0, and γ,α,β,ω∈ℂ.


## 2.2. Shifted Legendre Polynomials

The analytical form of the shifted Legendre polynomialsLjx of degree j is given by(9)Ljx=∑h=0j−1j+hj+h!j−h!h!2xh.A functionyx, which is square integrable in 0,1, can be expressed in terms of shifted Legendre polynomials as(10)yx=∑j=0∞cjLjx,where the coefficient cj is given by(11)cj=2j+1∫01yxLjxdx,j=0,1,2,….(11) can be written like that due to the orthogonality condition of shifted Legendre polynomials(12)∫01LixLjxdx=12i+1fori=j,0fori≠j.In this work, only the firstN+1 term of shifted Legendre polynomials is considered. So, we have(13)yx=∑j=0NcjPjx=CTLx,where(14)CT=c0,c1,c2,…,cN,Lx=L0x,L1x,L2x,…,LNxT.


## 2.3. Analytical Expression for Prabhakar Integral and Derivative forxp

In this section, we find the analytical expression for Prabhakar integration as well as Prabhakar derivative forfx=xp where p is positive integer.Theorem 1.
The Prabhakar integration forfx=xp where p is positive integer can be defined as follows:(15)Iα,β,ωγxp=∑k=0∞Γγ+kωkp!xαk+β+pΓγΓαk+β+p+1Γk+1.Proof.
For sake of simplicity, we leta=0 (that in equation (4)); by using Definition 3 and equations (5) and (6),(16)Iα,β,ωγxp=∫0xx−τβ−1∑k=0∞Γγ+kωkx−ταkΓγΓαk+βk!τpdτ=∑k=0∞Γγ+kωk∫0xx−τβ−1x−ταkτpdτΓγΓαk+βk!.
By using integration by parts, we obtain(17)∫0xx−ταk+β−1τpdτ=Γαk+βp!xαk+β+pΓαk+β+p+1.
Substituting (17) into (16) and after some algebra manipulation, we obtain(18)Iα,β,ωγxp=∑k=0∞Γγ+kωkΓαk+βp!xαk+β+p/Γαk+β+p+1ΓγΓαk+βk!=∑k=0∞Γγ+kωkp!xαk+β+pΓγΓαk+β+p+1Γk+1.
The expression in Theorem1 is equivalent to that in Lemma 4 [20]. If γ=α=1,β=2,ω=3,p=2, Theorem 1, Definition 3, or Lemma 4 in [20] gives the same result which is 2/81e3x−1/9x3−1/9x2−2/27x−2/81.Theorem 2.
The Prabhakar derivative of order0<γ,α,β<1 for fx=xp where p is positive integer can be defined as follows:(19)Dα,β,ωγxp=∑k=0∞Γ−γ+kωkΓp+1xαk−β+pΓ−γΓαk−β+p+1Γk+1,where γ is not equal to integer positive and α≠β.Proof.
Takingm=1 in Definition 4, for sake of simplicity, we let a=0 (that in equation (7)), and using equation (6) and applying the similar approach as in the proof of Theorem 1, we obtain the following results:(20)Dα,β,ωγxp=ddx∫0xeα,1−β−γx−τ;ωτpdτ=∫0x−x−τ−β−1β∑k=0∞Γ−γ+kωkx−ταkΓ−γΓαk−β+1k!τpdτ+∫0xx−τ−β∑k=0∞Γ−γ+kωkx−ταkkαx−τΓ−γΓαk−β+1k!τpdτ=−β∑k=0∞Γ−γ+kωk∫0xx−ταk−β−1τpdτΓ−γΓαk−β+1k!+α∑k=0∞Γ−γ+kωk∫0xx−ταk−β−1xpdτΓ−γΓαk−β+1Γk.
By using integration by parts, we obtain(21)∫0xx−ταk−β−1τpdτ=Γαk−βp!xαk−β+pΓαk−β+p+1.
Substituting (21) into (20) and after some algebra manipulation, we obtain(22)Dα,β,ωγxp=−β∑k=0∞Γ−γ+kωkΓp+1xαk−β+pΓ−γΓαk−β+p+1αk−βΓkk+α∑k=0∞Γ−γ+kωkΓp+1xαk−β+pΓ−γΓαk−β+p+1Γkαk−β,where α≠β. (22) can be further reduced as follows:(23)Dα,β,ωγxp=∑k=0∞Γ−γ+kωkΓp+1xαk−β+pΓ−γΓαk−β+p+1Γk+1.
In (22), since the denominator consists of αk−β, we must have α≠β. For m=1, Theorem 2 for xp is applicable for both Prabhakar derivative and regularized Prabhakar derivative. For regularized Prabhakar derivative [10], we have,(24)Dα,β,ωγCxp=0,p=0,1,…,m−1,m=β.
The result obtained here can be verified via the series representation of fractional calculus operators involving generalized Mittag-Leffler functions (see Theorem 2.1 in19).


## 3. Operational Matrix for Prabhakar Fractional Derivative

In this section, we will derive the new operational matrix based on shifted Legendre Polynomials for Prabhakar fractional derivative.Theorem 3.
SupposeLx is the shifted Legendre Polynomials vector(25)Lx=L0x,L1x,L2x,…,LNxT.
Let0<α<1, m=1. Then,(26)Dα,β,ωγyx=Pα,β,ωγLx,where Pα is N+1×N+1 operational matrix of fractional derivative of order α in Prabhakar sense and is defined as follows:(27)Pα,β,ωγ=00⋯0⋮⋮⋯⋮00⋯0∑h=ββθj,h,α,β,ω,0γ∑h=ββθj,h,α,β,ω,1γ⋯∑h=ββθj,h,α,β,ω,Nγ⋮⋮⋯⋮∑h=βjθj,h,α,β,ω,0γ∑h=βjθj,h,α,β,ω,1γ⋯∑h=βjθj,h,α,β,ω,Nγ⋮⋮⋯⋮∑h=βNθj,h,α,β,ω,0γ∑h=βNθj,h,α,β,ω,1γ⋯∑h=βNθj,h,α,β,ω,Nγ,where θj,h,α,β,ω,lγ is given by(28)θj,h,α,β,ω,lγ=∑h=αj−1j+hj+h!j−h!h!2∑k=0∞Γ−γ+kωkΓh+1Γ−γΓαk−β+h+1Γk+1cl,and cl can be obtained from inner product via (11) and ⌈⌉ is the ceiling function.Proof.
From (9), we can write the shifted Legendre polynomials in analytical form and its fractional derivative in Prabhakar sense is given as in the following equation:(29)Dα,β,ωγLx=∑h=0j−1j+hj+h!j−h!h!2Dα,β,ωγxh.
Dα,β,ωγxh can be calculated using Theorem 2.(30)Dα,β,ωγLx=∑h=0j−1j+hj+h!j−h!h!2∑k=0∞Γ−γ+kωkΓh+1xαk−β+hΓ−γΓαk−β+h+1Γk+1.
Letfx=xαk−β+h; by using truncated shifted Legendre polynomials, we have fx=∑l=0NclLlx. Substituting this in (30), we obtain(31)Dα,β,ωγLx=∑l=0N∑h=βj−1j+hj+h!j−h!h!2∑k=0∞Γ−γ+kωkΓh+1Γ−γΓαk−β+h+1Γk+1clLlx=∑l=0N∑h=βjθj,h,α,β,ω,lγLlx,where θj,h,α,β,ω,lγ is given in (28). Rewriting (31) in vector form, we have(32)Dα,β,ωγLx=∑h=βjθj,h,α,β,ω,0γ∑h=βjθj,h,α,β,ω,1γ…∑h=βjθj,h,α,β,ω,NγLx,where j=β…N. For j=0,1,…,β−1, we have(33)Dα,β,ωγLx=0,0,…0Lx,j=0,1,…,β−1.
Hence, by combining (32) and (33), the Legendre operational matrix for Prabhakar derivative is proved as in (27).
Since this is the first time that operational matrix for Prabhakar derivative is derived, to show that the operational matrix of Prabhakar derivative is correct, we find the exact solution of Prabhakar derivative withγ=1/2,α=1/2,β=1/4,ω=1/4 for first few terms of shifted Legendre polynomials. The corresponding operational matrix when N=2 is given by(34)0001.1416468810.8989013133−0.1256075578−0.89890131370.80360123681.256894671.
The comparison between the exact solution for Prabhakar derivative forL1x and L2x with the approximation using operational matrix as in (27) is shown in Figures 1 and 2. The accuracy of the approximation can be increased by using bigger N.Figure 1
Comparison of the exact solution and approximation forL1x.Figure 2
Comparison of the exact solution and approximation forL2x.


## 4. Proposed Scheme and Error Analysis

Basically, in order to solve fractional differential equations in Prabhakar sense as in equation (1) using the operational matrix method, the following procedure can be applied.Step 1.
Write each terms of fractional differential equations in terms of shifted Legendre polynomials, sayDα,β,ωγyx+yx=gx. Using (26), we have Dα,β,ωγyx=Pα,β,ωγLx. Also, we can have yx=∑k=0NckLx and gx=∑k=0NgkLx.Step 2.
FromPα,β,ωγLx+∑k=0NckLx=∑k=0NgkLx, we obtain a system of algebraic equation for the unknown variables ck. Solving the system to obtain the values for ck, the solution of fractional differential equations in Prabhakar sense is given by yx=∑k=0NckLx.
For the error estimation for the numerical scheme, let us consider the residual correction procedure which can be used to estimate the absolute error. From equation (1), i.e., Dα,β,ωγyx=gx,yx,(35)Dα,β,ωγyx−gx,yx=0.
IfN⟶∞, using the operational matrix via shifted Legendre polynomials and approximate gx,yx via shifted Legendre polynomials, we obtain(36)Pα,β,ωγLx−gx,yNx≈0,N⟶∞.
However, in practical,N is finite, say m term of shifted Legendre polynomials had been used. There will be a small error, em.(37)Dα,β,ωγyx−Pα,β,ωγLx2+gx,y∞x−gx,yNx2=em,N=m.
Letem∗ is the approximation solution of (1) obtained by the operational matrix method, if em−em∗ϵ are sufficiently small, then the absolute errors em can be estimated by em∗. Hence, the optimal value m (i.e. N) can be obtained.


## 5. Numerical Examples

In this section, some examples are presented to illustrate the applicability and accuracy of this new operational matrix for the Prabhakar derivative. All the computations are done by using the software Maple.Example 1.
Consider a simple fractional differential equation defined in Prabhakar derivative.(38)Dα,β,ωγyx=−32x7/4315Γ3/4π72xΓ3422F214,34;32,134;x−152F2−14,14;12,114;xπ,where γ=1/2,α=1/2,β=1/4,ω=1. The exact solution is given by yx=x2.
By using the proposed method, we obtain the following approximation as shown in Table1. Here, we compare it with the iteration method introduced recently in [17] with h=0.02. From the numerical result, by using very small terms, i.e., N=2, we were able to obtain a good result compared to the iteration method in [17].Table 1
Comparison of the absolute errors obtained by proposed method with the iteration method used in [17] for Example 1.
xExact solutionAbs. Error (iteration method)Abs. Error (proposed method)0.10.012.67306E-033.28100E-030.20.044.29289E-035.83288E-030.30.095.99314E-037.65566E-030.40.167.75504E-038.74932E-030.50.259.56803E-039.11388E-030.60.361.14254E-028.74932E-030.70.491.33229E-027.65566E-030.80.641.52571E-025.83288E-030.90.811.72261E-023.28100E-03Example 2.
Consider a simple fractional differential equation defined in Prabhakar derivative.(39)Dα,β,ωγyx+yx=2x2+322Γ3/4x5/45π3F3−16,16,12;13,23,94;x64−72x19/12133Γ7/123F316,12,56;23,43,3112;x64−9x23/12506Γ11/123F312,56,76;43,53,3512;x64,where γ=1/2,α=1/3,β=3/4,ω=1/4. The exact solution is given by yx=2x2.
By using the proposed method, we obtain the following approximation as shown in Table2. Here, we compare it with the iteration method introduced recently in [17] with smaller h=0.01. From the numerical result, by using very small terms, i.e., N=2, we were able to obtain a good result compared to the iteration method in [17].
As can be seen from both examples, the absolute error for iteration method introduced in [17] will increase when the iteration step is increasing. This will not occur for the operational matrix method.Table 2
Comparison of the absolute errors obtained by proposed method with the iteration method used in [17] for Example 2.
xExact solutionAbs. Error (iteration method)Abs. Error (proposed method)0.10.021.94868E − 032.68860E − 030.20.083.93195E − 034.51364E − 030.30.185.92093E − 035.47514E − 030.40.327.91259E − 035.57308E − 030.50.509.90584E − 034.80748E − 030.60.721.19001E − 023.17832E − 030.70.981.38952E − 026.85620E − 040.81.281.58908E − 022.67063E − 030.91.621.78870E − 026.89044E − 03


## 6. Conclusion

In this paper, we had successfully derived a new operational matrix based on shifted Legendre polynomials for solving fractional differential equations in Prabhakar sense. Only a few terms of shifted Legendre polynomials are needed to obtain a good approximation. Using the same process, the operational matrix can be derived using other types of polynomials. For future work, we hope to extend the existing fractional calculus problem such as in [21, 22] to Prabhakar fractional derivative.



---

*Source: 7220433-2022-05-25.xml*