Abstract

In this paper, the efficient combined method based on the homotopy perturbation Sadik transform method  (HPSTM) is applied to solve the physical and functional equations containing the Caputo–Prabhakar fractional derivative. The mathematical model of this equation of order with is presented as follows: where for and , equations are changed into the equal width and modified equal width equations, respectively. The analytical method which we have used for solving this equation is based on a combination of the homotopy perturbation method and Sadik transform. The convergence and error analysis are discussed in this article. Plots of the analytical results with three examples are presented to show the applicability of this numerical method. Comparison between the obtained absolute errors by the suggested method and other methods is demonstrated.

1. Introduction

The integrals and derivatives of fractional order in fractional calculus play an important role in many branches such as mathematics, engineering, and physics [16]. The differential systems of fractional order are regarded as an extension of the differential systems of integer order and have important application in expressing nonlinear phenomena in many scopes such as applied sciences [7, 8], engineering [9, 10], and physics [11, 12]. One of the most important nonlinear differential equations used in many fields of mathematics, engineering, and physics is the equal width (EW) differential equation of fractional order where the equation is a partial differential equation that describes physical behaviors such as water transfer in soils [13], crystal growth [14], and shallow water waves [15]. The main reason for choosing these types of equations is due to their application in mathematics and physics, plasma waves, fluid mechanics, solid-state physics, and chemical physics [16, 17]. In this paper, we state the time-fractional comprehensive EW equation containing Caputo–Prabhakar fractional derivative as follows [1820]:where , and is an integer. For equation (1), we consider two cases as follows:(1)If are considered, then equation (1) becomes the nonlinear equal width equation of fractional order . This type of equation is one of the most important equations that examine the nonlinear behaviors of physical phenomena such as hydromagnetic waves [21], optics [22], and biological systems [23].(2)If are considered, then equation (1) becomes the nonlinear modified equal width equation (MEW) of fractional order .

Moreover, in equation (1), function is the probability density function, is the spatial coordinate, and is the temporal coordinate. To find the solutions of these types of equations, the homotopy perturbation Sadik transform method is used, which is a generalization of the homotopy perturbation and Sumudu transform methods. The method presented in this article is similar to the methods discussed in [2431]. In equation (1), is the Caputo–Prabhakar fractional derivative of order which is defined bywhere is the Prabhakar fractional integral and is defined by [32]where and that is the Sobolev space. Also, is the three-parameter Mittag-Leffler function which is defined by [32]and is the Pochhammer symbol which is given by [33]

The main reason for choosing the three-parameter Mittag-Leffler function in this paper is related to its application in many models such as disordered materials and heterogeneous models [34], Havriliak–Negami models [35, 36], viscoelasticity models [37], stochastic models [38], probability models [39], spherical stellar models [40], Poisson models [41], and fractional models or integral models [4245]. Many research studies have been conducted in numerical fields to obtain numerical solutions of fractional differential equations such as the Adomian decomposition method [46], the q-homotopy analysis transform method (q-HATM) [47], the Laplace transform  method [48], the method based on the implementation of an iterative perturbation method [49], the numerical method based on the Petrov–Galerkin method [50], the Petrov–Galerkin finite element scheme [51], the local discontinuous Galerkin method [52], and other methods [24, 27, 30, 53]. Recently, a new integral transform named the Sadik transform was introduced by Sadikali Latif Shaikh in 2018 (see [54, 55]). The Sadik transform is nothing but an unification of the Laplace transform, Sumudu transform, Elzaki transform, and all those integral transforms whose kernels are of exponential type or similar to the kernel of the Laplace transform. The Sadik transform is a very powerful integral transform. In this paper, we apply the HPSTM to solve the nonlinear time-fractional EW equation, MEW equation, and VMEW equation. The key intention of the present work is to extend the utilization of the HPSTM to derive analytical and approximate solutions of the time-fractional EW equation, time-fractional MEW equation, and time-fractional VMEW equation.

This paper is organized as follows. Some definitions and mathematical preliminaries of the fractional calculus are presented in Section 2. In Section 3, we introduce an analytical method based on the HPSTM to the nonlinear nonhomogeneous partial differential equation. In Section 4, to show the validity of the suggested method in Section 3, three examples are simulated.

2. Important Notations

In this section, we express theorems and lemmas which are used in the next parts.where and .

Lemma 1. The Laplace transformation of function (4) is given by [33, 56]

Lemma 2 (see [3]). The following relation holds:where and .

Lemma 3 (see [57]). Suppose is a Sadik transform of and is a Sadik transform of . Then, Sadik transform of is given bywhere is a convolution. Also, and .

Theorem 1 (see [57]). Suppose and , are continuous. Then,

Also, the Sadik transform of integration of is defined by

Theorem 2. The Sadik transform of the Caputo–Prabhakar fractional derivative of order for is obtained bywhere and .

Proof. By applying equations (2), (6), (8), and (9), we obtainThe proof is proven.

3. The Proposed Scheme

This section focuses on an analytical method based on the homotopy perturbation and Sadik transform method for finding the solutions of the following fractional differential equations:where is a linear differential operator, is a nonlinear differential operator, and is an indicated source term. Using equation (11) and taking the Sadik transform on both sides of equation (13), we get

Now applying Sadik inverse transform on equation (14), we obtain

To get the solution of equation (15), we can state the solutions of equation (15) by the following infinite series:where , are known. Also, the nonlinear term can be displayed as the following infinite series:where are the Adomian polynomials which were introduced in [46]. Now, substituting equations (16) and (17) into (15), we getwhere . Therefore, we put the coefficients on powers of on both sides of relation (18) equally, which results in

Then, the analytical solution of equation (13) can be obtained as

Theorem 3. Let be the exact solution of equation (13) and be the approximate solution of equation (13) so that . Then, the series defined by (20) converges.

Proof. From , we obtainThen,Therefore,Because , as . The proof is proven.

Theorem 4. Let be the exact solution of equation (1) and be the best approximate solution of equation (1) so that . Then, the following relation for the error holds:

Proof. Since is the exact solution of equation (1) and is the best approximate solution of equation (1), then we haveUsing Lemma 2, equation (3), and , for , we getand by using the Lipschitz condition with constant for derivative of the first order with respect to the variable , we haveBy using the Lipschitz condition with constant for derivative of the first order with respect to the variable and theorem assumptions, from equation (25), we obtainand according to Theorem 2 introduced in [58], we havewhere , , that , given [58].

4. Numerical Results

In this section, we show the approximate solution obtained by the proposed method which is named the homotopy perturbation Sadik transform method with three examples to express its performance. In this section, we consider the absolute error as follows:where is the exact solution of equation (1) and is the numerical solution of equation (1).

4.1. Example 1

We consider the following fractional EW equation with :when , and the exact solution is . By applying the proposed method given in Section 3, we solve equation (31); then,where is the nonlinear . The polynomials are calculated as follows:

Comparing coefficients of same powers of in equation (32), we obtain

Then, the approximate solution of is obtained as

The comparison between the exact and the approximate solutions and the absolute error are expressed in Figure 1 for different values of when . Also, in Figure 2, the exact and the numerical solutions at are given. In Table 1, a comparison between the absolute error in [2] and the absolute error for HPSTM with different values of is given.

4.2. Example 2

Consider the following MEW equation with :for , and the exact solution is . Applying the HPSTM on equation (36), we obtainwhere the polynomials are represented as the nonlinear terms which are denoted as

The polynomials are calculated as follows:

Comparing coefficients of same powers of in equation (37), we obtain

Finally, the approximate solution of is given by

The comparison between exact solution and the approximate solution and absolute error are expressed in Figure 3 for different values of when . Figure 4 shows the comparison of exact solution with homotopy perturbation Sadik transform method solution at . In Table 2, a comparison between the absolute error in [2] and the absolute error for HPSTM with different values of is given.

In the example 2, we study the solution of a variant of fractional modified equal width equation (VMEW) which is obtained by the proposed method.

4.3. Example 3

Consider the following VMEW equation:and the exact solution is when . Applying the homotopy perturbation and Sadik transform method on equation (42), we havewhere are obtained asand are calculated as follows:

Comparing coefficients of same powers of in equation (43), we obtain

So, the approximate solution of is calculated as

Comparison between the exact and the approximate solutions and the absolute error are shown in Figure 5 for different values of when . Also, in Figure 6, comparison of the exact solution and the approximate solution at is illustrated. In Table 3, a comparison between the absolute error in [2] and the absolute error for homotopy perturbation Sadik transform method with different values of is given.

5. Conclusion

This paper describes a numerical method based on the homotopy perturbation Sadik transform method for solving the equations such as EW, MEW, and VMEW. The fractional derivatives are used in this paper in the Caputo–Prabhakar sense. This paper demonstrates that the HPSTM is sufficient, easy, and suitable for various other nonlinear models. It can be seen that as increases, it leads to variation in for waves in plasma at specific value of . All tables and figures show that the proposed method produces numerical solutions with more accuracy. Comparison between the obtained absolute errors by the suggested method and method presented in [2] is demonstrated. Discussion on the convergence and error analysis of the proposed method is presented. Three numerical examples are given to show the applicability of the suggested method.

Data Availability

No data were used to support this study.

Conflicts of Interest

The author declares that there are no conflicts of interest.