Abstract
In this research, we aim to propose a non-integer order mathematical model that illustrates the detrimental impact of illegal logging on forest biomass. The model is based on the idea that soil nutrients play a vital role in the growth of forestry biomass. However, the growth of forestry biomass is hindered by the harvesting of trees for industrial expansion. To develop the model, we use a fractional order approach, which allows us to explore the impact of conservation efforts on the extension of industries. If the amount of forestry biomass that is illegally harvested exceeds a certain threshold, the system becomes unstable. To investigate the proposed model's solution, we use fixed-point theory, which is highly effective in determining the existence and uniqueness of the solution. We also use Ulam-Hyers stability to analyse the stability of the proposed model. We have also employed new versions of fractional numerical techniques using power law and exponential decay. Our proposed model provides more details through numerical calculations with fractional orders different from integer ones. Moreover, we have addressed all significant solution behaviours in the fractal-fractional system and parameter values. The graphical representation offers a better understanding of the forestry biomass model, allowing us to identify the direction and stability of solutions. From various perspectives, we analysed the features of solutions to this fractal-fractional model for stopping illegal logging.
1. Introduction
An area of mathematics known as fractional calculus focuses on implementations for arbitrary order integrals and derivatives. Conversely, there was limited conversation surrounding fractional calculus, especially regarding its practical uses. Nonetheless, fractional calculus has emerged as one of the most sought-after academic fields in recent years, thanks to its theoretical progress and potential applications in engineering and scientific disciplines [Citation1]. A system's fundamental dynamics depend on the past behaviour of real-world mathematical modelling problems [Citation2]. As a result, the system's future depends on its current state and history of state and dependents. The integer order differential operator cannot capture the behaviour of such systems in hereditary terms. In contrast, the fractional order derivative, which is nonlocal, is a valuable tool for modelling nonlocal occurrences in both space and time. Non-integer order operators describe every step of the process study in the proposed model. Because of this, such operators are good tools for describing the systems' memory and hereditary characteristics. The nonlocal property of non-integer order operators contributed to the success of fractional calculus. Mathematically modelling real-world systems in terms of fractional-order derivatives generates several fractional-order differential equations frequently intractable by analytical means, leading to high computational expenses to render such solutions realistic [Citation3]. The Mittag-Leffler function plays a distinctive and challenging role in fractional calculus.
Researchers have devised numerical and approximative analytical methodologies for solving differential equations since analytical solutions are either unavailable or limited to complex mathematical functions [Citation4]. The range of numerical methods for differential equations is limited because of many of the challenges in their analysis, but the researcher has developed new numerical techniques. Classical initial conditions are used in fractional differential equations involving the Caputo fractional derivative, which are understandable and can be accurately measured. Researchers present several fractional derivatives, the most popular and often used non-classical derivatives of Caputo, Atangana-Baleanu, and Caputo-Fabrizio (C-F) [Citation5, Citation6]. There has been a recent development in kernels of fractional fractal derivatives in different forms, such as power law, exponential decay, and the Mittag-Leffler function [Citation7]. One advantage of these differential operators is that they depend on two orders: the fractional order and the fractal dimension. The fractal-fractional (F-F) operator, in the sense of Caputo and C-F, effectively conveys complexities that are difficult to depict using traditional fractional operators. The F-F approach has been used to solve numerous mathematical models. Many analysis fields, including banking [Citation8] and disease models [Citation9, Citation10], frequently employ the F-F operator [Citation11].
The Predictor-Corrector technique determines the numerical solution of fractional differential equations with several examples of this new fractional derivative [Citation12, Citation13]. This altered form of the Predictor-Corrector method non-uniform grid is an emphatic weapon for investigating fractional differential equations. Two dis-contiguous perspectives serve as the foundation for the methodologies employed in the numerical simulation of chaotic fractional order systems. ‘Time domain methods’ are techniques used to directly approximate the response of a fractional order system when solving differential equations numerically. The enhanced Adams-Bashforth Moulton algorithm, presented based on the predictor-corrector scheme, is one of the finest algorithms in this area [Citation14, Citation15]. The approach can effectively address the stability of non-linear systems of equations that easy-to-handle chaotic problems, which is highly significant. Various numerical methods that utilize fractional order operators to define these real-life issues differ and are customized.
There have been several numerical techniques developed to handle problems that involve fractional initial values. Mathematical modelling is a method used to depict reality or define present circumstances. It is also used to make predictions based on scientific principles about future possibilities. However, mathematical modelling has some constraints. It's challenging to estimate with many parameters and functions within constraints. The accuracy of mathematical modelling depends on how well it describes a particular issue [Citation16]. For instance, several variables can be ignored when simulating infectious diseases, such as weather and the health of each person in the community. Identifying the variables that impact modelling is essential to minimize errors and predict how the disease evolves.
Since the dawn of time, people have faced illnesses and combated them in their local communities. Science and technology have brought about rapid changes that make life more convenient and comfortable. Unfortunately, human progress has also led to significant deterioration of the natural environment in virtually every part of the world. Human growth in science and technology has caused many natural phenomena, such as cultural norms, modern agriculture, education, transportation, and others, to undergo significant changes [Citation17]. Advances in technology have made it possible to develop genetically changed food products. However, humans' overuse of natural resources is leading to the extinction of several plant and animal species. Attributing all the world's problems to pandemics when seasonal fluctuations significantly impact the dynamics of various species populations is unfair. How chronic wasting disease spreads among animals in a population and how environmental factors affect this illness [Citation18]. A mathematical model of non-integer predator and prey that includes two predators coexisting on a single prey in this article [Citation19].
Palm oil is the most commonly used vegetable oil worldwide and can be found in many packaged foods, stores, and cosmetics. It is an efficient and cost-effective oil that requires less land than other oils to produce the same amount of energy. However, the destruction of forests in some countries has led to widespread concern and government action because of the adverse effects on habitats and air quality. Illegal logging, which involves harvesting wood against the law, is a significant cause of deforestation, soil erosion, and biodiversity loss. This kind of logging can be carried out using dishonest means to acquire permission to use protected forests, cutting down protected species without a license, or exceeding the permitted amount of timber. These environmental problems can lead to even more severe catastrophes like climate change and ecological degradation. Forest degradation and deforestation are significant issues in developing countries, and the trade in illegally harvested wood products contributes to these problems. The soil's nutritional balance is crucial for maintaining productivity, as trees rely on nutrients from the soil for growth and maintenance. The paper's author proposes a mathematical framework for managing forest biomass, focussing on the example of Indonesia, which is affected by forest clearance because of mining and oil plantations [Citation20].
Most illicit cutting down, timber-prohibited activity, and the utilization of fossil fuels for cooking are to blame for the misappropriation of forestry biomass in underdeveloped countries. Domestic demand is one of the reasons that consume a sizable amount of illegal timber. An effort is made to build a mathematical model, stating tax as the control measure for lawful logging, taking into account the severity of the issue experienced by the prevalence of developing countries. Mathematical models are widely used when studying fisheries and forestry population conservation. The nineteenth century saw the introduction of research on forestry and fishing populations that considered economic factors. Taxes are a mechanism by which regulatory and governmental bodies modify the costs associated with access to forest resources. To encourage operators to maintain the value of the resources in the forest industry, keeping the tax rate in line with profitability incentives is essential. This will motivate them to improve their performance, even if they cannot meet the new requirements. In ecology, the functional response is crucial for illustrating the relationship between a predator and its prey. It has been demonstrated in the paper [Citation21] that the system undergoes Hopf-bifurcation when Holling type IV functional response and nonlinear age-selective prey harvesting utilizing the prey-generalist predator model are taken into account. Before now, many mathematical models for the utilization and preservation of forestry biomass were put out, taking into interpretation the impact of population, technology, and pollution. In this paper [Citation22], a theoretical framework and mathematical model that combines the population of humans with both reserved and unreserved forest biomass is put forth in response to the rise in atmospheric carbon dioxide.
A deeper understanding of complex systems, such as biomass dynamics, can be achieved using fractional derivatives, extending the concept of differentiation to non-integer orders. However, traditional calculus-based derivatives are still helpful as a reference point for comparison. In this study, we compare the effectiveness of these methods and investigate the dynamics of forestry biomass, building upon previous research. By doing so, we can better understand the role of forestry biomass in environmental sustainability, and the analysis yields new insights. F-F derivative provides a mathematical basis for characterizing the complex, self-similar nature of various elements in a forest, such as branches, leaves, and roots, at different observational levels. With the help of the fractional operator and the Lagrangian interpolation approach, forestry models can accurately depict the intricate configurations of biomass elements, leading to more reliable and precise biomass estimation. Using numerical methods also enables the measurement of these patterns across scales, enhancing the spatial accuracy of biomass models and accounting for the spatial heterogeneity, the variation in biomass distribution across different environmental conditions and locations, a defining characteristic of forest ecosystems. This approach offers a scalable and continuous way of describing the geographical distribution of biomass, particularly in complex forest environments. These methods provide flexibility to researchers, allowing them to adapt their methodologies to fit different forest types, observational scales, and modelling goals. The biomass model can use these techniques in big datasets and intricate spatial structures.
An ordered breakdown of the complete work is given in the following sections: concepts related to fractional calculus are defined in Section 2. Section 3 discusses the fractal-fractional model of forestry biomass. We address the model solution's stability, existence, and uniqueness in Sections 4 and 4.1. Numerical methods are presented in Section 5 for a fractional-order forestry biomass model. Numerical simulation and discussion are found in Section 6. The conclusion is found in Section 7.
2. Preliminaries
This part thoroughly looks at several definitions and theorems about fractional operators. There are also some fractal-fractional calculus-related outcomes. ψ denotes the fractional-order and φ represents the fractal-dimension.
Definition 2.1
[Citation23, Citation24]
We suppose that be an open interval, and be a continuous function and along with in the R-L derivative with power law kernel is (1) (1) with and
Definition 2.2
[Citation23, Citation24]
We suppose that be an open interval, and be a continuous function and along with in the R-L derivative with exponentially decaying kernel is (2) (2) and
Definition 2.3
[Citation23, Citation24]
We suppose that be an open interval, and be a continuous function and along with integral with power law kernel is (3) (3)
Definition 2.4
[Citation23, Citation24]
We suppose that be an open interval, and be a continuous function and along with integral with exponential decay kernel is (4) (4)
Theorem 2.1
[Citation25, Citation26]
Suppose that be an operator and completely continuous. Let The operator Θ has at least one fixed point or the set is unbounded.
3. The mathematical formulation of the forestry biomass system
Mathematical biology is a field that combines mathematical methods and biological research to create models. Although it's crucial to use mathematics in biology because of the increasingly quantitative nature of the field, the number of research groups and individuals working in this area is limited. Modern applications of mathematics in mathematical biology are fascinating and rapidly increasing. In the article [Citation27], we see author have suggested and studied some models while taking into interpretation various elements that take advantage of the natural population. We've taken into account the issue of illegal logging as one of the factors contributing to the absorption and degradation of forest biomass. When the author developed the model also took into reflection the approach put forth in this paper [Citation28]. They took into account alternative foods for predators. They found that in some circumstances is a shortage of prey, a rise in focal prey is seen if the predator uses fewer alternatives, and when the prey population approaches the carrying capacity, consumption of choices tends to zero. In developing the mathematical model, a similar idea of alternate resources in the development of forestry biomass is taken into account. In the proposed model, N denotes the nutrient density, and the forestry biomass is represented by P. E denotes the measure of effort intensity, and the industrial density is represented by the I. Regarding the nutrient-dependent growth rate, it is predicted that the biomass of the forest will increase logistically. The terms and represent legal and illicit tree-cutting and are directly affected by these trees cutting the loss of forestry biomass. Furthermore, it is considered that taxes levied by the government result in fewer industries, which help to preserve forest biomass. In addition, it is also believed that as taxes rise, the alternatives aid in the expansion of the industries denoted by the term ‘’ in the equation for I. Additionally, it is presumed that when there are alternatives, industries expand the fastest. Because the forestry population will be gotten near its carrying limit, alternate consumption tends to be zero. The state transition flow diagram for the model formulation is provided in Figure after carefully considering the assumptions. Q represents the constant input rate of nutrients, and β denotes the utilization rate of nutrients by forestry biomass. denotes the modified intrinsic growth rate of forestry biomass, and represents the washout rate of nutrients. K symbolizes the carrying capacity of forestry biomass, and γ signifies the new plantation rate. c denotes the fixed cost per unit of effort expended to harvest forestry biomass, and represents the illegal logging of forestry biomass. symbolizes the harvesting rate of forestry biomass, and signifies the reduction rate of industries. α means the rate of growth industries on account of efforts, and μ represents the maximum growth due to alternative industries. p and τ symbolize the fixed price and tax per unit of forestry biomass respectively.
The four differential equations listed below describe the model [Citation29]. (5) (5) the proposed model has the following initial conditions: where .
We apply the Fractal-Fractional (F-F) Caputo-Fabrizio (C-F) derivative [Citation10] in the model (Equation5(5) (5) ), we obtain (6) (6) the proposed model has the following initial conditions:
4. Existence and uniqueness of the solution
The existence and uniqueness theory [Citation30] for the forestry biomass model is constructed in this section. The system (Equation6(6) (6) ) is now rewrite, looking like this: (7) (7) where The system (Equation7(7) (7) ) rewrite as (8) (8) We employ the F-F integral in the C-F sense, we obtain (9) (9) where In the system (Equation7(7) (7) ), we employ the F-F integral in the C-F sense, we obtain (10) (10) We use the fixed point theory to prove the uniqueness and existence of the solution. A Banach space can be defined as follows = is space of all functions N, P, E, I respectively to a Banach space with the norm respectively to Consequently, the norm in the product space is created as Banach space is defined as We suppose the operator is defined in system (Equation10(10) (10) ), then we have (11) (11) where (12) (12)
Theorem 4.1
Let assume that the functions are continuous and constants are and such that ∀ where and n = 1, 2, 3, 4, then we get If the condition are fulfilled, then system (Equation6(6) (6) ) has a unique solution, where
Proof.
We suppose that and Let and we assume that be closed convex ball i.e We suppose that we have (13) (13) Similarly, we obtain (14) (14) By the use of Equations (Equation13(13) (13) ), (Equation14(14) (14) ) and Υ definition, we obtain (15) (15) when for any we obtain (16) (16) Similarly, we obtain (17) (17) use of Equations (Equation16(16) (16) ) and (Equation17(17) (17) ), we obtain (18) (18) Thus, is a contraction operator and Using the Banach contraction theorem and has a unique fixed point (i.e Leray-Schauder alternative theorem). Hence the proposed model (Equation6(6) (6) ) has a unique solution.
Theorem 4.2
We suppose that such that ∀ we get with and Other assumptions are where and then the system (Equation6(6) (6) ) has at least one solution.
Proof.
First, we suppose that be an operator and continuous. We can say that Θ operator is continuous because is continuous. Let be a bounded set and there ∃ constants such that, We have (19) (19) Similarly, we obtain (20) (20) Thus, using Equations (Equation19(19) (19) ) and (Equation20(20) (20) ), we proved that is uniformly bounded.
Thus, we show that Θ is equicontinuous. First, we suppose that we have (21) (21) Similarly, we obtain (22) (22) (23) (23) (24) (24) Therefore, is equicontinuous. Thus, is completely continuous.
Finally, we show that is a bounded. Let then When then and Then (25) (25) We reduce the complexity of the Equation (Equation25(25) (25) ), we have (26) (26) After using a similar procedure, we obtain (27) (27) The Equations (Equation26(26) (26) ) and (Equation27(27) (27) ) are now added, then we obtain (28) (28) Thus, we obtain Hence is bounded. Θ has at least one fixed point, by the theorem (2.1), therefore proposed model (Equation6(6) (6) ) has a solution.
4.1. Ulam-Hyer stability condtion
In this part, we define a few concepts and conditions of stability [Citation31] for the proposed model. First, we assume that a small perturbed parameter , and satisfy this condition Let (29) (29)
Lemma 4.1
The solution of the system if the following condition is true,
Theorem 4.3
We are using the systems (Equation7(7) (7) ) and (Equation8(8) (8) ) and Lemma 4.1 is U-H stable, if the solution of proposed model is U-H stable when where
Proof.
As we have shown, the proposed model has a unique solution. So first we assume that be solution and be any solution of proposed model, we have (30) (30) As a result, Equation (Equation30(30) (30) ) can be written as this indicates that the system's solution is stable.
5. Numerical methods for forestry biomass model
Managing nonlinearity in a biological model with fractional derivatives is a challenging task. Using a fractional model and dealing with nonlinearity is a problematic endeavour. Researchers have recently analysed biological models with a few novel numerical algorithms in recent years. These numerical methods help find an approximate solution to our problem. To discover an approximative solution, we first construct a fractional system.
5.1. Numerical method with F-F in Caputo sense
We apply the Fractal-Fractional Caputo derivative [Citation5, Citation31] in the model (Equation5(5) (5) ), we obtain (31) (31) We start with a power-law situation and use the model to use a numerical approach [Citation6]. Before the scheme begins, we formulate the proposed model regarding Volterra representation in the R-L sense. We reformulate the system (Equation31(31) (31) ) with the following results: (32) (32) After we apply the fractional integral to the system (Equation32(32) (32) ), we obtain (33) (33) For now, we solve the first equation of the system (Equation33(33) (33) ). Analogous solutions to the first equation are obtained for other equations as well. (34) (34) inserted in Equation (Equation34(34) (34) ), we obtain (35) (35) Reducing the complexity of Equation (Equation35(35) (35) ), we obtain (36) (36) by utilising the Lagrangian interpolation method to determine the approximate function in the interval into equation (Equation36(36) (36) ), we obtain (37) (37) We use the Equations (Equation37(37) (37) ) to (Equation36(36) (36) ), we obtain (38) (38) the following outcomes can be obtained by further solving equation (Equation38(38) (38) ), (39) (39) where and We apply the preceding analogous procedure to additional equations, and we obtain (40) (40) where and
5.2. Numerical method with F-F in C-F sense
The numerical technique [Citation6] in the C-F sense is being used with the following structure: (41) (41) After we apply the C-F fractional integral to the system (Equation41(41) (41) ), we obtain (42) (42) For now, we solve the first equation of the system (Equation42(42) (42) ). Analogous solutions to the first equation are obtained for other equations as well. (43) (43) inserted in Equation (Equation43(43) (43) ), we obtain (44) (44) Reducing the complexity of Equation (Equation44(44) (44) ), we obtain (45) (45) Applying the idea of Lagrange polynomial yields the following result: (46) (46) Reducing the complexity of Equation (Equation46(46) (46) ), we obtain (47) (47) We apply the preceding analogous procedure to additional equations, and we obtain (48) (48)
6. Numerical simulation and result discussion
Numerical simulations are presented to study the dynamic complexity of compartments of the proposed time-fractional order forestry biomass system. Our analysis used the Adams-Bashforth technique along with Caputo-Fabrizio and F-F Caputo operators. To demonstrate the effectiveness of the suggested strategy, simulations were performed using MATLAB [Citation32]. Table shows the parameter values used for simulations derived from the literature. Fractional derivatives and fractal dimensions in time can be used to examine various intrinsic properties of the forestry biomass model, compared to classical order models. As a first step, we simulate different fractional order values with fixed parameter values. As a next step, we investigate the effects of changing the sensitive parameter values. Considering different values of fractional order, it is crucial to investigate the impact of these key parameters on solution profiles for crucial decisions.
Figure shows the three-dimensional dynamic behaviour of the forestry biomass model (Equation31(31) (31) ) at different fractional orders and fractal dimensions using F-F Caputo derivative. Figure illustrates the two-dimensional dynamic behaviour of the forestry biomass model (Equation31(31) (31) ) at different fractional orders and fractal dimensions using F-F Caputo derivative. Figure represents the time series plots of the forestry biomass model (Equation31(31) (31) ) at different fractional orders and fractal dimensions using F-F Caputo derivative. We observed that the nutrient density fluctuated at varied rates for combinations of ψ and φ before reaching a stable equilibrium point during the simulation. In contrast, the forestry biomass showed a quick drop, a rise, and a slow reduction for different fractional orders and dimensions. Similar trends were seen in the effort intensity measure, which rose and fell over time on varying pathways depending on the fractional orders. The industrial density increases quickly and then stabilizes after some time. Further, when the fractal dimension is fixed , then we have analysed the behaviour of all state variables of proposed model. So, Figure shows the three-dimensional dynamic behaviour of the forestry biomass model (Equation31(31) (31) ) at different fractional orders. Figure illustrates the two-dimensional dynamic behaviour of the forestry biomass model (Equation31(31) (31) ) at different fractional orders. Figure represents the time series plots of the forestry biomass model (Equation31(31) (31) ) at different fractional orders. Additionally, as observed in this figure, the proposed model becomes rapidly stable when the fractional order ψ is changed compared to the integer model. Figure exhibits the time series diagrams of the forestry biomass model (Equation31(31) (31) ) at different fractal dimension orders and fractional order is fixed . Next, we analysed the forestry biomass model with the Caputo-Fabrizio derivative. Figure shows the three-dimensional dynamic behaviour of the forestry biomass model (Equation6(6) (6) ) at different fractional orders using Caputo-Fabrizio derivative. Figure illustrates the two-dimensional dynamic behaviour of the forestry biomass model (Equation6(6) (6) ) at different fractional orders using Caputo-Fabrizio derivative. Figure represents the time series plots of the forestry biomass model (Equation6(6) (6) ) at different fractional orders using Caputo-Fabrizio derivative. Because even a tiny adjustment to the sensitive parameter will result in a significant change, it is essential to estimate extremely delicate parameters correctly. Adjusting a less sensitive parameter won't significantly impact the system's behaviour and will take time to calculate. The investigation has led us to the sensitive factors, and . Moreover, we analysed the dynamic complexity of the different parameter values. So, Figure exhibits 3D dynamic behaviour of the forestry biomass model (Equation31(31) (31) ) at different illegal logging of forestry biomass . Figure illustrates the 2D dynamic behaviour of the forestry biomass model (Equation31(31) (31) ) at different illegal logging of forestry biomass Figure represents the time series plots of the forestry biomass model (Equation31(31) (31) ) at different illegal logging of forestry biomass . A decline in forestry biomass is observed when the rate of illegal logging increases; this relationship is illustrated in Figure . Figure exhibits 3D dynamic behaviour of the forestry biomass model (Equation31(31) (31) ) at different harvesting rate of forestry biomass . Figure illustrates the 2D dynamic behaviour of the forestry biomass model (Equation31(31) (31) ) at different harvesting rate of forestry biomass Figure represents the time series plots of the forestry biomass model (Equation31(31) (31) ) at different harvesting rate of forestry biomass . A decline in forestry biomass is observed when the rate of harvesting increases; this relationship is illustrated in Figure .
7. Conclusion
In this research, we have demonstrated how illegal logging and tax imposition affect the utilization and preservation of forestry biomass. The proposed model study also assesses the impact of alternative resources on the development of industries. This research aims to analyse forestry biomass system dynamics utilizing effective approximation techniques for the F-F derivatives in the Caputo and C-F sense. Among the most fundamental features of these versions is the concept of memory. The F-F idea can predict the system complexity in two orders more precisely. First, we prove the existence and uniqueness of the fractional-order system and then the numerical solution of it through numerical simulation. The numerical analysis demonstrates that forestry biomass can be preserved more economically when government organizations implement the tax. Numerical simulations show more powerful converging effects of the forestry biomass model for different fractional orders resembling integer orders. The fractal fractional operator has been used in two senses to obtain the numerical results for the non-integer model presented here. The analysis shows that forestry biomass increases more quickly in the presence of tax than in the absence of tax. We identify sensitive parameters in the proposed model and analyse those to examine the behaviour of the numerical solution and its sensitivity to different transmission parameters. It has been shown that even a slight change in the order of derivatives and the system parameters can significantly influence the behaviour of a graphical representation. The techniques presented in this research can be extended along with other operators. In the future, it can be beneficial if we resolve the proposed model with the help of other different numerical techniques and compare the numerical solutions. Finding the best way for forestry biomass conservation and the optimal control policy applied.
Acknowledgement
Researchers Supporting Project number (RSPD2024R526), King Saud University, Riyadh, Saudi Arabia.
Authors contributions
All authors contributed equally and significantly in writing this paper and typed,read, and approved the final manuscript.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Data availability statements
Data available on request from the authors.
References
- Hilfer R. Applications of fractional calculus in physics. Universität Mainz & Universität Stuttgart, Germany: World Scientific; 2000. p. 472. doi: 10.1142/3779
- Machado JT, Kiryakova V, Mainardi F. Recent history of fractional calculus. Commun Nonlinear Sci Numer Simul. 2011;16(3):1140–1153. doi: 10.1016/j.cnsns.2010.05.027
- Sabatier JATMJ, Agrawal OP, Machado JAT. Advances in fractional calculus. Vol. 4. Dordrecht: Springer; 2007. doi: 10.1007/978-1-4020-6042-7
- Das S. Functional fractional calculus. Vol. 1. Berlin Heidelberg: Springer-Verlag; 2011. doi: 10.1007/978-3-642-20545-3. ISBN 978-3-642-20544-6
- Atangana A, Araz S.I. New numerical approximation for chua attractor with fractional and fractal-fractional operators. Alex Eng J. 2020;59(5):3275–3296. doi: 10.1016/j.aej.2020.01.004
- Khan MA, Atangana A, Muhammad T, et al. Numerical solution of a fractal-fractional order chaotic circuit system. Rev Mex Fís. 2021;67(5):5. doi: 10.31349/RevMexFis.67.051401
- Atangana A. Fractal-fractional differentiation and integration: connecting fractal calculus and fractional calculus to predict complex system. Chaos Solitons Fract. 2017;102:396–406. doi: 10.1016/j.chaos.2017.04.027
- Wang W, Khan MA. Analysis and numerical simulation of fractional model of bank data with fractal–fractional atangana–baleanu derivative. J Comput Appl Math. 2020;369:112646. doi: 10.1016/j.cam.2019.112646
- Ghanbari B. On the modeling of the interaction between tumor growth and the immune system using some new fractional and fractional-fractal operators. Adv Differ Equ. 2020;2020(1):585. doi: 10.1186/s13662-020-03040-x
- Kumar P, Kumar A, Kumar S, et al. A fractional order co-infection model between malaria and filariasis epidemic. Arab J Basic Appl Sci. 2024;31(1):132–153.
- Khan MA. The dynamics of dengue infection through fractal-fractional operator with real statistical data. Alex Eng J. 2021;60(1):321–336. doi: 10.1016/j.aej.2020.08.018
- Diethelm K, Ford NJ, Freed AD. A predictor-corrector approach for the numerical solution of fractional differential equations. Nonlinear Dyn. 2002;29:3–22. doi: 10.1023/A:1016592219341
- Daftardar-Gejji V, Sukale Y, Bhalekar S. A new predictor–corrector method for fractional differential equations. Appl Math Comput. 2014;244:158–182.
- Ali Z, Rabiei F, Shah K, et al. Fractal-fractional order dynamical behavior of an hiv/aids epidemic mathematical model. Eur Phys J Plus. 2021;136(1):36. doi: 10.1140/epjp/s13360-020-00994-5
- Abro KA. Numerical study and chaotic oscillations for aerodynamic model of wind turbine via fractal and fractional differential operators. Numer Methods Partial Differ Equ. 2022;38(5):1180–1194. doi: 10.1002/num.v38.5
- Baleanu D, Diethelm K, Scalas E, et al. Fractional calculus: models and numerical methods. Vol. 3. Singapore: World Scientific; 2012. p. 428.
- Magin R. Fractional calculus in bioengineering, part 1. Crit Rev Biomed Eng. 2004;32(1):104. doi: 10.1615/critrevbiomedeng.v32.i1.10
- Kumar P, Kumar A, Kumar S. A study on fractional order infectious chronic wasting disease model in deers. Arab J Basic Appl Sci. 2023;30(1):601–625.
- Kumar A, Kumar S, Momani S, et al. A chaos study of fractal–fractional predator–prey model of mathematical ecology. Math Comput Simul. 2023. doi: 10.1016/j.matcom.2023.09.010. Corpus ID: 263718777
- Pratama MAA, Zikkah RN, Anggriani N, et al. A mathematical model to study the effects of population pressure on two-patch forest resources. In: AIP Conference Proceedings. Vol. 2264, AIP Publishing LLC; 2020. p. 050003.
- Datta J, Jana D, Upadhyay RK. Bifurcation and bio-economic analysis of a prey-generalist predator model with holling type iv functional response and nonlinear age-selective prey harvesting. Chaos Solitons Fract. 2019;122:229–235. doi: 10.1016/j.chaos.2019.02.010
- Devi S, Mishra RP. Preservation of the forestry biomass and control of increasing atmospheric CO 2 using concept of reserved forestry biomass. Int J Appl Comput Math. 2020;6(1):17. doi: 10.1007/s40819-019-0767-z
- Atangana A, Qureshi S. Modeling attractors of chaotic dynamical systems with fractal–fractional operators. Chaos Solitons Fract. 2019;123:320–337. doi: 10.1016/j.chaos.2019.04.020
- Z Liu ZL, Khan MA. Fractional investigation of bank data with fractal-fractional caputo derivative. Chaos Solitons Fract. 2020;131:109528. doi: 10.1016/j.chaos.2019.109528
- Ali Z, Shah K, Zada A, et al. Mathematical analysis of coupled systems with fractional order boundary conditions. Fractals. 2020;28(08):2040012. doi: 10.1142/S0218348X20400125
- Granas A, Dugundji J. Fixed point theory. Vol. 14. New York: Springer-Verlag; 2003. ISBN-10:0-387-00173-5
- Chaudhary M, Dhar J, Misra OP. A mathematical model for the conservation of forestry biomass with an alternative resource for industrialization: a modified leslie gower interaction. Model Earth Syst Environ. 2015;1:1–10. doi: 10.1007/s40808-015-0056-8
- Ghosh B, Kar TK. Sustainable use of prey species in a prey–predator system: jointly determined ecological thresholds and economic trade-offs. Ecol Modell. 2014;272:49–58. doi: 10.1016/j.ecolmodel.2013.09.013
- Chaudhary M, Dhar J, Misra OP. Deleterious effects of legal and illegal logging on forestry biomass using a dynamical model. In: AIP Conference Proceedings. Vol. 2451. AIP Publishing; 2022.
- Ali Z, Rabiei F, Hosseini K. A fractal–fractional-order modified predator–prey mathematical model with immigrations. Math Comput Simul. 2023;207:466–481. doi: 10.1016/j.matcom.2023.01.006
- Xu C, Saifullah S, Ali A, et al. Theoretical and numerical aspects of rubella disease model involving fractal fractional exponential decay kernel. ‘Res Phys. 2022;34:105287.
- The MathWorks Inc. Matlab version: 9.0 (r2016a). 2016.