Next Article in Journal
A New Formula for Calculating Uncertainty Distribution of Function of Uncertain Variables
Next Article in Special Issue
Numerical Validation of a Boost Converter Controlled by a Quasi-Sliding Mode Control Technique with Bifurcation Diagrams
Previous Article in Journal
The Multi-Compartment SI(RD) Model with Regime Switching: An Application to COVID-19 Pandemic
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Solutions of the Mathematical Models on the Digestive System and COVID-19 Pandemic by Hermite Wavelet Technique

1
Department of Mathematics, Bangalore University, Bengaluru, Karnataka 560056, India
2
Faculty of Education, Harran University, Sanliurfa 63050, Turkey
3
Departamento de Anatomía y Psicobiología, Universidad de Murcia, 30001 Murcia, Spain
*
Author to whom correspondence should be addressed.
Submission received: 9 November 2021 / Revised: 26 November 2021 / Accepted: 12 December 2021 / Published: 15 December 2021

Abstract

:
This article developed a functional integration matrix via the Hermite wavelets and proposed a novel technique called the Hermite wavelet collocation method (HWM). Here, we studied two models: the coupled system of an ordinary differential equation (ODE) is modeled on the digestive system by considering different parameters such as sleep factor, tension, food rate, death rate, and medicine. Here, we discussed how these parameters influence the digestive system and showed them through figures and tables. Another fractional model is used on the COVID-19 pandemic. This model is defined by a system of fractional-ODEs including five variables, called S (susceptible), E (exposed), I (infected), Q (quarantined), and R (recovered). The proposed wavelet technique investigates these two models. Here, we express the modeled equation in terms of the Hermite wavelets along with the collocation scheme. Then, using the properties of wavelets, we convert the modeled equation into a system of algebraic equations. We use the Newton–Raphson method to solve these nonlinear algebraic equations. The obtained results are compared with numerical solutions and the Runge–Kutta method (R–K method), which is expressed through tables and graphs. The HWM computational time (consumes less time) is better than that of the R–K method.

1. Introduction

Recently, with the help of newly developed computational methods, many experts extracted more deep features of real world problems such as symmetry, optical waves, gravitational potentials and so on. Experts extracted the symmetrical properties of such models by using mathematical norms. Thus, many experts developed some new operators and models. So, the coupled ordinary differential equations (ODEs) of integer and fractional order are used to describe many physical phenomena in symmetry, physics, fluid dynamics, mathematical biology, and bio-modeling. Many experts have considered these models over hundreds of years, and there are many well-developed solution procedures. Often, some models described by ODEs are so complex that a purely analytical solution to such equations is not tractable. For those sorts of complex models, numerical methods are helpful. When the analytical solution is not possible, then switching over to the numerical method is preferred. As a result, we proposed a novel approach called HWM to solve the system of coupled ODEs. The primary purpose is to present and explain a new numerical method for obtaining the approximate solution to the system of coupled equations of integer and fractional orders that cannot be solved exactly.
Wavelet theory is one of the recent emerging approaches in applied mathematics. It has a wide range of applications in the following fields as signal analyses, computer science, mathematical modeling, image processing, and applied sciences. Many mathematicians’ contributions toward wavelets-based numerical methods and fractional models are as follows: the Laguerre wavelets method for the Lane–Emden equation [1], Hermite wavelets method [2,3], Laguerre wavelets collocation method [4,5], numerical solution of some biological models [6,7,8,9,10], and numerical approach for fractional models [11,12,13,14,15,16,17,18,19,20,21,22,23,24], etc. According to the current literature survey, HWM is used to solve ODE, PDE, fractional ODE, and fractional PDEs. Still, mathematical models on the digestive system and fractional models on the COVID-19 pandemic have not been studied by any mathematicians using the wavelet technique. Using the information with the vaccination, experts may observe more profound properties of the vaccination by using mathematical tools that are guaranteed scientifically. This impels us to solve such equations via the Hermite wavelets method (HWM). The interest of the present work is to solve the system of ODEs for both integer and fractional order by using the HWM, and the obtained results are compared with the results of the numerical method.
A variety of exact, approximate and purely numerical methods are available to solve the system of ODEs. Most of these methods are computationally rigorous because they are trial-and-error in character or need complex computations. Some problems are often solved by the system of ordinary differential equations in reality, such as in modern physics, engineering, etc. Anticipating the exact solution for such equations is rigorous. So, it is necessary to use a numerical method to solve such a problem. This paper mainly used the R–K method and HWM to solve the integer and fractional order ODE system. Consider that the mathematical model on the digestive system is of the form [7]:
d T ( x ) d x = α F ( x ) β T ( 0 ) = A d F ( x ) d x = δ T ( x ) + γ M ( x ) F ( 0 ) = B d M ( x ) d x = δ T ( x ) + β M ( x ) d M ( 0 ) = C }
where α represents the quantity of food, and the minus sign indicates the smaller amount of food. β shows the sleep factor, δ is used for the tension rate, γ denotes the recovery term, and d shows the death rate. Moreover, T ( 0 ) = A , F ( 0 ) = B , M ( 0 ) = C are the initial conditions for all compartments of the given model. We considered four different cases to analyze the given model and solved it by the R–K method and HWM. Here, T ( x ) , F ( x ) , M ( x ) stands for tension, food, and medicine, respectively. We study to observe the numerical and symmetrical wave propagations of Equation (1) in graphs and tables. Another important fractional system of ODEs for the COVID-19 pandemic model [11] given as:
D 0 c x δ S = δ μ δ S β δ S ( E + I ) S ( 0 ) = 153 D 0 c x δ E = β δ S ( E + I ) π δ E ( μ δ + γ δ ) E E ( 0 ) = 55 D 0 c x δ I = π δ E σ δ I μ δ I I ( 0 ) = 79 D 0 c x δ Q = γ δ E + σ δ I ( θ δ + μ δ ) Q Q ( 0 ) = 68 D 0 c x δ R = θ δ Q μ δ R R ( 0 ) = 20 }
where 0 < δ < 1, and D 0 c x δ denotes the fractional derivative in the Caputo sense. S , E , I , Q , R represent the susceptible, exposed, infected, isolated, and recovered population, respectively. Λ = μN is recruitment rate, β is the rate at which those susceptible become infected and exposed, π is the rate at which the exposed population become infected, γ is the rate at which exposed people become isolated, σ is the rate at which infected people are added to isolated individuals, θ is the rate at which isolated persons become recovered, μ is the natural death rate and disease-related death rate. We solved this problem for different fractional-order and compared the HWM solution with the R–K method solution.
The R–K method is a well-known and universally accepted method. The main advantages of Runge–Kutta methods are that they are easy to implement, and they are very stable. The Hermite wavelet method is a new method, and here, we showed that it can be used to solve the biological models for both integer and fractional ordered coupled ODEs. The primary disadvantages of Runge–Kutta methods are that they require significantly more computer time than multi-step methods of comparable accuracy, and they do not easily yield good global estimates of the truncation error. However, this cannot happen with the Hermite wavelet method. For the particular domain, we get the solution with less computation time. As a result, the Hermite wavelet method is better than the R–K method. Therefore, we used these two methods in this article.
This paper’s organization is as follows: Section 2 is devoted to the progress of the Hermite wavelet operational matrix of integration, function approximation, and theorems on convergence analysis. Section 3 reveals the method of solution. In Section 4, we present the application of the HWM to the governing model. By using projected method, we observe the symmetrical and wave features of the governing models. In Section 5, we present some important results and discussion. Finally, this paper is completed by giving critical new findings in conclusion in Section 6.

2. Hermite Wavelet Operation Matrix of Integration, the Definition of Fractional Derivatives, and Some Results on Convergence Analysis

Definition 1.
The Riemann–Liouville’s fractional integral of f   ϵ   C μ of the order δ 0 is defined as [15]
J s δ f ( s ) = { f ( s ) i f   δ = 0 1 Γ ( δ ) 0 s ( s t ) δ 1 f ( t ) d t i f   δ > 0 .
Here, Γ denotes the gamma function.
Definition 2.
The Caputo fractional derivative of f ( s )   ϵ   C μ is defined as [15]:
δ f ( s ) s δ = 1 Γ ( m δ ) 0 s ( s t ) m δ 1 f ( m ) ( t )   d t
for m 1 < δ m , m is any positive integer, s > 0 ,   f ( s ) ϵ C μ m ,   μ 1 .
The Hermite wavelet is one of the continuous polynomial basis wavelets studied in [4]. Now, we approximate the solution y ( x ) of the nonlinear differential equation under the Hermite wavelet space as follows:
y ( x ) = n = 1 m = 0 C n , m φ n , m ( x ) ,
where φ n , m ( x ) is given in [4]. We approximate y ( x ) by truncating the series as follows
y ( x ) n = 1 2 k 1 m = 0 M 1 C n , m φ n , m ( x ) = A T φ ( x ) ,
where A and φ ( x ) are 2 k 1 M × 1 matrix,
A T = [ C 1 , 0 , C 1 , M 1 , C 2 , 0 C 2 , M 1 , C 2 k 1 , 0 , C 2 k 1 , M 1 ] ,
φ ( x ) = [ φ 1 , 0 , φ 1 , M 1 , φ 2 , 0 , φ 2 , M 1 , φ 2 k 1 , 0 , φ 2 k 1 , M 1 ] T .
Let { ϕ i , j } be the sequence of Hermite wavelets, n = 1, 2, … and m = 0, 1, … For every fixed n, there is a Hermite space spanned by the elements of the sequence { ϕ i , j } . That is, L ( { ϕ i , j } ) = H 2 [ 0 , 1 ) is Banach space.
Theorem 1
([17]). Let { ϕ i , j ( t ) } be the sequence of Hermite wavelets in C ( [ a ,   b ] ) in t on [ a ,   b ] converges to the function ϕ ( t ) in C ( [ a ,   b ] ) uniformly in t on [ a ,   b ] . Then, ϕ ( t ) is continuous in C ( [ a ,   b ] ) in t on [ a ,   b ] .
Theorem 2
([19]). A continuous function y ( x ) in H 2 [ 0 , 1 ) defined on [ 0 , 1 ) be bounded; then, the Hermite wavelets expansion of y ( x ) converges to it.
Theorem 3
([4]). Suppose that y ( x ) = C m [ 0 , 1 ] and C T ϕ ( x ) is the approximate solutionusing Hermite wavelets. Then, the error bound would be given by
E ( x ) 2 m ! 4 m 2 m ( k 1 ) | y m ( x ) | x [ 0 , 1 ] m a x .
Theorem 4
([16]). Let ( S ( 0 ) , E ( 0 ) , I ( 0 ) , Q ( 0 ) ) be any initial data belonging to R + 4 and ( S ( x ) , E ( x ) , I ( x ) , Q ( x ) ) be the solution corresponding to the initial data. Then, the set R + 4 is a positively invariant set of the model. Furthermore, we have
lim sup x S ( x ) S = δ μ δ ,
lim sup x E ( x ) E = δ π δ + μ δ + γ δ ,
lim sup x I ( x ) I = π δ E σ δ + μ δ ,
lim sup x Q ( x ) Q = γ δ E + σ δ I μ δ + α δ ,
lim sup x R ( x ) θ δ μ δ .
Theorem 5
([19]). Suppose y C p [ 0 , 1 ) is an p times continuously differentiable function such that y = n = 1 2 k 1 y n ( x ) and { ϕ n , m } be a sequence of Hermite wavelets, where n = 1 , , 2 k 1 and m = 0 , , M 1 , k is any positive integer. Let Y n = L ( { φ n . m } ) be the linear space spanned by { ϕ n , m } . If C n T H n ( x ) is the best approximation to y n , then C T H ( x ) approximates y with the following error bound:
| | y C T H ( x ) | | 2 k ( 2 p + 1 ) 2 ( k 1 ) ( p + 1 2 ) ,   K = max y p ( ξ ) , ξ [ n 1 2 k 1 , n 2 k 1 ) .

Operational Matrix of Integration

Hermite wavelet: Wavelets constitute a family of functions constructed from the dilation and translation of a single function called the mother wavelet. When the dilation parameter a and translation parameter b vary continuously, we have the following family of continuous wavelets:
φ a , b ( x ) = | a | 1 / 2 φ ( x b a ) , a , b R , a 0 .
If we restrict the parameters a and b to discrete values as a = a 0 k , b = n b 0 a 0 k , a 0 > 1 , b 0 > 0 , we have the following family of discrete wavelets: φ k , n ( x ) = | a | 1 / 2 φ ( a 0 k x n b 0 ) , a , b R , a 0 , where φ k , n form a wavelet basis for L 2 ( R ) . In particular, when a 0 = 2 and b 0 = 1 , then φ k , n ( x ) forms an orthonormal basis. Hermite wavelets are defined as,
φ n , m ( x ) = { 2 k + 1 2 π h m ( 2 k x 2 n + 1 ) , n 1 2 k 1 x < n 2 k 1 , 0 , otherwise .
where m = 0 , 1 , , M 1 and n = 1 , 2 , , 2 k 1 , k is assumed as any positive integer. Here, h m ( x ) are the Hermite polynomials of degree m concerning the weight function W ( x ) = 1 x 2 on the real line R and satisfy the following recurrence formula,
h 0 ( x ) = 1 , h 1 ( x ) = 2 x , h m + 2 ( x ) = 2 x h m + 1 ( x ) 2 ( m + 1 ) h m ( x ) ,
where m = 0 , 1 , 2 ,
Consider six Hermite wavelet bases at k = 1 and M = 6 = N; then, we get
φ 1 , 0 ( x ) = 2 π φ 1 , 1 ( x ) = 1 π ( 8 x 4 ) φ 1 , 2 ( x ) = 1 π ( 32 x 2 32 x + 4 ) φ 1 , 3 ( x ) = 1 π ( 128 x 3 192 x 2 + 48 x + 8 ) φ 1 , 4 ( x ) = 1 π ( 512 x 4 1024 x 3 + 384 x 2 + 128 x 40 ) φ 1 , 5 ( x ) = 1 π ( 2048 x 5 5120 x 2 + 2560 x 3 + 1280 x 2 800 x + 16 )
where
φ ( x ) = [ φ 1 , 0 ( x ) , φ 1 , 1 ( x ) , φ 1 , 2 ( x ) φ 1 , 3 ( x ) φ 1 , 4 ( x ) φ 1 , 5 ( x ) ] .
Now, integrate the above basis concerning x limit from 0 to x; then, express as a linear combination of Hermite wavelet basis as
0 x ϕ 1 , 0 ( x ) = [ 1 2 1 4 0 0 0 0 ] ϕ 6 ( x ) 0 x ϕ 1 , 1 ( x ) = [ 1 4 0 1 8 0 0 0 ] ϕ 6 ( x ) 0 x ϕ 1 , 2 ( x ) = [ 1 3 0 0 1 12 0 0 ] ϕ 6 ( x ) 0 x ϕ 1 , 3 ( x ) = [ 5 4 0 0 0 1 16 0 ] ϕ 6 ( x ) 0 x ϕ 1 , 4 ( x ) = [ 2 5 0 0 0 0 1 12 ] ϕ 6 ( x ) 0 x ϕ 1 , 5 ( x ) = [ 23 3 0 0 0 0 0 ] ϕ 6 ( x ) + 1 24 ϕ 1 , 6 ¯ ( x ) .
Therefore,
0 x ϕ ( x ) d x = P 6 × 6 ϕ 6 ( x ) + ϕ 6 ¯ ( x )
where
P 6 × 6 = [ 1 2 1 4 0 0 0 0 1 4 0 1 8 0 0 0 1 3 0 0 1 12 0 0 5 4 0 0 0 1 16 0 2 5 0 0 0 0 1 20 23 3 0 0 0 0 0 ] , φ 6 ¯ ( x ) = [ 0 0 0 0 0 1 24 φ 1 , 6 ¯ ( x ) ] .

3. Hermite Wavelets Method

Here, we would like to express the solution of a system of ODE in terms of Hermite wavelets.

3.1. Hermite Wavelets Method for Digestive Model

Assume that
T = A T ϕ ( x ) .
Integrate Equation (5) concerning x, limit from 0 to x, we get
T = T ( 0 ) + A T [ P ϕ ( x ) + ϕ ( x ) ] ¯
T = 0.9999 + A T [ P ϕ ( x ) + ϕ ( x ) ] . ¯
Again, assume
F = B T ϕ ( x ) .
Integrate Equation (7) concerning x, limit from 0 to x, we get
F = F ( 0 ) + B T [ P ϕ ( x ) + ϕ ( x ) ] ¯
F = 0.0001 + B T [ P ϕ ( x ) + ϕ ( x ) ] ¯ .
Again, choose
M = C T ϕ ( x ) .
Integrate Equation (9) concerning x, limit from 0 to x.
M = M ( 0 ) + C T [ P ϕ ( x ) + ϕ ( x ) ] ¯
M = C T [ P ϕ ( x ) + ϕ ( x ) ] ¯
Now, substitute (5) to (10) in (1); then, we get
A T ϕ ( x ) = α [ 0.001 + B T [ P ϕ ( x ) + ϕ ( x ) ] ¯ ] β
B T ϕ ( x ) = δ [ 0.9999 + A T [ P ϕ ( x ) + ϕ ( x ) ] ¯ ] + γ C T [ P ϕ ( x ) + ϕ ( x ) ] ¯
C T ϕ ( x ) = δ [ 0.9999 + A T [ P ϕ ( x ) + ϕ ( x ) ] ¯ ] + β [ C T [ P ϕ ( x ) + ϕ ( x ) ] ¯ ] d .
Collocate the above three equations by the following collocation points, x i = 2 i 1 2 N ,   i = 1 ,   2 ,   3 ,   ,   N , where N represents the size of the matrix. Then, we get a system containing 3N algebraic equations and solve these equations by the Newtons–Raphson method that yields Hermite wavelet unknown coefficient values. Then, substitute these unknown coefficient values in (6), (8), and (10); we get the corresponding Hermite wavelet-based numerical solutions.

3.2. Hermite Wavelets Method for Fractional COVID-19 Pandemic Model

Assume that
S = A T ϕ ( x ) ,   E = B T ϕ ( x ) ,   I = C T ϕ ( x ) ,   Q = D T ϕ ( x ) ,   R = F T ϕ ( x ) .
Integrate the equations in (11) concerning x, limit from 0 to x, along with the initial conditions in (2).
S = 153 + A T [ P φ ( x ) + φ ( x ) ¯ ] E = 55 + B T [ P φ ( x ) + φ ( x ) ¯ ] I = 79 + C T [ P φ ( x ) + φ ( x ) ¯ ] Q = 68 + D T [ P φ ( x ) + φ ( x ) ¯ ] R = 20 + F T [ P φ ( x ) + φ ( x ) ¯ ] }
Now, differentiate (12) fractionally concerning x using the Caputo derivative definition; then, we get
D 0 c x δ S = D 0 c x δ ( 153 + A T [ P φ ( x ) + φ ( x ) ¯ ] ) D 0 c x δ E = D 0 c x δ ( 55 + B T [ P φ ( x ) + φ ( x ) ¯ ] ) D 0 c x δ I = D 0 c x δ ( 79 + C T [ P φ ( x ) + φ ( x ) ¯ ] ) D 0 c x δ Q = D 0 c x δ ( 68 + D T [ P φ ( x ) + φ ( x ) ¯ ] ) D 0 c x δ R = D 0 c x δ ( 20 + F T [ P φ ( x ) + φ ( x ) ¯ ] ) } .
Substitute (13) and (12) in (2) and collocate the obtained equations by the following collocation points, x i = 2 i 1 2 N ,   i = 1 ,   2 ,   3 ,   ,   N , where N represents the size of the matrix. Then, we get a system containing 5N nonlinear algebraic equations and solve these equations by the Newtons–Raphson method that yields the Hermite wavelet unknown coefficients’ values. Then, substitute these unknown coefficient values in (12), and we get the corresponding Hermite wavelet-based fractional numerical solutions.

4. Applications of the Proposed Method

This section discusses the given system by assuming different values for the sleep factor, tension, food rate, death rate, and medicine term. We solve this system using the Hermite wavelet method, and the results are compared with the numerical method. Table 1, Table 2 and Table 3 show the different values of the T ( x ) , F ( x ) , and M ( x ) in the interval [0,2] in which fixed initial values have been taken for all three cases. The numerical calculations have been achieved by using the HWM.
Example 1.
Here, we solved Equation (1) by the Hermite wavelet method. Figure 1, Figure 2 and Figure 3 represent the tension, food, and medicine on varying food intake. Table 1 justifies the numerical values by HWM with the R–K method at the constant values as follows: α = 0.1 ,   β = 0.2 , δ = 0.3 , γ = 0.4 , d = 0.5 , with initial conditions T ( 0 ) = 0.9999 ,   F ( 0 ) = 0.0001 , M ( 0 ) = 0 . In model (1), the negative sign indicates a smaller amount. Figure 1 shows that too little food consumption affects the nervous digestive system that leads to increased tension. Figure 2 shows that continuing the consumption process of too little food leads to other complications to the body; this automatically increases medicine consumption. Figure 3 shows the effect of α on food. Table 2 justifies the numerical values by HWM with the R–K method at the constant values as follows, α = 1 , β = 0.2 , δ = 1 , γ = 0.4 , d = 0.1 , with same initial conditions. Tension is directly affected by sleep; when tension increases, it causes a decrease in sleep. So, Figure 4 reflects that the sleep factor is inversely proportional to the tension. Tension increases while the tension factor δ increases; this can be seen in Figure 5. When tension increases automatically, to control the high tension, medicines are preferred. The consumption of medicine is directly proportional to the tension that reflects Figure 6. Table 3 justifies the numerical values by HWM with the R–K method at the constant values as follows: α = 2 , β = 0.2 , δ = 2 , γ = 0.4 , d = 0.01 , with the same initial conditions. Figure 7 shows that increasing the recovery factor leads to decreasing the consumption of medicine. Any disease that causes more death for such disease consumption of medicine is also increased. The death rate and medicine consumption due to the nervous digestive system increase simultaneously, as shown in Figure 8.
Example 2.
Consider a fractional ordered COVID-19 pandemic model given in the Equation (2). Here, we solved this model by using HWM with following parameter values are as follows, Λ = 0.145; μ = 0.000411; β = 0.00038; π = 0.00211; γ = 0.0021; σ =0.0169; θ = 0.0181, with total population N = 375. Obtained results are compared with R–K method through graphs (9–13). Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13 reflects the HWM solution at δ = 1 ,   0.9 ,   0.8 ,   0.7 for the Susceptible, Exposed, Infected, Isolated, and Recovered population.

5. Results and Discussion

From Figure 1, Figure 3, Figure 5 and Figure 7, it can be observed that as the tension rate increases simultaneously, Medication also increases, which leads to the stomach’s upset. Consumption of too little (too much) food also impacts the nervous-digestive system. The death rate is decreasing by increasing the consumption of Medicine. Also, from Figure 1, too little food consumption causes the high secretion of digestive juice that disturbs the stomach. The mathematical Equation (2) represents the COVID-19 pandemic model with isolated classes in fractional order. A system of fractional-order differential equations describes this model includes the following factors, susceptible, exposed, infected, isolated, and recovered. The construction and analysis of the model (2) help us understand the mechanism of the transmission and the characteristics of diseases. Therefore, we can propose effective strategies to predict, prevent, and restrain diseases (COVID-19) and protect population health. In this model, the exposed and infected population is linked to the susceptible population. Also, we assume that the natural death rate includes the disease death rate. If there is no disease symptom, the exposed class moves with a specific rate to the isolated class, but it moves to the infected class when symptoms are developed. In Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13, we plot numerical solutions of the model (2) obtained by the HWM and the R–K method when N = 375 for α = 1, 0.9, 0.8, 0.7. From these figures, the results obtained using the HWM algorithm match the results of the R–K method well, which implies that the presented method can predict the behavior of these variables accurately in the region under consideration. The increasing number of days susceptible populations decrease, as shown in Figure 9, because the susceptible population gets exposed and infected. Figure 13 shows that when a susceptible population gets more infected, the death rate increases linearly. Figure 12 shows that when a susceptible population gets more infected, the home quarantine of the population also increases. The susceptible population decreases on increasing exposed population it can be seen in Figure 11.

6. Conclusions

This paper has successfully applied the HWM to the digestive model, which is assumed for sleep, tension, food, and death rates. We have also observed the results numerically in detail. Analysis of this model has been shown through figures. Moreover, we have introduced the numerical values in the tables. From the tables, it may be concluded that approximate numerical results are very close to the numerical solutions (R–K method) of the governing model. It is also used to comment that these results are very close the symmetrical wave simulations of the governing model of Equation (1). Thus, considered problems are introduced to test the proposed method’s efficiency, accuracy, and validity in symmetrical aspect. Also, Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 reveals that how these parameters impact the digestive system.
Secondly, we observed one another model described by a fractional-order differential equation system with different terms: susceptible, exposed, infected, isolated, and recovered in the population. We have applied the HWM and Runge–Kutta method to simulate the proposed model. Moreover, it can also be presented that this method may also be applied for obtaining numerical solutions of the other differential equation arising in health problems in extracting features of the governing models. The method used in this paper helps to predict, prevent, and restrain diseases (COVID-19) and protect population health. Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13 reveal how the susceptible, exposed, infected, isolated, and recovered populations in symmetrical sense are interrelated.

Author Contributions

Conceptualization, Y.G.S.; formal analysis, H.M.B.; investigation and writing—original draft preparation, K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This paper has been partially founded by Ministerio de Ciencia, Innovación y Universidades (No. PGC2018-097198-B-I00) and Fundación Séneca de la Región de Murcia (No. 20783/PI/18).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data is available from the corresponding author with a reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yousefi, S.A. Legendre wavelet method for solving differential equations of Lane–Emden type. Appl. Math. Comput. 2006, 181, 1417–1422. [Google Scholar] [CrossRef]
  2. Verma, A.K.; Rawani, M.K.; Cattani, C. A numerical scheme for a class of generalized Burgers’ equation based on Haar wavelet nonstandard finite difference method. Appl. Numer. Math. 2021, 168, 41–54. [Google Scholar] [CrossRef]
  3. Hosseininia, M.; Heydari, M.H.; Cattani, C. A wavelet method for nonlinear variable-order time fractional 2D Schrödinger equation. Discret. Contin. Dynam. Systems-S 2021, 14, 2273–2295. [Google Scholar] [CrossRef]
  4. Shiralashetti, S.C.; Kumbinarasaiah, S. Theoretical study on continuous polynomial wavelet bases through wavelet 288 series collocation method for nonlinear lane-Emden type equations. Appl. Math. Comput. 2017, 315, 591–602. [Google Scholar]
  5. Shiralashetti, S.C.; Kumbinarasaiah, S. Laguerre wavelets collocation method for the numerical solution of the Benjamina–Bona–Mohany equations. J. Taibah Univ. Sci. 2019, 13, 9–15. [Google Scholar] [CrossRef] [Green Version]
  6. Agarwal, P.; Singh, R.; Rehman, A. Numerical solution of a hybrid mathematical model of dengue transmission with relapse and memory via Adam–Bashforth–Moulton predictor-corrector scheme. Chaos Solitons Fractals 2021, 143, 110564. [Google Scholar] [CrossRef]
  7. Sánchez, Y.J.; Sabir, Z.; Günerhan, H.; Baskonus, H.M. Analytical and approximate solutions of design of a novel biological nervous stomach mathematical model. Discret. Dynam. Nat. Soc. 2020, 2020, 9. [Google Scholar]
  8. Praveen, A.; Ram, S. Modelling of transmission dynamics of Nipah virus (Niv): A fractional-order Approach. Phys. A Stat. Mech. Appl. 2020, 547, 124243. [Google Scholar]
  9. Günay, B.; Agarwal, P.; Guirao, J.L.G.; Momani, S. A Fractional Approach to a Computational Eco-Epidemiological Model with Holling Type-II Functional Response. Symmetry 2021, 13, 1159. [Google Scholar] [CrossRef]
  10. Attiqul, R.; Ram, S.; Praveen, A. Modeling, analysis, and prediction of new variants of COVID-19 and dengueco-infection on complex network. Chaos Solitons Fractals 2021, 150, 111008. [Google Scholar]
  11. Goyal, M.; Baskonus, H.M.; Prakash, A. Regarding new positive, bounded, and convergent numerical solution of nonlinear time-fractional HIV/AIDS transmission model. Chaos Solitons Fractals 2020, 139, 110096. [Google Scholar] [CrossRef]
  12. Yel, G.; Akturk, T. A New Approach to (3+1) Dimensional Boiti–Leon–Manna–Pempinelli Equation. Appl. Math. Nonlinear Sci. 2020, 5, 309–316. [Google Scholar] [CrossRef]
  13. Gao, W.; Veeresha, P.; Prakasha, D.G.; Baskonus, H.M.; Yel, G. New approach for the model describing the deathly disease in pregnant women using Mittag Leffler function. Chaos Solitons Fractals 2020, 134, 109696. [Google Scholar] [CrossRef]
  14. Singh, J.; Kumar, D.; Hammouch, Z.; Atangana, A. A fractional epidemiological model for computer viruses pertaining to a new fractional derivative. Appl. Math. Comput. 2018, 316, 504–515. [Google Scholar] [CrossRef]
  15. Kumbinarasaiah, S.; Rezazadeh, H. Numerical solution for the fractional-order one-dimensional telegraph equation via wavelet technique. Int. J. Nonlinear Sci. Numer. Simul. 2020, 22, 767–780. [Google Scholar]
  16. Zhang, Z.; Zeb, A.; Egbelowo, O.F.; Vedat, S.E. Dynamics of a fractional-order mathematical model for COVID-19 epidemic. Adv. Differ. Equ. 2020, 420, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Vladimirov, V.S. Equations of Mathematical Physics; Marcel Dekker Incorporated: New York, NY, USA, 1971. [Google Scholar]
  18. Gupta, A.K.; Saha Ray, S. An investigation with Hermite Wavelets for accurate solution of Fractional Jaulent–Miodek equation associated with energy-dependent Schrödinger potential. Appl. Math. Comput. 2015, 270, 458–471. [Google Scholar] [CrossRef]
  19. Shiralashetti, S.C.; Kumbinarasaiah, S. New generalized operational matrix of integration to solve nonlinear singular boundary value problems using Hermite wavelets. Arab J. Basic Appl. Sci. 2019, 26, 2591–2600. [Google Scholar] [CrossRef] [Green Version]
  20. Lanbaran, N.M.; Celik, E.; Yigider, M. Evaluation of Investment Opportunities with Interval-Valued Fuzzy Topsis Method. Appl. Math. Nonlinear Sci. 2020, 5, 461–474. [Google Scholar] [CrossRef]
  21. Ihan, E.; Kiymaz, I.O. A generalization of truncated M-fractional derivative and applications to fractional differential equations. Appl. Math. Nonlinear Sci. 2020, 5, 171–188. [Google Scholar]
  22. Durur, H.; Yokus, A. Exact solutions of (2+1)-Ablowitz-Kaup-Newell-Segur equation. Appl. Math. Nonlinear Sci. 2021. [CrossRef]
  23. Li, Y.M.; Rashid, S.; Hammouch, Z.; Baleanu, D.; Chu, Y.M. New Newton’s Type Estimates Pertaining to Local Fractional Integral via Generalized p-Convexity with Applications. Fractals 2021, 29, 2140018-747. [Google Scholar] [CrossRef]
  24. Aksoy, N.Y. The Solvability of First Type Boundary Value Problem for a Schrödinger Equation. Appl. Math. Nonlinear Sci. 2020, 5, 211–220. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Influence of α on tension when β = 0.2 , δ = 0.3 , γ = 0.4 , d = 0.5 .
Figure 1. Influence of α on tension when β = 0.2 , δ = 0.3 , γ = 0.4 , d = 0.5 .
Symmetry 13 02428 g001
Figure 2. Influence of α on medicine when β = 0.2 , δ = 0.3 , γ = 0.4 , d = 0.5 .
Figure 2. Influence of α on medicine when β = 0.2 , δ = 0.3 , γ = 0.4 , d = 0.5 .
Symmetry 13 02428 g002
Figure 3. Influence of α on food when β = 0.2 , δ = 0.3 , γ = 0.4 , d = 0.5 .
Figure 3. Influence of α on food when β = 0.2 , δ = 0.3 , γ = 0.4 , d = 0.5 .
Symmetry 13 02428 g003
Figure 4. Influence of β on tension when α = 0.1 , δ = 0.3 , γ = 0.4 , d = 0.5 .
Figure 4. Influence of β on tension when α = 0.1 , δ = 0.3 , γ = 0.4 , d = 0.5 .
Symmetry 13 02428 g004
Figure 5. Influence of δ on tension when α = 0.1 , β = 0.2 , γ = 0.4 , d = 0.5 .
Figure 5. Influence of δ on tension when α = 0.1 , β = 0.2 , γ = 0.4 , d = 0.5 .
Symmetry 13 02428 g005
Figure 6. Influence of δ on medicine when α = 0.1 , β = 0.2 , γ = 0.4 , d = 0.5 .
Figure 6. Influence of δ on medicine when α = 0.1 , β = 0.2 , γ = 0.4 , d = 0.5 .
Symmetry 13 02428 g006
Figure 7. Influence of γ on medicine when α = 0.1 , β = 0.2 ,   δ = 0.2 , d = 0.5 .
Figure 7. Influence of γ on medicine when α = 0.1 , β = 0.2 ,   δ = 0.2 , d = 0.5 .
Symmetry 13 02428 g007
Figure 8. Influence of d on medicine when α = 0.1 , β = 0.2 ,   δ = 0.2 , γ = 0.4 .
Figure 8. Influence of d on medicine when α = 0.1 , β = 0.2 ,   δ = 0.2 , γ = 0.4 .
Symmetry 13 02428 g008
Figure 9. Graphical representation of the susceptible population at different values of δ by HWM and the R–K method.
Figure 9. Graphical representation of the susceptible population at different values of δ by HWM and the R–K method.
Symmetry 13 02428 g009
Figure 10. Graphical representation of infected population at different values of δ by HWM and the R–K method.
Figure 10. Graphical representation of infected population at different values of δ by HWM and the R–K method.
Symmetry 13 02428 g010
Figure 11. Graphical representation of the exposed population at different values of δ by HWM and the R–K method.
Figure 11. Graphical representation of the exposed population at different values of δ by HWM and the R–K method.
Symmetry 13 02428 g011
Figure 12. Graphical representation of the isolated population at different values of δ by HWM and the R–K method.
Figure 12. Graphical representation of the isolated population at different values of δ by HWM and the R–K method.
Symmetry 13 02428 g012
Figure 13. Graphical representation of the recovered population at different values of δ by HWM and the R–K method.
Figure 13. Graphical representation of the recovered population at different values of δ by HWM and the R–K method.
Symmetry 13 02428 g013
Table 1. Numerical comparison of HWM at α = 0.1 , β = 0.2 , δ = 0.3 , γ = 0.4 , d = 0.5 with numerical results.
Table 1. Numerical comparison of HWM at α = 0.1 , β = 0.2 , δ = 0.3 , γ = 0.4 , d = 0.5 with numerical results.
xHWMNumerical
Result
HWMNumerical
Result
HWMNumerical
Result
T(x)T(x)F(x)F(x)M(x)M(x)
0.00.99990.99990.00010.00010.00000
0.10.97980.97990.02820.0282−0.0805−0.0804
0.20.95930.95950.05240.0525−0.1620−0.1618
0.30.93870.93900.07280.0729−0.2445−0.2442
0.40.91790.91830.08920.0893−0.3281−0.3277
0.50.89690.89740.10160.1018−0.4127−0.4122
0.60.87590.87650.11000.1103−0.4985−0.4979
0.70.85480.85550.11420.1147−0.5852−0.5845
0.80.83360.83440.11440.1150−0.6732−0.6724
0.90.81250.81340.11040.1112−0.7622−0.7613
1.00.79140.79240.10220.1032−0.8524−0.8514
1.10.77050.77160.08970.0909−0.9438−0.9427
1.20.74960.75080.07290.0743−1.0364−1.0352
1.30.72900.73030.05170.0534−1.1302−1.1289
1.40.70860.71000.02620.0281−1.2254−1.2240
1.50.68850.6900−0.0038−0.0015−1.3218−1.3203
1.60.66870.6703−0.0383−0.0357−1.4195−1.4179
1.70.64930.6510−0.0773−0.0744−1.5187−1.5170
1.80.63030.6321−0.1208−0.1176−1.6193−1.6175
1.90.61170.6136−0.1690−0.1654−1.7213−1.7194
2.00.59370.5957−0.2218−0.2178−1.8248−1.8228
Table 2. Numerical comparison of HWM at α = 1 , β = 0.2 , δ = 1 , γ = 0.4 , d = 0.1 with numerical results.
Table 2. Numerical comparison of HWM at α = 1 , β = 0.2 , δ = 1 , γ = 0.4 , d = 0.1 with numerical results.
xHWMNumerical
Result
HWMNumerical
Result
HWMNumerical
Result
T(x)T(x)F(x)F(x)M(x)M(x)
0.00.99990.99990.00010.00010.00000
0.10.98490.98490.09850.0984−0.0401−0.0401
0.20.97950.97950.19420.1941−0.0808−0.0808
0.30.98360.98370.28820.2882−0.1222−0.1222
0.40.99710.99730.38140.3815−0.1648−0.1646
0.51.01991.02020.47470.4749−0.2088−0.2085
0.61.05211.05250.56900.5693−0.2545−0.2541
0.71.09381.09430.66510.6655−0.3022−0.3017
0.81.14521.14580.76390.7644−0.3523−0.3516
0.91.20671.20750.86620.8670−0.4051−0.4043
1.01.27861.27960.97310.9740−0.4610−0.4600
1.11.36151.36271.08541.0865−0.5204−0.5192
1.21.45591.45731.20411.2054−0.5837−0.5822
1.31.56251.56421.33021.3318−0.6512−0.6495
1.41.68221.68421.46481.4667−0.7236−0.7216
1.51.81581.81811.60911.6112−0.8013−0.7990
1.61.96441.96701.76411.7666−0.8847−0.8822
1.72.12912.13201.93131.9341−0.9746−0.9717
1.82.31122.31452.11192.1150−1.0715−1.0683
1.92.51222.51582.30742.3109−1.1761−1.1725
2.02.73352.73752.51922.5232−1.2890−1.2850
Table 3. Numerical comparison of HWM at α = 2 , β = 0.2 , δ = 2 , γ = 0.4 , d = 0.01 with numerical results.
Table 3. Numerical comparison of HWM at α = 2 , β = 0.2 , δ = 2 , γ = 0.4 , d = 0.01 with numerical results.
xHWMNumerical
Result
HWMNumerical
Result
HWMNumerical Result
T(x)T(x)F(x)F(x)M(x)M(x)
0.00.99990.99990.00010.000100
0.10.96010.96020.19610.1960−0.0308−0.0307
0.20.88230.88250.37920.3790−0.0605−0.0603
0.30.76980.77010.54190.5416−0.0878−0.0875
0.40.62730.62770.67810.6777−0.1118−0.1114
0.50.46070.46120.78240.7819−0.1316−0.1311
0.60.27680.27740.85080.8502−0.1465−0.1459
0.70.08290.08360.88080.8801−0.1560−0.1553
0.8−0.1129−0.11210.87140.8706−0.1597−0.1589
0.9−0.3030−0.30210.82330.8224−0.1576−0.1567
1.0−0.4798−0.47880.73850.7375−0.1499−0.1489
1.1−0.6364−0.63530.62080.6197−0.1370−0.1359
1.2−0.7670−0.76580.47510.4739−0.1195−0.1183
1.3-0.8674−0.86610.30750.3062−0.0982−0.0969
1.4−0.9351−0.93370.12490.1235−0.0741−0.0727
1.5−0.9700−0.9685−0.0652−0.0667−0.0484−0.0469
1.6−0.9751−0.9735−0.2549−0.2565−0.0222−0.0206
1.7−0.9565−0.9548−0.4369−0.43860.00330.0050
1.8−0.9244−0.9226−0.6043−0.60610.02690.0287
1.9−0.8934−0.8915−0.7516−0.75350.04740.0493
2.0−0.8831−0.8811−0.8752−0.87720.06400.0660
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Srinivasa, K.; Baskonus, H.M.; Guerrero Sánchez, Y. Numerical Solutions of the Mathematical Models on the Digestive System and COVID-19 Pandemic by Hermite Wavelet Technique. Symmetry 2021, 13, 2428. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13122428

AMA Style

Srinivasa K, Baskonus HM, Guerrero Sánchez Y. Numerical Solutions of the Mathematical Models on the Digestive System and COVID-19 Pandemic by Hermite Wavelet Technique. Symmetry. 2021; 13(12):2428. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13122428

Chicago/Turabian Style

Srinivasa, Kumbinarasaiah, Haci Mehmet Baskonus, and Yolanda Guerrero Sánchez. 2021. "Numerical Solutions of the Mathematical Models on the Digestive System and COVID-19 Pandemic by Hermite Wavelet Technique" Symmetry 13, no. 12: 2428. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13122428

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