In this manuscript, we study a nonlinear fractional-order predator-prey system while considering uncertainty in initial values. We derive the feasibility region and the boundness of the solution. The suggested model’s equilibrium points and the basic reproduction number are calculated. The stability of equilibrium points is presented. We use the metric fixed point theory to study the existence and uniqueness results concerning the solution of the model. We use the notion of UH-stability to show that the model is Ulam–Hyres type stable. To attain the approximate solution of the proposed model, we construct a method that uses the fuzzy Laplace transform in collaboration with the ADM (Adomian decomposition method). Finally, we simulate our theoretical results using MATLAB to show the dynamics of the considered model.

1. Introduction

A mathematical model is a mathematical representation of a process. The different phenomenon of the real world has been modelled by DEs or by a system of DEs [1]. Various techniques such Taylor series method, variational iteration method, modified variational iteration algorithms, and homotopy perturbation method have been discussed to solve different models [25]. Here, we consider the simplest model named as “Lotka–Volterra” or Predator-Prey, presented independently by “Lotka” in (1925) and “Volterra” in (1926) [6]. It refers to the relationship between predator and prey controls, as well as the population development of both species. The model describes the evolution of a biological Predator-Prey relationship through a system of two nonlinear DEs, generalized to a more complex and realistic phenomenon [69]. The general system for the dynamic phenomenon is as follows:where is the prey (rabbits), is the predator (foxes) population at a time , and are its functions, respectively. The growth rate and the density at time have a direct relationship. This is presumed that the populations live in a community, where the age factor is not considered. The prey is abundant in natural resources, and the only threat is the specialist predator, whose growth depends on prey catches. For the prey model, the rate of growth of prey is presumed to be constant, and the particular rate of growth is reduced by an amount equal to the population of the predator [6]. The proposed model of the classical order is governed by

From the century, the subject of fractional-order integrals and derivatives has a substantial effect on modelling and simulation due to its nonlocal nature [10, 11]. It is seen that models involving integrals and derivatives of fractional order are more accurate, rather than classical models because fractional operators have long memory and heredity properties [12]. Das and Gupta [13] studied system (2) under fractional differential equations with constant coefficients [13]. Ahmad et al. [14] analyzed the human liver model by using hybrid fractional operators [15]. Alderremy et al. [16] analyzed multi space-fractional Gardner equation in [16]. For more applications of different fractional operators, see [1720]. It should be noted that, in many branches of mathematics, DEs have been expanded to fill the study gape. Since we live an uncertain environment, therefore, we cannot model any phenomenon accurately. So, fuzzy operators are the best tools to handle the uncertain situation. Therefore, classical calculus is generalized to fractional and fuzzy Calculus, while DEs are generalized to fractional and fuzzy DEs [2123]. Fuzzy DEs can model a physical phenomenon more accurately by considering fuzziness or uncertainty in initial conditions. To solve such problems, they utilized different techniques [24, 25]. Recently, fuzzy DEs have been utilized to analyze various models that occur in biology and physics [15, 26, 27]. As in the classical model, the coefficients are constant; therefore, we extend the model given in (2) to fuzzy fractional-order operator with uncertain initial data. For order, and .where

Here, represents fractional derivative in the Caputo sense, and are fuzzy valued functions, and are fuzzy numbers.

The goal of this research is to look at the model under consideration by including fuzziness in the initial values. The model’s basic dynamical features are studied. Fixed point theory and nonlinear analysis are used to derive the existence and stability conclusions of model (3). LADM is used to find the solution to the presented model. To demonstrate the model’s dynamics, the obtained results are simulated for various fractional-order values.

The structure of the paper is as follows: Section 1 deals with the introduction and motivation part of the paper. Section 2 provides the basic concepts of fuzzy fractional calculus. Section 3 deals with stability and feasibility of the equilibrium points. Section 4 provides the existence and uniqueness theory of the proposed model. Section 5 discusses the Ulam–Hyres stability via nonlinear functional analysis. Section 6 gives the general procedure for the solution of the considered model via fuzzy Laplace transform. Section 7 demonstrates the model through 2D and 3D simulations. Section 8 includes conclusion of the manuscript.

2. Preliminaries

Definition 1. (see [28]). Consider a fuzzy set of real line, i.e., . Then, is called a fuzzy number if it fulfills the following features:(1)The closure of the set is compact(2) is upper semicontinuous on (3) is convex (4) is normal for some ;

Definition 2. (see [28]). If is a fuzzy number, then the parametric form is represented by such that and fulfills the properties which are given as follows:(1) is bounded, left continuous, and increasing function in and right continuous at 0(2) is bounded, right continuous, and decreasing in and right continuous at 0(3)Also, for crisp case, we have .

Definition 3. (see [25]). Let contain all fuzzy numbers. We take into consideration a mapping and take and as two fuzzy numbers. The Hausdorff distance between and is presented as follows:In , the metric satisfies the properties which is given as follows:(1) (2) (3) (4) is a complete metric space

Definition 4. (see [28]). Let be a fuzzy set valued function. Then, is called continuous, if for each and fixed value of , the following relation holds:

Definition 5. (see [28]). Let be a continuous fuzzy function on , then for , a fuzzy Riemann–Liouville fractional integral is given byAdditionally, if , where and are the spaces that contains fuzzy Lebesgue integrable functions and fuzzy continuous, respectively, thenwhere

Definition 6. (see [28]). Consider a fuzzy function with parametric form , then the Caputo fuzzy fractional-order derivative is presented bywhereso that the right side integration converges and is bracket function of .

Definition 7. (see [28]). The fuzzy Laplace transform (LT) for continuous and fuzzy Reimann-integrable function on is defined as follows:Using parametric form of , we have for :Hence,

Theorem 1 (see [21]). The LT of Caputo fractional derivative of is presented as follows:

3. Feasibility and Stability Analysis

In this part, the feasibility and stability of the points of equilibrium will be discussed. To discuss the model’s boundedness and feasibility, we have the following theorem.

Theorem 2. The solution of the proposed model is bound to the feasible region given by

Proof. By adding both equation of (2) and considering , one getsOn using simple integration, we haveThis implies that . Consequently, the result is obtained.

For stability, the equilibrium points for (3) must be found as

We have two points of equilibrium which are as follows: if and if .

Theorem 3. The reproduction number for (1) is .

Proof. To find the reproduction number, consider equation of (1) as ,such that , where and represent the linear and nonlinear terms, respectively. The next step is to determine the next generation matrix , wherethenThe leading eigen value is equal to at free equilibrium point in , which can be computed asSimilarly, we can compute the reproduction number for equation of (1). Hence, .

Theorem 4. If , then free equilibrium points of (1) are locally asymptotically stable.

Proof. Proof is obvious.

4. Existence Theory

Here, we will establish the existence theory for the model by using fixed point theory. Let us take into consideration the right hand sides of the given system as follows:

Thus, the given system (3) can be written as follows: for and ,where

On using the corresponding integral , we have

We denote Banach space as under the fuzzy norm as follows:

Then, we reachwhere and .

We take some assumption on nonlinear function as follows.

C-1 . For each such that

C-2 . There exist constants and such that

Theorem 5. Let assumption C-2 holds, then model (25) possesses at least one solution.

Proof. Let us consider a closed and convex fuzzy set . Let us define a mapping such thatFor any , we haveThe last relation shows that . From the last inequality, we have . It follows that is bounded operator. Next, the next step is to show that the operator is completely continuous. For this let such that , thenThe last inequality implies that

Theorem 6. Suppose that C-1 holds, then the considered model (25) possesses the unique solution if .

Proof. Let , thenThis implies that . Hence, fulfills the contraction condition. Thus, by “Banach fixed point theorem” system (25) possesses the unique solution.

5. Ulam–Hyres Stability

We take a small alteration such that depends only on the solution as follows:(i)(ii)

Lemma 1. The solution of the perturbed problemfulfills the relation which is given as follows:

Proof. The proof is straightaway.

Theorem 7. Assuming along with (38), the solution of equation (37) is UH-stable. Thus, the numerical results of the model being considered are UH-stable if .

Proof. Let be the solution of equation (37) and be the unique solution of equation (3), thenFrom (39), we getHence, we infer from (39) that the solution of (38) is UH-stable. Consequently, the model (3) is UH-stable.

6. Derivation of General Procedure for Solution

Here, we derive an algorithm for the solution of the proposed model by fuzzy LADM. Consider system (3) with initial conditions as follows:where

Then, we reach

The infinite series solution is presented as follows:

Thus, we reach

Comparison of terms on both sides gives

Taking inverse Laplace transform, we get

Therefore, we get

7. Results and Simulations

In this part, we illustrate the proposed method through simulations; consider the system above under some specific values of parameters as follows:where

After applying the proposed method, we get

Similarly, the next terms are computed as

We can obtain other terms of the series, i.e.,and so on. The unknown terms in the above equations are given as

Now, we simulate our obtained results to investigate the dynamics of the proposed model for the uncertain initial conditions. In Figures 1 and 2, we plot the fuzzy approximate results up to three terms for the given examples under two different uncertainty values as follows:

We displayed the approximate fuzzy solutions for the considered model for particular fuzzy initial conditions and at the given uncertainty against different fractional orders in Figures 1 and 2. The fuzzy fractional derivative produced global dynamics of the interaction of the two species under fuzzy concept that under taking some values for uncertainty . Increasing the uncertainty, the dynamics of interaction is also affected as shown in Figures 1 and 2. Such situation is natural and may be observed in daily life. Next, in Figure 3 and Figure 4, we attempt on 3D plots of fuzzy solutions for first three terms as follows: by surfaces plots given in Figures 3 and 4, we have graphed the solutions for different fractional order over the interval of uncertainty for the given problem 15.

8. Conclusion

In this research, we have used a fuzzy fractional derivative in Caputo sense to extend the fractional predator-prey model. We have discussed the equilibrium points and their stability. Using metric fixed point theory, we were able to derive the existence and uniqueness results for the nonlinear fuzzy fractional predator-prey model. We have shown that the proposed model is HU stable through nonlinear functional analysis. Besides, we have developed a general algorithm for obtaining an approximate solution to the proposed model using an efficient method (fuzzy LADM). If we replace in the obtained numerical results, then we recover the results obtained in the fractional-order model. Thus, our proposed model is the generalization of the fractional model of the predator-prey system. Finally, we simulated numerical results for different fractional order and 0.5 and 0.9 uncertainty values via MATLAB. From figures, we have found that fractional calculus can be combined with the fuzzy theory and thus glorify the global dynamics of species interaction.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding this research work.