Novel Approaches for Solving Fuzzy Fractional Partial Differential Equations

: In this paper, we present a comparison of several important methods to solve fuzzy partial differential equations (PDEs). These methods include the fuzzy reduced differential transform method (RDTM), fuzzy Adomian decomposition method (ADM), fuzzy Homotopy perturbation method (HPM), and fuzzy Homotopy analysis method (HAM). A distinguishing practical feature of these techniques is administered without the need to use discretion or restricted assumptions. Moreover, we investigate the fuzzy ( n + 1 ) -dimensional fractional RDTM to obtain the solutions of fuzzy fractional PDEs. The much more distinctive element of this method is that it requires no predetermined assumptions, and reduces the computational effort. We apply the suggested techniques to a set of initial valued problems and get approximate numerical solutions for linear and nonlinear time-fractional PDEs. It is demonstrated that the fuzzy ( n + 1 ) -dimensional fractional RDTM is both accurate and simple to use. The methods are based on gH-differentiability and fuzzy fractional derivatives. Some illustrative numerical examples are given to demonstrate the effectiveness of our proposed methods. The results show that the methods are powerful mathematical tools for solving fuzzy partial differential equations.


Introduction
One of the most important areas of study in the fuzzy analysis is the differential and integral theory of fuzzy valued function, which is grounded in the idea of fuzzy number space. In particular, the fuzzy differential and integral equations, that are extensively used in engineering technology and social science, have piqued the interest of scholars from a variety of disciplines. The study of fuzzy differential equations is mostly based on the following three approaches; the first is based on the H-derivative and the generalized derivative of Bede. The second is considered under Zadeh's extension principle. The third is predicated on differential inclusion theory and fuzzy differential equations theory. These three explanations are different from one another.
In this work, we consider the H-derivative and the generalized derivative of Bede. We summarize the contributions and novelty as follows: • We present the comparison for a fuzzy (n + 1)-dimensional RDTM, ADM, VIM [1], and fuzzy HPM [2] demonstrates that even though the results of these approaches when implemented to the fuzzy wave-like and heat-like equations are the same. But, the fuzzy (n + 1)-dimensional RDTM, like fuzzy HPM, does not require specific algorithms and complex calculations such as fuzzy ADM or construction of correction functionals using general Lagranges multipliers in the fuzzy variational iteration method. In particular, the fuzzy RDTM and HPM are simple to apply and represent two successful techniques to obtain the solution of fuzzy PDEs.

•
We investigated the comparison of fuzzy (n + 1)-dimensional RDTM, ADM, HPM, and fuzzy HAM to obtain the solutions of fuzzy wave-like, heat-like and Zakharov-Kuznetsov equations. Although the results of these methods are the same when applied to problems. Moreover, the fuzzy (n + 1)-dimensional RDTM, HPM, and HAM don't require complex techniques and computations as fuzzy ADM. The results recall that the fuzzy RDTM, HPM, and HAM are easy to use for solving fuzzy partial differential equations.

•
We propose the solutions of fuzzy fractional wave-like, heat-like, and Zakharov-Kuznetsov equations using (n + 1)-dimensional fuzzy fractional RDTM. The method is flexible and can solve problems without calculating complicated Adomian polynomials or making unrealistic assumptions about nonlinear behavior. The provided technique is thus an influential way of solving fuzzy fractional PDEs and fractional order problems in physics, engineering, and other areas.
Fuzzy analysis and fuzzy differential equations have been proposed to deal with uncertainty due to incomplete information that appears in several mathematical or computer models of certain deterministic real-world phenomena. This theory has developed a large number of applications in which fuzzy fractional differential equations and fractional differential equations have emerged as important topics. Stefanini and Bede [3] proposed the generalized Hukuhara differentiation of interval-valued functions and interval differential equations. Also, Bede and Stefanini [4] introduced the generalized differentiation of fuzzy-valued functions. Gomes and Barros [5] discussed the generalized difference and the generalized differentiability. Hong et al. [6] presented an exhaustive review of various modern fractional calculus applications.
The fuzzy partial differential equations (FPDEs) have attracted great interest because of their practical applications in many fields such as physics, social science, and other areas of science and engineering. The FPDEs have been studied by many authors using different methods. Keshavarza et al. [19] presented the fuzzy solution to the mathematical model of a cancer tumor under Caputo-generalized Hukuhara partial differentiability by using fuzzy integral transforms. Keshavarz and Allahviranloo [20] studied the fuzzy fundamental triangular solution of the fractional diffusion equation under Caputo generalized Hukuhara partial differentiability by using the fuzzy Laplace transform and the fuzzy Fourier transform. Furthermore; see [1,[21][22][23][24]. The authors [25,26] presented the various transport/diffusion problem and an overview of the corresponding numerical solution approaches.
The differential transform method (DTM) was originally discussed by Zhou [27] in 1986, this technique adopts an analytic solution in polynomial form, which is different from the traditional higher-order Taylor formula technique. After that, many researchers have proposed this method to solve many problems [19,24,[28][29][30]. To overcome the demerits of complex computation of DTM, the RDTM was introduced by Keskin et al. [31,32] the method is based on reputable semi-analytical technique and can be applied to find approximate solutions of PDEs, also there are several significant implementations employing RDTM; see [32][33][34][35][36][37][38][39][40][41].
The Adomian decomposition method (ADM) is a well-known and effective approach for solving any type of problem. It is efficient not just for linear but also for nonlinear issues. This technique is famous for fast convergence and achieving the desired appropriate precision in just a few iterations. Several authors have already contributed their works via this technique; for example, see [1,24,[42][43][44].
He [45][46][47] is considered as the pioneer of HPM by combining HAM [48,49] and the perturbation method [50]. This method has been used to solve a wide range of problems with forwarding. Kashkari et al. [51] studied dissipative nonplanar solitons in an electronegative complex plasma by using the HPM. The HPM is used to solve both linear and nonlinear higher-order boundary value problems numerically by Kanth and Aruna [52]. This method was used by Biazar et al. [53] to solve nonlinear systems of integro differential equations. Osman et al. [24] compared the fuzzy HPM and other techniques applied to solving the fuzzy (1 + n)-dimensional Burgers equation. Xu [54] proposed a perturbational approach to construct analytical approximations based on the double-parameter transformation perturbation expansion method. Ahmad et al. [55] studied the nonlinear fractional order KdV and Burger equation with exponential-Decay Kernel using HPM.
The HAM [56,57] was introduced by Liao in 1992. HAM was further developed and improved by Liao in various subjects [58][59][60]. Several researchers have applied the HAM for solving differential equations. Saratha et al. [61] studied the notion of a fractional generalized integral transform under a modified Riemann-Liouville derivative with the Mittag-Leffler function as a kernel. Li et al. [62] presented the time-delay feedback control of a cantilever beam with concentrated mass based on the HAM. Naika et al. [48] studied the estimating an approximate analytical solution of the HIV viral dynamic model via HAM.
This paper is structured as follows. In Section 2, we recall some basic definitions. In Section 3, we applied the fuzzy (n + 1)-dimensional RDTM, ADM, HPM, and fuzzy HAM to obtain the solutions of fuzzy partial differential equations. In Section 4, we present the solution of fuzzy fractional partial differential equations via fuzzy (n + 1)-dimensional fractional RDTM. Finally, a conclusion is given in Section 5.

Preliminaries
In this paper, we will denote the set of fuzzy numbers by E 1 , that are, normal, fuzzy convex, upper semi-continuous and compactly supported fuzzy sets defined over the real line. [63,64]). Let u, v ∈ E 1 and k ∈ R, the addition and scalar multiplication are defined as The triangular fuzzy number defined as a fuzzy set in E 1 , determined by u = (a, b, c) ∈ R and a ≤ b ≤ c such that u σ = a + (b − a)σ and u σ = c − (c − b)σ are the endpoints of σ-level sets for all σ ∈ [0, 1]. A support of fuzzy number u is given as where cl is the closure of set {ϑ ∈ R|u(ϑ) > 0}.
The Hausdorff distance D : between fuzzy numbers is defined as in [65] where d H is the Hausdorff metric.
The metric space E 1 , D is complete, locally compact and the following properties from [65] for metric D are valid , as long as u v, and w e ∀ u, v, w, e ∈ E 1 , where is the H-difference, it means that w v = u if and only if u ⊕ v = w.
Definition 1 ([4,66]). The gH-difference between two fuzzy numbers u, v ∈ E 1 is defined as In terms of σ−levels, and if the H-difference exists, then u v = u gH v; the conditions for the existence of e = u gH v ∈ E 1 are It is easy to show that (i) and (ii) are both valid if and only if e is a crisp number.

Proposition 1 ([67]
). Let u, v ∈ E 1 are two fuzzy numbers. Then whenever the statement on the right exists, especially, u gH u = u u = 0. • If u gH v exists in sense (i), then v gH u exists in sense (ii) and vice versa.

Definition 3 ([3]).
We say that a point ϑ 0 ∈ (a, b) is a switching point for the differentiability of a function f if in any neighborhood V of ϑ 0 there exist points ϑ 1 < ϑ 0 < ϑ 2 such that • type I at ϑ 1 (4) holds while (5) does not hold and at ϑ 2 (5) holds and (4) does not hold, or • type II at ϑ 1 (5) holds while (4) does not hold and at ϑ 2 (4) holds and (5) does not hold.

Definition 4 ([63]
). Let f : [a, b] → E 1 and f gH (ϑ) be gH-differentiable at ϑ 0 ∈ (a, b) and there is no switching point on (a, b), with f (ϑ; σ) and f (ϑ; σ) are both differentiable at ϑ 0 . Then -differentiable whenever the type of gH-differentiability f (ϑ) and f gH (ϑ) is the same: Definition 5 ([68]). Let us suppose a function f : [a, b] → E 1 be fuzzy Riemann integrable in I ∈ R F if for any ε > 0 there exists δ > 0 such that for any division where ∑ * p denotes the fuzzy summation and I indicates b a f (ϑ)dx.

Fuzzy Partial Differential Equations
In this section, we present the solution of fuzzy partial differential equations. We considered the following fuzzy (n + 1)-dimensional reduced differential transform.

Fuzzy (n + 1)-Dimensional Reduced Differential Transform
We propose the fuzzy (n + 1)-dimensional reduced differential transform for solving fuzzy partial q-differential equations, the theory of (n + 1)-dimensional RDTM with uncertainty represented by using fuzzy concepts is explained as follows.
Thus, if ϑ ς (t) be (i)-differentiable, then ϑ ς (t) can be expressed as: and if ϑ ς (t) be (ii)-differentiable, then ϑ ς (t) can be expressed as: The mentioned equations are known as the inverse transformation of X(l; σ), which can be defined as , l is odd, , l is odd, and X (l; σ) = P(l) ∂ l (ϑ ς (t; σ)) ∂t l t=0 , l is even, , l is even, when ϑ ς (t) is (ii)-differentiable, then, the function ϑ ς (t) can be expressed as: when ϑ ς (t) are (i)-differentiable, and if ϑ ς (t) be (ii)-differentiable, we obtain where P(l) > 0, P(l) denoted the weighting factor. In this work P(l) = C l l! is applied, where C is the time horizon on interest. Consequently, if ϑ ς (t) be (i)-differentiable, then and if ϑ ς (t) be (ii)-differentiable, then and Unitizing the fuzzy (n + 1)-dimensional RDTM, the fuzzy PDEs in the particular domain is transformed into an algebraic equation in the domain K, and ϑ ς (t) is provided as the finite-term Taylor series plus a reminder as: when ϑ ς (t) is (i)-differentiable and In this section, we present the solution of fuzzy PDEs at the equally spaced grid points [t 0 , t 1 , ..., t n ] where t ς = a + ςl * for each (ς = 0, 1, 2, ...n), and l * = b−a n . That is, the domain of interest are proved to n is sub-domain, and the fuzzy approximation functions in each sub-domain are ϑ ς (t; σ) for ς = 0, 1, 2, ..., n − 1, respectively.

Proposition 2.
Let u(X , t) and v(X , t) are fuzzy-valued functions and their fuzzy (n + 1)dimensional reduced differential transformations denoted by U l (X ) and V l (X ), respectively. Then provided the generalized Hukuhara difference (gH-difference) exists.
Proof. By using definition (9), the proof is obvious.

Proposition 3.
Let us consider the fuzzy-valued function w ∈ E 1 and f (X , t) = ∂w(X ,t) ∂t , then we can obtain F l (X ) = (l+1)! l! W l+1 (X ), l ≥ 1 where F l (X ) and W l (X ) are the fuzzy (n + 1)dimensional reduced differential transformations of fuzzy-valued functions f and w, respectively.

Lemma 1.
Suppose w ∈ E 1 and f (X , t) = ∂w(X ,t) ∂ϑ ς , then we can obtain F l (X ) = ∂W l (X ) ∂ϑ ς , l ≥ 1 where F l (X ) and W l (X ) are the fuzzy (n + 1)-dimensional reduced differential transformations of fuzzy-valued functions f and w, respectively.
Proof. Using definition (9), we can obtain the following equation for σ ∈ [0, 1] Similarly, in view of definition (9) the fuzzy RDTM function can be written as: We achieve the result by differentiating the right side of the preceding equality with consideration to ϑ ς , Hence, the proof is completed by achieving our desired result. , l ≥ n where F l (X ) and W l (X ) are the fuzzy (n + 1)-dimensional reduced differential transformations of fuzzy-valued functions f and w, respectively.
Proof. Using definition (9), we obtain for 0 ≤ σ ≤ 1 From the calculus, one can obtain Consequently, the fuzzy (n + 1)-dimensional RDTM of fuzzy-valued function the proof is completed.
Theorem 1. Let W l (X ) and G l (X ) are the (n + 1)-dimensional fuzzy RDTM of w(X , t) is a positive real-valued function and g(X , t) is a fuzzy-valued function. Also let us suppose that if f (X , t) = w(X , t)g(X , t), then Proof. Using definition (9), we get In general, we obtain and from the definition of (n + 1)-dimensional RDTM, we obtain This completes our required proof.

Proof. From definition
This means .., ϑ ı n n , the required proof is completed.

Theorem 2.
Let us consider the real-valued function w ∈ R and f (X , t) = w(X ) · g(X , t), then F l (X ) = w(X ) · G l (X ), where F l (X ) and G l (X ) are (n + 1)-dimensional RDTM of real-valued functions f and g, respectively.

Applications
In this section, we propose some examples in [1,2] to illustrate the applicability of the alternative approach of fuzzy (n + 1)-dimensional RDTM to obtain the solutions of fuzzy heat-like and wave-like equations with variable coefficients.

Example 1.
We consider the following fuzzy (2 + 1)-dimensional heat-like equation [1,2] ∂w ∂t with the initial condition where n = 1, 2, 3, ... Applying the fuzzy reduced differential transform to (36), we get Similarly, applying fuzzy reduced differential transformation on the initial condition (37) to achieve Putting Equations (40) into (38), we obtain thus, we can achieve the solution of w(ϑ, θ, t; σ) as follows: Example 2. Consider the following fuzzy (3 subject to the initial condition where Applying the fuzzy (n + 1)-dimensional reduced differential transform on (41) to get Using the initial condition (42), we obtain Substituting (46) into (43), we obtain the series solution as we can obtain the exact solution as:

Using the fuzzy RDTM for
From initial conditions (48), we obtain Substituting (52) into (49), we get the series solution as: We can obtain the exact solution as:

Example 4. Consider the following fuzzy
with the initial conditions where n = 1, 2, 3, ... Applying (53), we get Taking Equation (54) yields Using (58) into (55), we get the series solution as: We can find the exact solution as: When this method is compared to other methods in [1,2], it shows that when these methods are used to solve fuzzy heat-like and wave-like equations, they all lead to the same proposed solution. In addition, fuzzy (n + 1)-dimensional RDTM like HPM doesn't always involve specific algorithms and complex calculations like fuzzy ADM or the development of correction functionals utilizing general Lagranges multipliers in the fuzzy VIM. So, the fuzzy (n + 1)-dimensional RDTM is a better way to solve fuzzy partial differential equations and is also simple and easy to use.

Fuzzy Zakharov-Kuznetsov Equations
In this part, we present the nonlinear fuzzy Zakharov-Kuznetsov equations as follows: (59) subject to the initial condition where Υ 1 , Υ 2 , Υ 3 are the arbitrary constants and m, n, l are integrals.

Fuzzy Adomian Decomposition Method
Consider the following formal nonlinear fuzzy differential equation as: where L is a linear differential operator, R denotes the linear operator's remainder, and N w denotes the nonlinear terms. We can obtain (61) using the inverse operator L −1 on both sides Firstly, (59) can be represented as and Suppose that L −1 and an integral operator defined by Using the integral operator L −1 on both sides of (59), we get The fuzzy decomposition method assumes a series solution forw(ϑ, θ, t; σ) given by an infinite sum of components as: where w 0 , w 1 , w 2 , ... are obtained sequentially. The nonlinear terms are decomposed into three infinite polynomial series where A l , B l , and C l are Adomian polynomials, which can be used to determine all types of nonlinearities using fuzzy Adomian's techniques. The analytical formulae for Adomian polynomials are: For the nonlinear operators (69), we provide the first few Adomian polynomials and and Using (70) into (68), we obtain We use the recursive relation to identifying the components w l (ϑ, θ, t), l ≥ 0, as We assume that all of the components w ς (ϑ, θ, t; σ) are calculated in light of (75) into (71) and Convergence analysis of the fuzzy ADM can be found in (Theorem 3.3, [24]).

The Fuzzy Homotopy Perturbation Method
We consider the following general nonlinear fuzzy differential equation under the boundary condition where B denotes the boundary operator, ∂Φ denotes the boundary of the domain Φ,w(℘) denotes the analytical function, and A is a general differential operator. The fuzzy operator A can be broken into fuzzy linear L and nonlinear N parts. Hence, Equation (76) can be rewritten as: We generate a homotopy using the fuzzy homotopy technique: where ∈ [0, 1] denote the embedding parameter, and forw 0 (℘) denote the initial approximation to (76) which satisfies the boundary conditions. Clearly, from (80), we obtain and and the changing process of from zero to unity is just thatṽ(℘, ; σ) from w 0 (℘; σ) to w(℘; σ). Applying the Homotopy parameter as an extending parameter to obtain As a result of → 1, the approximate solution of (76) is given as Convergence analysis of the fuzzy HPM can be found in (Theorem 3.4, [24]).

The Fuzzy Homotopy Analysis Method
We consider the following fuzzy differential equation as: where0∈ E 1 , N is a nonlinear operator, ℘ and t were independent variables, and w(℘, t; σ) denote the unknown fuzzy-valued function, respectively. For simplicity, we disregard all boundary or initial conditions, that can be handled in a similar manner. Constructions for the so-called zero-order deformation equation are made possible through the generalization of the classical homotopy technique.
for σ ∈ [0, 1] denotes the fuzzy number, ∈ [0, 1] denotes the embedding parameter, h = 0 denotes a non-zero auxiliary parameter, H(℘, t) = 0 denotes the non-zero auxiliary function, and L denotes the auxiliary linear operator with the follows: w 0 (℘, t; σ) shows an initial guess forw(℘, t; σ), andw(℘, t, ; σ) = [w(℘, t, ; σ), w(℘, t, ; σ)] presents an unknown fuzzy-valued function. It the important to note that HAM provides a large amount of flexibility in choosing auxiliary items. Clearly, this is accurate for = 0 and = 1, when the quantity increases from 0 to 1, the solutionφ(℘, t, ), changes from the initial guesses, Taylor series can be extended with respect to : where If such auxiliary linear operator, the initial approximation, the auxiliary parameterh, and the auxiliary fuzzy-valued function are all appropriately determined, and the series (96) and (97) converges at = 1. Then, we obtain the following result: According to (98) and (99), the governing equation can be deduced from the zero-order deformation Equations (90) and (91). Define the vector The m th order deformation equation is obtained by differentiating Equations (90) and (91) times with respect to parameter at = 0 where and

Applications
In this section, we present examples 5 and 6 to illustrate the discussed methods for effectiveness by solving Zakharov-Kuznetsov equations.

(119)
The solution in a series form as Similarly, the series solution of w(ϑ, θ, t; σ) on the Formula (116) can be determined as follows: Thus, we have obtained the exact solution w(ϑ, θ, t; σ) of (111) as

Case [B]. Fuzzy Homotopy perturbation method.
Applying the fuzzy HPM, we construct a homotopy as follows We consider the initial approximation that satisfies the initial condition Substituting (85) and (86), with (122), and equating the terms of identical powers of p is
Similarly, we can represent the series solution of w(ϑ, θ, t; σ) in Equation (148) as: Using Taylor series into (150) and (151), we obtained the closed form solution

(158)
Next, we can get the series solutions Similarly, the series solution of w(ϑ, θ, t; σ) on Formula (155) can be derived as follows: According to Taylor series into (158), we obtain
Substituting (85) and (86) into Equations (161) and (162) and equating the terms with identical powers of p, we have and
In Figure 1, we plotted 2D and 3D graphs of the ZK(2, 2, 2) equation. Figure 1a shows that for ϑ = 30, θ = 45 and ρ = 1 using n = 1 at t = 0.001 the ZK(2, 2, 2) equation is bounded and closed. Furthermore, the blue + sign shows increasing functions and red * presents decreasing functions on the σ-level set of w. To discuss the concept of the σ-level set, one can see Figure 2a, which shows that the σ-level set of ZK(2,2,2) equation is bounded and closed for ϑ = 30 and 0 < θ ≤ 2π. Similarly, in Figure 2, we can observe the same explanation of σ-level set closedness and boundedness for example 6.

Fuzzy Fractional Partial Differential Equations
In this section, we present the solution of fuzzy fractional partial differential equations via fuzzy (n + 1)-dimensional fractional RDTM.

Fuzzy Fractional Calculus
We regard to C F [a, b] as the space of all continuous fuzzy-valued functions on [a, b]. Also, we denote the space of all Lebesgue integrable fuzzy-valued functions on the bounded interval [a, b] ⊂ R by L F [a, b], refs. [71].
Additionally, using Taylor series for ϑ ς (t ς ; σ), the solution on the grid points t ς+1 can be obtained as: The Properties of Fuzzy (N + 1)-Dimensional Fractional Reduced Differential Transform We investigate some mathematical operations of fuzzy (n + 1)-dimensional fractional reduced differential transform.

Lemma 5.
Let us consider u(X , t) and v(X , t) are fuzzy-valued functions and their fuzzy (n + 1)dimensional fractional RDTM denoted by U αl (X ) and V αl (X ), respectively. Then proposed the generalized Hukuhara difference (gH-difference) exists.
Proof. According to Definition (13), the proof is obvious.
Proof. Using Definition (13), we obtain for 0 ≤ σ ≤ 1 Using definition of fuzzy fractional RDTM, we obtain the proof is completed.
Proof. From definition (13), we obtain for 0 ≤ σ ≤ 1 The (n + 1)-dimensional fuzzy fractional RDTM function is written as: Using differentiating the right side of the mentioned equality with respect to ϑ ς , we obtain the proof is completed. , l ≥ n where F αl (X ) and G αl (X ) are the fuzzy (n + 1)-dimensional fractional reduced differential transformations of fuzzy-valued functions f and g, respectively.
Note: Assuming η = nα, then the expression above can be represented as follows:
Theorem 5. Let us consider u ∈ R and f (X , t) = u(X )g(X , t), then F αl (X ) = u(X )G αl (X ), where F αl (X ) and G αl (X ) are the fuzzy (n + 1)-dimensional fractional reduced differential transformations of real-valued functions f and g, respectively.

Examples
We propose some examples to illustrate this method is a powerful mathematical tool for solving fuzzy fractional partial differential equations.
Using the properties of fuzzy (n + 1)-dimensional fractional RDTM, we have and W α(l+n) (ϑ, θ, φ; σ) = Γ(αl + 1) Γ(α(l + n) + 1) ϑ 2 + θ 2 + φ 2 δ(αl) Taking the initial conditions (210), we have and W 0 =0, for β = 1.5, i.e., n = 3, α = 1 2 , and l = 0, 1, 2, 3, ... (214) into (211), we have and For β = 2, i.e., n = 2, α = 1, and l = 0, 1, 2, 3, ... (214) into (211), we have thus, we can obtained the exact solution as: The results corresponding to example 7 are shown in Figure 3 at different values of β. But, if we compare it with others methods in [1,2] shows that although the result of these methods implemented the same at β = 2. But, unlike fuzzy ADM or the generation of correction functionals using general Lagranges multiplication in fuzzy VIM. The fuzzy (n + 1)-dimensional fractional RDTM does not call for additional algorithms and complicated calculations. Table 1 shows the error term between exact and approximate solutions of example 7 for σ between 0 and 1. We have also checked and verified the convergence for time t in this example, which shows that example 7 exhibit convergent solutions till time t = 709 and as the value of t exceeds 709, the solutions tend to infinity and show divergence. In Figure 3a, we have compared solutions of fuzzy wave-like equations based on integer as well as fractional order derivatives. It can be seen that red and blue colored are for exact solution using β = 2, while orange and purple colored dashed-dotted lines are for fractional order at β = 1.5.
For specific values of ϑ = 0.02, θ = 0.002, φ = 0.03 the solution of fuzzy fractional wave-like equations at β = 1.5 and 2 are same. Therefore, for a detailed study, we plot a three-dimensional Figure 3b in which we fix all the parameters except θ. Here, one can observe in detail that at the start there exists an error in the exact and approximate solution which reduces time and finally the approximate solution overlaps the exact solution.
Similar to previous examples, here we have also checked the convergence for time t, which shows that example 8 exhibit convergent solutions till time t = 355 and as the value of t exceeds 355, the solutions tend to infinity and show divergence.
In Figure 4, we plotted 2D and 3D graphs of the ZK(2, 2, 2) equation but with different initial condition. Figure 4a shows that for ϑ = 0.0001, θ = 0.05, φ = 0.6 and ρ = 1 using n = 1 at t = 0.07 the ZK(2, 2, 2) equation become bounded and closed. Furthermore, the pink colored sign shows increasing functions and blue colored presents decreasing functions on the σ-level set of w.
To discuss the concept of the σ-level set, one can see Figure 4b, which shows that the σ-level set of ZK(2,2,2) equation is bounded and closed for ϑ = 0.0001, 0 < θ < 1 and φ = 0.6.
Finally, the convergence for example 9 shows that their solutions are convergent till time t = 4254. Figure 5 also satisfies the condition of σ-level set in both (two and three dimensional) cases for example 9.

Conclusions
In this paper, we have successfully compared (n + 1)-dimensional fuzzy RDTM, ADM, HPM, and fuzzy HAM to obtain the solutions of fuzzy heat-like and wave-like equations, and fuzzy Zakharov-Kuznetsov equations. Furthermore, we investigated the fuzzy (n + 1)-dimensional fractional RDTM to apply the solution of fuzzy fractional heatlike and wave-like equations, and fuzzy Zakharov-Kuznetsov equations. The RDTM is applied in an uncomplicated approach, without discretization or limiting assumptions. Previous numerical studies demonstrated that the RDTM is occasionally more effective than other techniques. We demonstrated that the suggested methods are highly accurate and efficient by applying them to some of the initial value problems. Hence, we have obtained several new results to solve the above problems when these methods have been applied. Moreover, we observed that our methods are strong mathematical tools for solving PDEs and issues in physics, engineering, and other fields. In future, we are trying our best to present new techniques for solving fuzzy fractional diffusion equations, and the numerical technique for solving fuzzy fractional Cauchy reaction-diffusion equations as well.