#### 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 of where 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.

#### 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, , with another symbol, which is an operational matrix, . 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 (where 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:subject to the initial condition .

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 Section 2. Section 2.3 presents analytical expression for Prabhakar integral and derivative for . 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:

*Definition 2. *The three-parameter generalization of Mittag-Leffler function is given by

*Definition 3. *For , Prabhakar fractional integral is defined bywhere , , and .

Normally, the definition in (4) can be written aswhereUsing (6), we have the following Prabhakar fractional derivative and regularized Prabhakar derivative.

*Definition 4. *For , , and Prabhakar fractional derivative (in Riemann–Liouville sense) is defined bywhere , , and .

*Definition 5. *For , , the regularized Prabhakar derivative (in Caputo sense) is defined bywhere , , and .

##### 2.2. Shifted Legendre Polynomials

The analytical form of the shifted Legendre polynomials of degree is given by

A function , which is square integrable in , can be expressed in terms of shifted Legendre polynomials aswhere the coefficient is given by

(11) can be written like that due to the orthogonality condition of shifted Legendre polynomials

In this work, only the first term of shifted Legendre polynomials is considered. So, we havewhere

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

In this section, we find the analytical expression for Prabhakar integration as well as Prabhakar derivative for where is positive integer.

Theorem 1. *The Prabhakar integration for where is positive integer can be defined as follows:*

*Proof. *For sake of simplicity, we let (that in equation (4)); by using Definition 3 and equations (5) and (6),By using integration by parts, we obtainSubstituting (17) into (16) and after some algebra manipulation, we obtainThe expression in Theorem 1 is equivalent to that in Lemma 4 [20]. If , Theorem 1, Definition 3, or Lemma 4 in [20] gives the same result which is .

Theorem 2. *The Prabhakar derivative of order for where is positive integer can be defined as follows:where is not equal to integer positive and .*

*Proof. *Taking in Definition 4, for sake of simplicity, we let (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:By using integration by parts, we obtainSubstituting (21) into (20) and after some algebra manipulation, we obtainwhere . (22) can be further reduced as follows:In (22), since the denominator consists of , we must have . For , Theorem 2 for is applicable for both Prabhakar derivative and regularized Prabhakar derivative. For regularized Prabhakar derivative [10], we have,The result obtained here can be verified via the series representation of fractional calculus operators involving generalized Mittag-Leffler functions (see Theorem 2.1 in 19).

#### 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. *Suppose is the shifted Legendre Polynomials vector**Let , . Then,where is operational matrix of fractional derivative of order in Prabhakar sense and is defined as follows:where is given byand 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: can be calculated using Theorem 2.Let ; by using truncated shifted Legendre polynomials, we have . Substituting this in (30), we obtainwhere is given in (28). Rewriting (31) in vector form, we havewhere . For , we haveHence, 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 for first few terms of shifted Legendre polynomials. The corresponding operational matrix when is given byThe comparison between the exact solution for Prabhakar derivative for and 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 .

#### 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, say . Using (26), we have . Also, we can have and .

*Step 2. *From , we obtain a system of algebraic equation for the unknown variables . Solving the system to obtain the values for , the solution of fractional differential equations in Prabhakar sense is given by .

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., ,If , using the operational matrix via shifted Legendre polynomials and approximate via shifted Legendre polynomials, we obtainHowever, in practical, is finite, say term of shifted Legendre polynomials had been used. There will be a small error, .Let is the approximation solution of (1) obtained by the operational matrix method, if are sufficiently small, then the absolute errors can be estimated by . Hence, the optimal value (i.e. ) 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.where . The exact solution is given by .

By using the proposed method, we obtain the following approximation as shown in Table 1. Here, we compare it with the iteration method introduced recently in [17] with . From the numerical result, by using very small terms, i.e., , we were able to obtain a good result compared to the iteration method in [17].

*Example 2. *Consider a simple fractional differential equation defined in Prabhakar derivative.where . The exact solution is given by .

By using the proposed method, we obtain the following approximation as shown in Table 2. Here, we compare it with the iteration method introduced recently in [17] with smaller . From the numerical result, by using very small terms, i.e., , 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.

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

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of the paper.