Next Article in Journal
Evolution of Generalized Brans–Dicke Parameter within a Superbounce Scenario
Previous Article in Journal
A Transition Model in f(R,T) Theory via Observational Constraints
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Post-Pandemic Sector-Based Investment Model Using Generalized Liouville–Caputo Type

by
Muath Awadalla
1,*,
Muthaiah Subramanian
2,*,
Prakash Madheshwaran
2 and
Kinda Abuasbeh
1
1
Department of Mathematics and Statistics, College of Science, King Faisal University, Hofuf 31982, Al Ahsa, Saudi Arabia
2
Department of Mathematics, KPR Institute of Engineering and Technology, Coimbatore 641407, India
*
Authors to whom correspondence should be addressed.
Symmetry 2023, 15(4), 789; https://doi.org/10.3390/sym15040789
Submission received: 6 March 2023 / Revised: 21 March 2023 / Accepted: 22 March 2023 / Published: 24 March 2023
(This article belongs to the Section Mathematics)

Abstract

:
In this article, Euler’s technique was employed to solve the novel post-pandemic sector-based investment mathematical model. The solution was established within the framework of the new generalized Caputo-type fractional derivative for the system under consideration that serves as an example of the investment model. The mathematical investment model consists of a system of four fractional-order nonlinear differential equations of the generalized Liouville–Caputo type. Moreover, the existence and uniqueness of solutions for the above fractional order model under pandemic situations were investigated using the well-known Schauder and Banach fixed-point theorem technique. The stability analysis in the context of Ulam—Hyers and generalized Ulam—Hyers criteria was also discussed. Using the investment model under consideration, a new analysis was conducted. Figures that depict the behavior of the classes of the projected model were used to discuss the obtained results. The demonstrated results of the employed technique are extremely emphatic and simple to apply to the system of non-linear equations. When a generalized Liouville–Caputo fractional derivative parameter ( ρ ) is changed, the results are asymmetric. The current work can attest to the novel generalized Caputo-type fractional operator’s suitability for use in mathematical epidemiology and real-world problems towards the future pandemic circumstances.

1. Introduction

Investments from stock markets and international companies are very important to the growth and strength of an economy. Leading countries around the world make many offers to show what startups can do for investors. Multi-dimensional sector investments are famous around the world due to the minimal sectional downturn of particular sectors. Pandemic situations such as severe acute respiratory syndrome (SARS CoV-1), Middle Eastern respiratory syndrome (MERS), and the ongoing SARS CoV-2 (COVID-19), followed by influenza viruses causing vnfluenza AH1N1 pdm 2009 significantly affected normal situations and made them worse. Currently, the world is in a post-pandemic era. Recovering from the downturned situation with new strategies and adaptations makes it happen. The intention of this work is to study the fractional investment model with regards to the future pandemic situation and attempt to derive some helpful clues that can be strategically used in the future. Saudi Arabia holds the 34th place in average development, with a rate of 6.29% [1]. Saudi Arabia holds a wide range of investment opportunities in numerous companies and industries under various sectors [2]. Recently, Javid et al. [3] investigated the determinants of short- and long-run private investment behaviour in Saudi Arabia for eight non-oil sectors. A moderating structural equation modeling-based model was studied on Saudi Arabia by Mohamed Ali Shabeeb Ali et al. [4]. This could be the first attempt in using fractional derivatives in an investing model according to the author’s knowledge. Researchers are familiar with the limitations of integer order calculus, particularly when used to analyse phenomena related to diffusion, hereditary traits, long-range waves, history-based phenomena, etc. The theory of fractional calculus has been around for more than 300 years, but it is still relevant today for solving problems in the real world. In research on fractional calculus, many fractional derivatives have been found. Caputo, Caputo–Fabrizio, and Atangana–Baleanu are the fractional derivatives that are used most often in many different fields. Generalized Caputo derivatives, a recently proposed derivative with properties resembling those of Caputo derivatives, are being used to analyze the findings in the current work [5]. Several investigations have used fractional order derivatives, and at the moment, mathematics-based epidemiology is one of their main fields of application. Numerous domains, including engineering, physics, signal and image processing, mechanics and dynamical systems, biology, control theory, and environmental sciences, have presented a number of non-integer order derivative models where non-locality plays a significant role [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]. Because the Caputo derivative is appropriate for initial value problems (IVPs) and shares many traits with integer-order derivatives, it has been used in applications of fractional calculus to mobilize a large number of physical issues. Feng et al. [21] presented a study on the general behaviour of the Maxwell mechanical model using the combined Caputo fractional derivative. The generalized Caputo-type fractional derivative’s nature shares several characteristics with the Caputo derivative. The projected model is solved in the current work using the Euler scheme. The value of the parameters has a significant impact on the behaviour of the generalized fractional integral operator, making it a useful tool for manipulating and creating mathematical models in applications of fractional calculus. The novel generalized Caputo-type fractional derivative includes additional qualities when compared to existing fractional derivatives such as the Caputo, Caputo–Fabrizio, and Atangana–Baleanu fractional derivatives. There is another parameter ρ in addition to the fractional order parameter that is notably helpful in graphical simulations with regard to real data. By changing the parameter value ρ , we can see more types of graphs. In other words, when classical derivatives are used, the classical derivatives are local. A fractional derivative is nonlocal. Because of this property, these derivatives are suitable for modeling more physical phenomena, such as earthquake vibrations and polymers, among many others. The structure of this document is as follows. Preliminary definitions are provided in Section 2. The generalised Liouville–Caputo interpretation of the fractional order model for sector-based investment is developed in Section 3. The model’s existence and uniqueness of solutions are established in Section 4. The stability analysis of the model’s solution in the Ulam–Hyers and generalised Ulam–Hyers segments is provided in Section 5. The numerical approach and numerical simulations that demonstrate the theoretical results are presented in Section 6. Lastly, Section 7 provides a conclusion.

2. Preliminaries

For our investigation, we recollect some rough definitions of generalized fractional derivatives and fractional integrals.
The space of all complex-valued Lebesgue measurable functions ψ on ( a , b ) equipped with the norm is denoted by X c z ( a , b ) :
ψ X c z = a b x c ψ ( x ) z d x x 1 z < , c R , 1 z .
Let L 1 ( a , b ) represent the space of all Lebesgue measurable functions ϱ on ( a , b ) endowed with the norm:
ϱ L 1 = a b ϱ ( x ) d x < .
We further recall that A C n ( H , R ) = { x : H R : x , x , , x ( n 1 ) C ( H , R ) and x ( n 1 ) } is absolutely continuous. For 0 ε < 1 , we define C ε , ρ ( H , R ) = { g : H R : ( t ρ a ρ ) ε g ( t ) C ( H , R ) endowed with the norm g C ε , ρ = ( t ρ a ρ ) ε g ( t ) C . Moreover, we define the class of functions g that have an absolute continuous γ n 1 derivative, denoted by A C γ n ( H , R ) , as follows: A C γ n ( H , R ) = { g : H R : γ n 1 g AC ( H , R ) , γ = t 1 ρ d d t } , which is equipped with the norm g C γ , ε n = k = 0 n 1 γ k g C + γ n g C ε , ρ and is defined by
C γ , ε n ( H , R ) = g : H R : γ n 1 g C ( H , R ) , γ n g C ε , ρ ( H , R ) , γ = t 1 ρ d d t .
Notice that C γ , 0 n = C γ n .
Definition 1
([22]). The generalized fractional integrals of g X d p ( a , b ) of order υ > 0 and ρ > 0 for < a < t < b < are defined as follows:
( ρ I a + υ g ) ( t ) = ρ 1 υ Γ ( υ ) a t σ ρ 1 ( t ρ σ ρ ) 1 υ g ( σ ) d σ ,
( ρ I b υ g ) ( t ) = ρ 1 υ Γ ( υ ) t b σ ρ 1 ( t ρ σ ρ ) 1 υ g ( σ ) d σ ,
Definition 2
([23]). The generalized fractional derivatives (GFDs) that are associated with GFIs (1) and (2) for 0 a < t < b < are defined as follows:
( ρ D a + υ g ) ( t ) = t 1 ρ d d t n ( ρ I a + n υ g ) ( t ) = ρ υ n + 1 Γ ( n υ ) t 1 ρ d d t n a t σ ρ 1 ( t ρ σ ρ ) υ n + 1 g ( σ ) d σ ,
( ρ D b υ g ) ( t ) = t 1 ρ d d t n ( ρ I b n υ g ) ( t ) = ρ υ n + 1 Γ ( n υ ) t 1 ρ d d t n t b σ ρ 1 ( t ρ σ ρ ) υ n + 1 g ( σ ) d σ ,
if the integrals exist.
Definition 3
([24]). The above GFDs define the left- and right-sided generalized Liouville–Caputo-type fractional derivatives of g A C δ n [ a , b ] of order υ 0
C ρ D a + υ g ( x ) = ρ D a + υ g ( t ) j = 0 n 1 δ j g ( a ) j ! t ρ a ρ ρ j ( x ) , δ = x 1 ρ d d x ,
C ρ D b υ g ( x ) = ρ D b υ g ( t ) j = 0 n 1 ( 1 ) j δ j g ( b ) j ! b ρ t ρ ρ j ( x ) , δ = x 1 ρ d d x ,
when n = [ υ ] + 1 .
Lemma 1
([24]). Let g AC δ n [ a , b ] or C δ n [ a , b ] and υ R . Then,
ρ I a + υ C ρ D a + υ g ( x ) = g ( x ) j = 0 n 1 δ j g ( a ) j ! x ρ a ρ ρ j ,
ρ I b υ C ρ D b υ g ( x ) = g ( x ) j = 0 n 1 ( 1 ) j δ j g ( b ) j ! b ρ x ρ ρ j .
In particular, for 0 < υ 1 , we have
ρ I a + υ C ρ D a + υ g ( x ) = g ( x ) g ( a ) , ρ I b υ C ρ D b υ g ( x ) = g ( x ) g ( b ) .

3. Formulation of the Model

Let an entire sector’s populace be S ( t ) , which may be listed under one of three categories. S N ( t ) , or sectors that were not affected, were not directly connected to the pandemic and were thus not affected during the pandemic. Even the chosen alternative options to sustain these sectors proved to be more beneficial than before the pandemic for this category of sectors, which includes IT, social media, and online education. S D ( t ) , or downturned sectors, were classified based on the serious impact of pandemic on these sectors, which made their situation fundamentally more difficult. Examples such as tourism, real estate, manufacturing, stock exchange, self-employment, etc. fall into this group. S N D ( t ) , or not-downturned sectors, were vigorously influenced by the pandemic, but this influence worked in a positive way and led to significant growth. Examples of this include hospitals, pharmaceutical companies, and telecommunications companies. A diagram of the dynamics of the sector-based investment model is shown in Figure 1. This model was made for the generalized Liouville–Caputo type. The classical Caputo type has been utilized for many physical interpretations in mathematical modeling problems. Considering this, in order to justify our model, comparison graphs (Figure 2, Figure 3, Figure 4 and Figure 5) have been included in the classical Caputo sense ( υ = 0.95 , ρ = 1 ) and the generalized Liouville—Caputo sense υ = 0.95 , ρ = 0.92 , 0.93 , 0.94 , 0.95 ) for S ( t ) , S D ( t ) , S N D ( t ) , and S N ( t ) , which exhibit excellent correlation and give us confidence to proceed with this model. This was a fresh attempt towards creating an investment model in fractional perspectives; the proposed investment model consists of a system of fractional order differential equations, which were inspired and improvised from the SEIR epidemic model for COVID-19 transmission by the Caputo derivative of fractional order [25]. The detailed correlations of factors influencing sector-based investments are explained in the form of a flow chart.
The dynamics of the sectors were represented through the system of fractional order differential equations.
c ρ D 0 + υ S ( t ) = Φ ( η 1 S N D ( t ) + η 2 S D ( t ) ) S ( t ) ζ ^ S ( t ) , c ρ D 0 + υ S D ( t ) = λ S D ( t ) + ( η 2 S D ( t ) ) S ( t ) + ζ ^ S N D ( t ) ζ ^ S D ( t ) χ S D ( t ) , c ρ D 0 + υ S N D ( t ) = λ S N D ( t ) + ( η 1 S N D ( t ) ) S ( t ) + ζ ^ S D ( t ) ζ ^ S N D ( t ) , c ρ D 0 + υ S N ( t ) = χ S D ( t ) ζ ^ S N ( t ) ,
where the initial conditions are S ( 0 ) = S 0 > 0 , S D ( 0 ) = S D 0 > 0 , S N D ( 0 ) = S N D 0 > 0 , and S N ( 0 ) = S N 0 > 0 .

Descriptions of Constraints

  • Φ ( n × N ) —enrollment rate;
  • N—total number of sectors and enrollment rate, where n is the sector-introducing rate;
  • ζ ^ —usual loss rate of the sector;
  • η 1 —rate of transfer of affected sectors not down turned;
  • η 2 —rate of transfer of affected and sectors down turned;
  • χ —rate of recovery from the pandemic;
  • λ —rate of transfer of the sector.
Consequently, the proposed model has the following structure:
c ρ D 0 + υ ψ ( t ) = G ( t , ψ ( t ) ) , t H : = [ 0 , T ] , 0 < υ 1 , ρ > 0 , ψ ( 0 ) = ψ 0 0 ,
on the condition that
ψ ( t ) = ( S , S D , S N D , S N ) T , ψ ( 0 ) = ( S 0 , S D 0 , S N D 0 , S N 0 ) T , G ( t , ψ ( t ) ) = ( Λ i ( S , S D , S N D , S N ) ) T , i = 1 , , 5 ,
where ( · ) T represents the transpose operation. Given Lemma 1, the integral form of the problem (8), which is identical to the model (7), is provided by
ψ ( t ) = ψ 0 + ρ I 0 + υ G ( t , ψ ( t ) ) = ψ 0 + ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ( σ ) ) d σ .

4. Existence and Uniqueness Results

Let K = C ( [ 0 , T ] , R ) denote the Banach space of all continuous functions from [ 0 , T ] to R endowed with the norm defined by
ψ = sup t H ψ ( t ) ,
where
ψ ( t ) = S ( t ) , S D ( t ) , S N D ( t ) , S N ( t ) ,
and S , S D , S N D , S N C ( [ 0 , T ] , R ) .
Theorem 1.
Suppose that the function G C ( H , R ) maps a bounded subset of H × R 4 into relatively compact subsets of R . In addition, there exists a constant M G > 0 such that
( F 1 ) G ( t , ψ 1 ( t ) ) G ( t , ψ 2 ( t ) ) M G ψ 1 ( t ) ψ 2 ( t ) for all t H and each ψ 1 , ψ 2 C ( H , R ) . Then, Problem (8) has a unique solution provided that Ψ M G < 1 , where
Ψ = T ρ υ ρ υ γ ( υ + 1 ) .
Proof. 
Consider the operator Q : K K defined by
( Q ψ ) ( t ) ) = ψ 0 + ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ( σ ) ) d σ .
Obviously, the operator Q is well defined, and the unique solution of model (7) is just the fixed point of Q . Indeed, let us take sup t H G ( t , 0 ) = P 1 and ω ψ 0 + Ψ P 1 . Thus, it is enough to show that Q B ω B ω , where the set B ω = { ψ K : ψ ω } is closed and convex. Now, for any ψ B ω , it yields
( Q ψ ) ( t ) ψ 0 + ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ( σ ) ) d σ ψ 0 + ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ( σ ) ) G ( σ , 0 ) G ( σ , 0 ) d σ ψ 0 + ( M G ω + P 1 ) ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ d σ ψ 0 + Ψ ( M G ω + P 1 ) ω .
Hence, the results follow. Additionally, given any ψ 1 , ψ 2 , we obtain
( Q ψ 1 ) ( t ) ( Q ψ 2 ) ( t ) ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ 1 ( σ ) ) G ( σ , ψ 2 ( σ ) ) d σ M G ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ ψ 1 ( σ ) ψ 2 ( σ ) d σ Ψ M G ψ 1 ( σ ) ψ 2 ( σ ) ,
which implies that ( Q ψ 1 ) ( Q ψ 2 ) Ψ M G ψ 1 ψ 2 . As a result of the Banach fixed-point theorem [26], the proposed model (7) possesses a unique solution. Using the concept of Schauder’s fixed-point theorem [26], we then demonstrate the existence of solutions for Problem (8) that are equivalent to the suggested model (7). Hence, the following premise is required.
( F 2 ) Suppose that there exist κ 1 , κ 2 K such that
G ( t , ψ ( t ) ) κ 1 ( t ) + κ 2 ψ ( t ) , f o r a n y ψ K , t H ,
such that κ 1 ^ = sup t H κ 1 ( t ) , κ 2 ^ = sup t H κ 2 ( t ) < 1 .
Lemma 2.
The operator Q defined in (11) is completely continuous.
Proof. 
The continuity of the function G gives the continuity of the operator Q . Therefore, for any ψ B ω , where B ω is defined above, we obtain
( Q ψ ) ( t ) ψ 0 + ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ( σ ) ) d σ ψ 0 + ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ( σ ) ) d σ ψ 0 + ( κ 1 ^ + κ 2 ^ ψ ) ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ d σ ψ 0 + Ψ ( κ 1 ^ + κ 2 ^ ψ ) < + .
Therefore, the operator Q is uniformly bounded. Next, we prove the equicontinuity of Q . To do so, we let sup ( t , ψ ) H × B ω G ( t , ψ ( t ) ) = G . Then, for any t 1 , t 2 H such that t 2 t 1 , it gives
( Q ψ 1 ) ( t ) ( Q ψ 2 ) ( t ) = ρ 1 ( υ ) G Γ ( υ + 1 ) 0 t 1 σ ρ 1 ( t 2 ρ σ ρ ) 1 υ σ ρ 1 ( t 1 ρ σ ρ ) 1 υ d θ + t 1 t 2 σ ρ 1 ( t 2 ρ σ ρ ) 1 υ d θ G ρ υ Γ ( υ + 1 ) [ 2 ( t 2 ρ t 1 ρ ) υ ( t 2 ρ υ t 1 ρ υ ) ] 0 as t 2 t 1 .
Hence, the operator Q is equicontinuous and is thus relatively compact on B ω . Therefore, as a consequence of the Arzela–Ascoli theorem [26], Q is completely continuous. □
Theorem 2.
Suppose that the function G : H × R 4 R is continuous and satisfies the assumption ( F 2 ) . Then, Problem (8). which is equivalent to the proposed Model (7), has at least one solution.
Proof. 
We define a set V = { ψ K : ψ = ξ ( Q ψ ) ( t ) , 0 < ξ < 1 } . Clearly, in view of Lemma 2, the operator Q : V K as defined in (11) is completely continuous. Now, for any assumption ( F 2 ) , it yields
( ψ ) ( t ) = ξ ( Q ψ ) ( t ) ψ 0 + ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ( σ ) ) d σ ψ 0 + ( κ 1 ^ + κ 2 ^ ψ ) ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ d σ ψ 0 + Ψ ( κ 1 ^ + κ 2 ^ ψ ) < + .
Thus, the set V is bounded. Thus, the operator Q has at least one fixed point, which is the proposed model’s solution (7). □

Asymmetric Result

Remark 1.
The generalized Caputo-type investment model reduces to a Caputo sense when ρ = 1 .
c D 0 + υ S ( t ) = Φ ( η 1 S N D ( t ) + η 2 S D ( t ) ) S ( t ) ζ ^ S ( t ) , c D 0 + υ S D ( t ) = λ S D ( t ) + ( η 2 S D ( t ) ) S ( t ) + ζ ^ S N D ( t ) ζ ^ S D ( t ) χ S D ( t ) , c D 0 + υ S N D ( t ) = λ S N D ( t ) + ( η 1 S N D ( t ) ) S ( t ) + ζ ^ S D ( t ) ζ ^ S N D ( t ) , c D 0 + υ S N ( t ) = χ S D ( t ) ζ ^ S N ( t ) .
where the initial conditions are S ( 0 ) = S 0 > 0 , S D ( 0 ) = S D 0 > 0 , S N D ( 0 ) = S N D 0 > 0 , and S N ( 0 ) = S N 0 > 0 .
Consequently, the proposed model has the following structure:
c D 0 + υ ψ ( t ) = G ( t , ψ ( t ) ) , t H : = [ 0 , T ] , 0 < υ 1 , ψ ( 0 ) = ψ 0 0 ,
on the condition that
ψ ( t ) = ( S , S D , S N D , S N ) T , ψ ( 0 ) = ( S 0 , S D 0 , S N D 0 , S N 0 ) T , G ( t , ψ ( t ) ) = ( Λ i ( S , S D , S N D , S N ) ) T , i = 1 , , 5 ,
where ( · ) T represents the transpose operation. Given Lemma 1, the integral form of Problem (13), which is identical to Model (12), is provided by
ψ ( t ) = ψ 0 + I 0 + υ G ( t , ψ ( t ) ) = ψ 0 + 1 Γ ( υ ) 0 t ( t σ ) υ 1 G ( σ , ψ ( σ ) ) d σ .

5. Stability Results

This section derives the stability of the suggested model (7) in terms of Ulam—Hyers stability and generalised Ulam—Hyers stability. Ulam introduced the concept of Ulam stability in [27,28]. The aforementioned stability was then explored in various research publications on classical fractional derivatives, for example, [29,30,31,32]. Additionally, as the generalised stability of the suggested model (7) is essential for an approximate solution, we work to employ nonlinear functional analysis on Ulam–Hyer stability. As a result, the following definitions are required. Consider the following inequality if θ > 0 :
D 0 + υ c ρ ψ ^ ( t ) G ( t , ψ ^ ( t ) ) θ , t H ,
where θ = max ( θ j ) T , j = 1 , , 5 .
Definition 4.
The proposed problem (8), which is equivalent to Model (7), is Ulam—Hyers stable if there exists E G > 0 such that, for every θ > 0 and for each solution ψ ^ K satisfying Inequality (16), there exists a solution ψ K of Problem (8), with
ψ ^ ( t ) ψ ( t ) E G θ , t H , w h e r e E G = max ( E G j ) T .
Definition 5.
Problem (8), which is equivalent to Model (7), is referred to as being generalized Ulam—Hyers stable if there exists a continuous function δ G : R + R + , with δ G ( 0 ) = 0 such that, for each solution ψ ^ K of Inequality (16), there exists a solution ψ K of Problem (8) such that
ψ ^ ( t ) ψ ( t ) δ G θ , t H , w h e r e δ G = max ( δ G j ) T .
Remark 2.
A function ψ ^ K is a solution of Inequality (16) if and only if there exists a function g K with the following property:
  • (i) g ( t ) θ , g = max ( g j ) T , t H ;
  • (ii) D 0 + υ c ρ ψ ^ ( t ) = G ( t , ψ ^ ( t ) ) + g ( t ) , t H .
Lemma 3.
Assume that ψ ^ K satisfies Inequality (16); then, ψ ^ satisfies the integral inequality described by
ψ ^ ( t ) ψ ^ 0 ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ^ ( σ ) ) d σ Ψ θ .
Proof. 
With the help of Remark 1, c ρ D 0 + υ ψ ^ ( t ) = G ( t , ψ ^ ( t ) ) + g ( t ) , and Lemma 1 gives
ψ ^ ( t ) = ψ ^ 0 + ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ^ ( σ ) ) d σ + ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ g ( σ ) d σ .
Using (i) of Remark 1, we obtain
ψ ^ ( t ) ψ ^ 0 ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ^ ( σ ) ) d σ ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ g ( σ ) d σ Ψ θ .
The proof is completed. □
Theorem 3.
Suppose that G : H × R 4 R is continuous for every ψ K and assumption ( F 1 ) holds with 1 Ψ M G > 0 . Thus, Problem (8), which is equivalent to Model (7) is Ulam—Hyers and, consequently, generalized Ulam—Hyers stable.
Proof. 
Suppose that ψ ^ K satisfies Inequality (16) and ψ K is a unique solution of Problem (8). Thus, for any θ > 0 , t H , and Lemma 3, it gives
ψ ^ ( t ) ψ ( t ) = max t H ψ ^ ( t ) ψ 0 ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ( σ ) ) d σ ψ ^ ( t ) ψ ^ 0 ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ^ ( σ ) ) d σ + max t H ρ 1 υ Γ ( υ ) 0 t σ ρ 1 ( t ρ σ ρ ) 1 υ G ( σ , ψ ^ ( σ ) ) G ( σ , ψ ( σ ) ) d σ Ψ θ + Ψ M G ψ ^ ( t ) ψ ( t ) .
Therefore,
ψ ^ ψ E G θ ,
where
E G = Ψ 1 Ψ M G .
Setting δ G ( θ ) = E G θ such that δ G ( 0 ) = 0 , we conclude that the proposed Problem (8) is both Ulam—Hyers and generalized Ulam—Hyers stable. □

6. Numerical Simulation

We consider System (7) in a compact form as follows:
c ρ D t υ z ( t ) = φ ( t , z ( t ) ) , z ( 0 ) = z 0 , 0 t T < ,
where z = ( S , S D , S N D , S R ) R + 4 , z 0 = ( S 0 , S D 0 , S N D 0 , S R 0 ) is the initial vector, and z ( t ) R is a continuous vector function satisfying the Lipschitz condition
φ ( z 1 ( t ) ) φ ( z 2 ( t ) ) γ z 1 ( t ) z 2 ( t ) , γ > 0 ,
Applying a fractional integral operator corresponding to the generalized Caputo derivative to Equation (17), we obtain
z ( t ) = z 0 + ρ I υ φ ( z ( t ) ) , 0 T < .
Set h = T ρ N and t n = n h , where t [ 0 , T ] , N is a natural number, and n = 0 , 1 , 2 , , N . Let z n be the approximation of z ( t ) at t = t n . The generalized Liouville–Caputo fractional derivative operator’s governing model has the following numerical technique:
S n + 1 = S 0 + h υ ρ υ γ ( υ + 1 ) k = 0 n ( ( n + 1 k ) υ ( n k ) υ ) ϑ 1 ( t k , z k ) , S D n + 1 = S D 0 + h υ ρ υ γ ( υ + 1 ) k = 0 n ( ( n + 1 k ) υ ( n k ) υ ) ϑ 2 ( t k , z k ) , S N D n + 1 = S N D 0 + h υ ρ υ γ ( υ + 1 ) k = 0 n ( ( n + 1 k ) υ ( n k ) υ ) ϑ 3 ( t k , z k ) , S N n + 1 = S N 0 + h υ ρ υ γ ( υ + 1 ) k = 0 n ( ( n + 1 k ) υ ( n k ) υ ) ϑ 4 ( t k , z k ) ,
where
ϑ 1 ( t , z ( t ) ) = Φ ( η 1 S N D ( t ) + η 2 S D ( t ) ) S ( t ) ζ ^ S ( t ) , ϑ 2 ( t , z ( t ) ) = λ S D ( t ) + ( η 2 S D ( t ) ) S ( t ) + ζ ^ S N D ( t ) ζ ^ S D ( t ) χ S D ( t ) , ϑ 3 ( t , z ( t ) ) = λ S N D ( t ) + ( η 1 S N D ( t ) ) S ( t ) + ζ ^ S D ( t ) ζ ^ S N D ( t ) , ϑ 4 ( t , z ( t ) ) = χ S D ( t ) ζ ^ S N ( t ) .
The key factors of this computational work were the choice of parameters. According to a report, in Saudi Arabia, around N = 4358 [33] sector-based companies were invested in, and this figure is consistently growing at a rate of n = 0.54%, which acts as the attracting point to future investors. Hence, the enrolment rate for the computational work was Φ = ( 4358 × 0.54 % ) = 102 ,413, and the loss rate was considered as ζ ^ = 1.823 × 10 6 [34], corresponding to the inflation rate of Saudi Arabia. λ = 2.9 365 × 4358 = 1.8 × 10 5 was the sector-affecting rate calculated for one calendar year. Here, values such as S ( 0 ) = 76 , 810 ; S D ( 0 ) = 49 , 160 ; S N D ( 0 ) = 52 ,230; and S N ( 0 ) = 12 ,290 were deliberated for computational convenience. We let η 1 = 1.9 × 10 6 be the rate of transfer of sectors that were affected but not downturned because of some advantages that work towards it. Simultaneously, η 2 = 2 × 10 7 was the rate of sector transfer for sectors that were affected by the pandemic and were downturned. χ = 0.032 [35] represents the recovery rate of the sectors after the pandemic in Saudi Arabia.

Computational Results & Discussion

The fractional investment model was solved using numerical computation, and results are presented in the form of graphs from Figure 2, Figure 3, Figure 4 and Figure 5 for various generalized cases.
Pandemic situations such as COVID-19 and SARS CoV-1 in leading countries that are already in the growing face of attracting investments in various sectors will be crucial, especially around the first month of the impact. Figure 2 clearly depicts the fact that, after or at the end of first month of any pandemic impact, the rate of sectors attracting investments is reduced from the peaking trend, and it becomes saturated thereafter.
Behaviour sectors based on tourism, real estate, manufacturing, stock exchange, and self-employerment are classified as downturned sectors, ( S D ) , during the pandemic, and the impact they experienced is showcased in Figure 3. Interestingly, a hike in the graph can be noted for the first 3 weeks duration, which may be due to factors such as precautionary measures towards lockdown fear. Later, the impact of the pandemic could hit these kinds of sectors and downturn them within three to four months. Lockdown plays a key role in these sectors, which mostly rely on man power and transportation.
Figure 4 clarifies the fortunate sectors that were not downturned even after the pandemic hit, such as medical and pharmaceutical-based sectors. These could initially face the direct impact of the pandemic through such things as accommodation issues, a lack of medicine, doctor sufficiency, and the need for early-stage treatment strategies. Such scenarios were clearly noted in Figure 4 with a slight dip in the curve. However, once this early chaotic situation was settled, these sectors seemed to be continuously growing better than other sectors.
Sectors that belong to the not-affected category are quick learners in finding alternatives, and these sectors really benefited due to the pandemic. Figure 5 portrays the growing trend of not-affected sectors such as IT, social media, and online forums. It can be clearly noted that, without any flaws, these sectors attained growth from even day one of the pandemic, and they exponentially gained along with the pandemic. This growing phase could last for the first six months of the pandemic. In the second half, their growth looks vital.

7. Conclusions

We investigated a system of four nonlinear fractional order equations in the generalized Liouville–Caputo sense to analyze the significance of various sectors in future pandemics and to identify a good era for executing future strategic working plans. The fixed-point theorems of Schauder and Banach, respectively, were used to demonstrate the existence and uniqueness of solutions to the suggested investment model in a pandemic situation. Ulam–Hyers and generalized Ulam–Hyers were established for stability analysis. The fractional form of the model under discussion via the generalized Liouville–Caputo fractional operator has been numerically simulated using the fractional Euler method, a first-order convergent numerical methodology. Each variable’s illustration and dynamical perspective with distinct fractional order values were investigated. We have highlighted the topic’s asymmetries in the remarks. This study based on a fractional investment model was analysed in view of exploring the behaviours of various sectors in a future pandemic situation and suggests a favourable period for executing strategic working plans for the future. An effective Euler’s computing technique was employed through MATLAB. Based on the data available along with some logical assumptions for computational convenience, the results were obtained in the form of graphs, and their factors were analysed to derive these conclusive notes. This model was performed for the generalized Liouville–Caputo type. The classical Caputo type has been utilized for many physical interpretations in mathematical modeling problems. In light of this, to defend our model, comparison graphs (Figure 2, Figure 3, Figure 4 and Figure 5) have been included in the classical Caputo sense ( υ = 0.95 , ρ = 1 ) and the generalized Liouville–Caputo sense υ = 0.95 , ρ = 0.92 , 0.93 , 0.94 , 0.95 ) for S ( t ) , S D ( t ) , S N D ( t ) , and S N ( t ) , which exhibit excellent correlation and give confidence to proceed with this model. During future pandemic situations:
  • Any investment plans are suggested to be put on hold for the first few weeks immediately after the pandemic.
  • Investments in unaffected category sectors such as IT, social media, and online forums look healthier throughout the year.
  • Investors in key sectors that are not downturned, such as the pharmaceutical and medical sectors, are suggested to have precautionary plans for the early dip to explore back-end benefits.
  • Both the above-mentioned sectors are recommendable for investments even during the pandemic.
  • Investors have to wait for the saturation period for sectors such as tourism, real estate, and self-employment.

Author Contributions

Conceptualization, M.S. and P.M.; formal analysis, M.A., M.S. and P.M.; methodology, M.S., P.M. and K.A. All authors have read and agreed to the published version of the manuscript.

Funding

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number INST201.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saudi Arabia Country Report 2022. Available online: https://bti-project.org/en/reports/country-report/SAU (accessed on 9 February 2023).
  2. Investment Sectors. Available online: https://www.investsaudi.sa/en/sectors-opportunities (accessed on 9 February 2023).
  3. Javid, M.; Hasanov, F.J.; Bollino, C.A.; Galeotti, M. Sectoral investment analysis for Saudi Arabia. Appl. Econ. 2022, 54, 4486–4500. [Google Scholar] [CrossRef]
  4. Ali, M.A.S.; Ammer, M.A.; Elshaer, I.A. Determinants of Investment Awareness: A Moder-ating Structural Equation Modeling-Based Model in the Saudi Arabian Context. Mathematics 2022, 10, 3829. [Google Scholar] [CrossRef]
  5. Odibat, Z.; Baleanu, D. Numerical simulation of initial value problems with generalized Caputo-type fractional derivatives. Appl. Numer. Math. 2020, 156, 94–105. [Google Scholar] [CrossRef]
  6. Hohenberg, P.C.; Halperin, B.I. Theory of dynamic critical phenomena. Rev. Mod. Phys. 1977, 49, 435–479. [Google Scholar] [CrossRef]
  7. Datsko, B.; Gafiychuk, V. Complex nonlinear dynamics in subdiffusive activator–inhibitor systems. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 1673–1680. [Google Scholar] [CrossRef]
  8. Datsko, B.; Gafiychuk, V. Complex spatio-temporal solutions in fractional reaction-diffusion systems near a bifurcation point. Fract. Calc. Appl. Anal. 2018, 21, 237–253. [Google Scholar] [CrossRef] [Green Version]
  9. Klafter, J.; Lim, S.; Metzler, R. Fractional Dynamics: Recent Advances; World Scientific: Singapore, 2012. [Google Scholar]
  10. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Amsterdam, The Netherlands, 1998. [Google Scholar]
  11. Valerio, D.; Machado, J.T.; Kiryakova, V. Some pioneers of the applications of frac-tional calculus. Fract. Calc. Appl. Anal. 2014, 17, 552–578. [Google Scholar] [CrossRef] [Green Version]
  12. Machado, J.T.; Kiryakova, V.; Mainardi, F. Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1140–1153. [Google Scholar] [CrossRef] [Green Version]
  13. Kilbas, A.A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science Limited: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  14. Ciegis, R.; Bugajev, A. Numerical approximation of one model of bacterial self-organization. Nonlinear Anal. Model. Control 2012, 17, 253–270. [Google Scholar] [CrossRef] [Green Version]
  15. Zarin, R.; Khaliq, H.; Khan, A.; Ahmed, I.; Humphries, U.W. A Numerical Study Based on Haar Wavelet Collocation Methods of Fractional-Order Antidotal Computer Virus Model. Symmetry 2023, 15, 621. [Google Scholar] [CrossRef]
  16. Alanzi, A.R.; Jamal, F.; Tahir, M.H.; Chesneau, C.; Kanwal, S.; Sami, W. A New Detec-tion Function Model for Distance Sampling Based on the Burr XII Model. Symmetry 2023, 15, 620. [Google Scholar] [CrossRef]
  17. Wang, K.J. The fractal active low-pass filter within the local fractional derivative on the Cantor set. COMPEL-Int. J. Comput. Math. Electr. Electron. Eng. 2023. [Google Scholar] [CrossRef]
  18. Wang, K. Fractal Traveling Wave Solutions For The Fractal-Fractional Ablow-itz–Kaup–Newell–Segur Model. Fractals 2022, 30, 1–9. [Google Scholar] [CrossRef]
  19. He, J.H. Fractal calculus and its geometrical explanation. Results Phys. 2018, 10, 272–276. [Google Scholar] [CrossRef]
  20. Wang, K.J. Diverse wave structures to the modified Benjamin–Bona–Mahony equation in the optical illusions field. Mod. Phys. Lett. B 2023, 2350012. [Google Scholar] [CrossRef]
  21. Feng, Y.Y.; Yang, X.J.; Liu, J.G. On overall behavior of Maxwell mechanical model by the combined Caputo fractional derivative. Chin. J. Phys. 2020, 66, 269–276. [Google Scholar] [CrossRef]
  22. Katugampola, U.N. New approach to a generalized fractional integral. Appl. Math. Comput. 2011, 218, 860–865. [Google Scholar] [CrossRef] [Green Version]
  23. Katugampola, U.N. A new approach to generalized fractional derivatives. Bull. Math. Anal. Appl. 2014, 6, 1–15. [Google Scholar]
  24. Jarad, F.; Abdeljawad, T.; Baleanu, D. On the generalized fractional derivatives and their caputo modification. J. Nonlinear Sci. Appl. 2017, 10, 2607–2619. [Google Scholar] [CrossRef] [Green Version]
  25. Rezapour, S.; Mohammadi, H.; Samei, M.E. SEIR epidemic model for COVID-19 trans-mission by Caputo derivative of fractional order. Adv. Differ. Equ. 2020, 2020, 1–19. [Google Scholar] [CrossRef]
  26. Granas, A.; Dugundji, J. Fixed Point Theory; Springer Science & Business Media: New York, NY, USA, 2013. [Google Scholar]
  27. Ulam, S.M. A Collection of Mathematical Problems; Interscience: New York, NY, USA, 1960; Volume 8. [Google Scholar]
  28. Ulam, S.M. Problems in Modern Mathematics; Courier Corporation: Los Alamos, NM, USA, 2004. [Google Scholar]
  29. Awadalla, M.; Subramanian, M.; Abuasbeh, K. Existence and Ulam–Hyers Stability Results for a System of Coupled Generalized Liouville–Caputo Fractional Langevin Equations with Multipoint Boundary Conditions. Symmetry 2023, 15, 198. [Google Scholar] [CrossRef]
  30. Subramanian, M.; Aljoudi, S. Existence and Ulam–Hyers Stability Analysis for Coupled Differential Equations of Fractional-Order with Nonlocal Generalized Conditions via Generalized Liou-ville–Caputo Derivative. Fractal Fract. 2022, 6, 629. [Google Scholar] [CrossRef]
  31. Subramanian, M.; Manigandan, M.; Zada, A.; Gopal, T.N. Existence and Hyers–Ulam stability of solutions for nonlinear three fractional sequential differential equations with nonlocal boundary conditions. Int. J. Nonlinear Sci. Numer. Simul. 2022. [Google Scholar] [CrossRef]
  32. Ahmad, B.; Aljoudi, S. Investigation of a Coupled System of Hilfer–Hadamard Frac-tional Differential Equations with Nonlocal Coupled Hadamard Fractional Integral Boundary Conditions. Fractal Fract. 2023, 7, 178. [Google Scholar] [CrossRef]
  33. Saudi Arabia Sees 54% Jump in Investment Licenses to 4358 in 2022. Available online: https://www.arabnews.com/node/2246391/business-economy (accessed on 18 February 2023).
  34. Saudi Arabia Inflation Rate. Available online: https://bit.ly/3L5aCf0 (accessed on 15 February 2023).
  35. Saudi Arabia: The Fastest Growing Economy Fighting for Its Future—November 2022. Available online: https://www.mfat.govt.nz/en/trade/mfat-market-reports/market-reports-middle-east/saudi-arabia-the-fastest-growing-economy-fighting-for-its-future-november-2022/ (accessed on 14 February 2023).
Figure 1. The diagram for the proposed model of sector-based investment.
Figure 1. The diagram for the proposed model of sector-based investment.
Symmetry 15 00789 g001
Figure 2. Sectors’ performance during early days of pandemic ( S ) . ( υ = 0.95, ρ = 0.92, 0.93, 0.94, 0.95)—generalized Lioville–Caputo type. ( υ = 0.95 , ρ = 1 ) —Caputo type.
Figure 2. Sectors’ performance during early days of pandemic ( S ) . ( υ = 0.95, ρ = 0.92, 0.93, 0.94, 0.95)—generalized Lioville–Caputo type. ( υ = 0.95 , ρ = 1 ) —Caputo type.
Symmetry 15 00789 g002
Figure 3. Performance of downturned sectors ( S D ) . ( υ = 0.95, ρ = 0.92, 0.93, 0.94, 0.95)—generalized Liouville–Caputo type. ( υ = 0.95, ρ = 1)—Caputo type.
Figure 3. Performance of downturned sectors ( S D ) . ( υ = 0.95, ρ = 0.92, 0.93, 0.94, 0.95)—generalized Liouville–Caputo type. ( υ = 0.95, ρ = 1)—Caputo type.
Symmetry 15 00789 g003
Figure 4. Performance of not-downturned sectors ( S N D ) . ( υ = 0.95, ρ = 0.92, 0.93, 0.94, 0.95)—generalized Lioville–Caputo type. ( υ = 0.95, ρ = 1)—Caputo type.
Figure 4. Performance of not-downturned sectors ( S N D ) . ( υ = 0.95, ρ = 0.92, 0.93, 0.94, 0.95)—generalized Lioville–Caputo type. ( υ = 0.95, ρ = 1)—Caputo type.
Symmetry 15 00789 g004
Figure 5. Performance of unaffected sectors ( S N ) . ( υ = 0.95, ρ = 0.92, 0.93, 0.94, 0.95)—generalized Lioville–Caputo type. ( υ = 0.95, ρ = 1)—Caputo type.
Figure 5. Performance of unaffected sectors ( S N ) . ( υ = 0.95, ρ = 0.92, 0.93, 0.94, 0.95)—generalized Lioville–Caputo type. ( υ = 0.95, ρ = 1)—Caputo type.
Symmetry 15 00789 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Awadalla, M.; Subramanian, M.; Madheshwaran, P.; Abuasbeh, K. Post-Pandemic Sector-Based Investment Model Using Generalized Liouville–Caputo Type. Symmetry 2023, 15, 789. https://doi.org/10.3390/sym15040789

AMA Style

Awadalla M, Subramanian M, Madheshwaran P, Abuasbeh K. Post-Pandemic Sector-Based Investment Model Using Generalized Liouville–Caputo Type. Symmetry. 2023; 15(4):789. https://doi.org/10.3390/sym15040789

Chicago/Turabian Style

Awadalla, Muath, Muthaiah Subramanian, Prakash Madheshwaran, and Kinda Abuasbeh. 2023. "Post-Pandemic Sector-Based Investment Model Using Generalized Liouville–Caputo Type" Symmetry 15, no. 4: 789. https://doi.org/10.3390/sym15040789

APA Style

Awadalla, M., Subramanian, M., Madheshwaran, P., & Abuasbeh, K. (2023). Post-Pandemic Sector-Based Investment Model Using Generalized Liouville–Caputo Type. Symmetry, 15(4), 789. https://doi.org/10.3390/sym15040789

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop