Next Article in Journal
MSCF-Net: Attention-Guided Multi-Scale Context Feature Network for Ship Segmentation in Surveillance Videos
Previous Article in Journal
Towards Discriminability with Distribution Discrepancy Constrains for Multisource Domain Adaptation
Previous Article in Special Issue
Multi-Soliton, Soliton–Cnoidal, and Lump Wave Solutions for the Supersymmetric Boussinesq Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fibonacci Wavelet Collocation Method for Solving Dengue Fever SIR Model

1
Department of Mathematics, Atma Ram Sanatan Dharma College, University of Delhi, New Delhi 110021, India
2
Department of Mathematics, Jamia Millia Islamia, New Delhi 110025, India
3
Department of Mathematics, Zakir Husain Delhi College, University of Delhi, New Delhi 110002, India
*
Author to whom correspondence should be addressed.
Submission received: 13 July 2024 / Revised: 5 August 2024 / Accepted: 11 August 2024 / Published: 20 August 2024

Abstract

:
The main focus in this manuscript is to find a numerical solution of a dengue fever disease model by using the Fibonacci wavelet method. The operational matrix of integration has been obtained using Fibonacci wavelets. The proposed method is called Fibonacci wavelet collocation method (FWCM). This biological model has been transformed into a system of nonlinear algebraic equations by using the Fibonacci wavelet collocation scheme. Afterward, this system has been solved by using the Newton–Raphson method. Finally, we provide evidence that our results are better than those obtained by various current approaches, both numerically and graphically, demonstrating the method’s accuracy and efficiency.

1. Introduction

Mathematics plays an important role in the betterment of human life. To solve real-life problems, a mathematical model has been proved to be one of the powerful tools; for example, one can discuss the severity of a disease in a particular region by using models. Over the years, mathematicians have done a tremendous job of revealing biological models. Usually, various biological models can be expressed as a system of ordinary differential equations. Most of these biological models are complex in nature, and only a few of them have analytical solutions. Consequently, a reasonable approximate solution is required. Numerical approaches are usually the best option for these systems since they lead us to answers and characterize the dynamic behavior of such systems.
One of the most dangerous and common diseases transmitted, due to mosquitoes in the world, is dengue [1,2]. About 390 million new cases of dengue are reported in a year, 96 million of which are symptomatic, according to recent data [3,4]. Asia accounts for the majority of newly diagnosed cases each year (about 70% of 390 million cases), with infection distributions of roughly 16.4%, 13.8%, and 0.2% shared by Africa, the Americas, and Oceania, respectively [3]. Aedes albopictus and Aedes aegypti are the primary vectors that spread dengue [5]. There are four different but closely related serotypes of the Flavivirus genus known as dengue virus (DENV): DENV-1, DENV-2, DENV-3, and DENV-4. Homologous immunity refers to the fact that an individual infected with any one of the types will never become infected with the same type again. However, a person loses immunity to the other varieties within around 12 weeks of contracting any of the aforementioned types.
Dengue fever’s mathematical modeling is [6,7]
d S ( t ) d t = μ h ( 1 S ( t ) ) α S ( t ) R ( t )
d I ( t ) d t = α S ( t ) R ( t ) β I ( t )
d R ( t ) d t = γ ( 1 R ( t ) ) I ( t ) δ R ( t )
where S ( t ) is the susceptible population (potential victims), I ( t ) is the infected population, R ( t ) represents recovered population from dengue at time ( t ) , μ h is the natural death rate of the human population, α is the vaccine efficiency coefficient, δ is the number of deaths among the susceptible population, γ is the recovery rate, and β is the rate of infection in the human population. This model is a variation of the classical SIR (Susceptible–Infected–Recovered) model. Over the years, various numerical techniques have been applied by researchers for the SIR model of dengue fever, and few of them are as follows: M. Khalid et al. [8] proposed a new approach called the Perturbation Iteration Algorithm (PIM), Rangkuti et al. [9] implemented the Homotopy Perturbation Method and Variational Iteration method, Mungkasi [10] recommended an Improved Variational Iteration Method for the SIR model of dengue, Umar et al. [11] established a stochastic numerical computing scheme with artificial neural networks, Lede et al. [12] presented Euler and Heun methods, and Fei et al. [13] suggested a numerical collocation method using Bernoulli wavelets (BWCM).
Wavelet techniques emerged in the mid-1980s, drawing inspiration from concepts in both applied and pure mathematics. Based on an analysis function, Joseph Fourier created a technique in 1807 for expressing signals with a set of coefficients. He established the mathematical foundation on which the wavelet theory was created. Alfred Haar made the initial reference to wavelets in his doctoral dissertation in 1909. Currently, wavelet techniques have emerged as an effective approach for solving differential equations. Wavelet-based approaches are the most popular among scientists and engineers because of their smoothness, good localization, acceptability, and orthonormality. The wavelet theory has been an area of extensive research during the last 20 years. Numerous types of nonlinear differential equations have been solved by researchers using Haar, Legendre, and Chebyshev wavelets [14]. Other wavelets that have gained popularity among researchers include the Taylor wavelets [15], Müntz wavelets [16], Bernoulli wavelets [17], CAS wavelets [18], and Jacobi wavelets [19]. Specifically, wavelet-based collocation techniques have become popular in the field of numerical analysis, basically due to their straightforward methodology, easy computation, and fast convergence. It is important to note that different wavelet families, such as the Haar, Legendre, Chebyshev, Bernoulli, Euler, Hermite, Gegenbauer, and Fibonacci wavelets, have been used to solve a wide range of linear and nonlinear differential equations.
In the light of the interesting features and numerous advantages of wavelet-based numerical techniques over other currently used numerical methods, we are highly motivated to solve the Dengue Fever SIR model using the Fibonacci wavelet collocation method (FWCM). As far as we are aware, this wavelet technique has not yet been implemented for obtaining numerical solutions in the case of the dengue model. We make use of this to provide a novel numerical approach to discuss the above model based on compactly supported Fibonacci wavelets.
An overview of this article is as follows: Section 2 includes some basic definitions used throughout this article. Section 3 contains details about function approximation by Fibonacci wavelets. Section 4 is devoted to the construction of the operational matrix of integration based on Fibonacci wavelets for different values of M. Section 5 deals with the solution of the Dengue Fever SIR model by the proposed method. In Section 6, we state some theorems on Error Analysis. A comparison of solutions and Error Analysis through graphs and tables is shown in Section 7. This paper ends with a conclusion in Section 8.

2. Fundamental Definitions

In 1982, a French geophysical engineer, Jean Morlet, introduced the idea of a wavelet transform that provides a new mathematical procedure for seismic wave investigation [20]. A family of functions ψ a , b constructed by the translation and dilation of a single function ψ (the mother wavelet) is called wavelet, as follows:
ψ a , b ( x ) = 1 | a | ψ x b a , a , b R , a 0
where a and b are the dilation and translation parameters, respectively [21].

2.1. Fibonacci Polynomials

The following general formula is used to define Fibonacci polynomials [22], as follows:
F n ( x ) = 1 , n = 0 , x , n = 1 , x F n 1 ( x ) + F n 2 ( x ) , n > 1 .
Furthermore, the power form of these polynomials is as follows [22]:
F n ( x ) = i = 0 n 2 n i i x n 2 i , n 0 .

2.2. Fibonacci Wavelets

Fibonacci wavelets are defined on [0, 1] as follows [23]:
ψ n , m ( x ) = 1 w m 2 ( k 1 ) / 2 F m ( 2 k 1 x n + 1 ) , n 1 2 k 1 x < n 2 k 1 0 , otherwise ,
with
w m = 0 1 F m 2 ( x ) d x ,
where m = 0 , 1 , , M 1 is the degree of the Fibonacci polynomial F m ( x ) , and k = 1 , 2 , and n = 1 , 2 , , 2 k 1 denote the levels of the resolution and translation parameters, respectively. The coefficient 1 w m is referred to as normalization factor.
For example, for k = 1 and M = 10 , we have Fibonacci wavelets as
ψ 1 , 0 ( x ) = 1 , 0 x < 1 0 , otherwise ψ 1 , 1 ( x ) = 3 x , 0 x < 1 0 , otherwise ψ 1 , 2 ( x ) = 1 2 15 7 1 + x 2 , 0 x < 1 0 , otherwise ψ 1 , 3 ( x ) = 105 239 2 x + x 3 , 0 x < 1 0 , otherwise ψ 1 , 4 ( x ) = 3 35 1943 1 + 3 x 2 + x 4 , 0 x < 1 0 , otherwise ψ 1 , 5 ( x ) = 3 4 385 2582 3 x + 4 x 3 + x 5 , 0 x < 1 0 , otherwise ψ 1 , 6 ( x ) = 3 5005 1,268,209 1 + 6 x 2 + 5 x 4 + x 6 , 0 x < 1 0 , otherwise ψ 1 , 7 ( x ) = 3 5005 2,827,883 4 x + 10 x 3 + 6 x 5 + x 7 , 0 x < 1 0 , otherwise ψ 1 , 8 ( x ) = 3 2 85,085 28,195,421 1 + 10 x 2 + 15 x 4 + 7 x 6 + x 8 , 0 x < 1 0 , otherwise ψ 1 , 9 ( x ) = 3 1,616,615 5,016,284,989 5 x + 20 x 3 + 21 x 5 + 8 x 7 + x 9 , 0 x < 1 0 , otherwise ψ 1 , 10 ( x ) = 3 1,616,615 11,941,544,471 1 + 15 x 2 + 35 x 4 + 28 x 6 + 9 x 8 + x 10 , 0 x < 1 0 , otherwise .

3. Function Approximation

Every function f L 2 [ 0 , 1 ) approximately expanded via Fibonacci wavelets defined by (7) can be expressed as
f ( x ) n = 1 m = 0 c n , m ψ n , m ,
where c n , m = f ( x ) , ψ n , m ( x ) denotes Fibonacci wavelet coefficients. Considering the truncated series in (10), we obtain
f ( x ) = n = 1 2 k 1 m = 0 M 1 c n , m ψ n , m = C T Ψ ( x ) ,
where C is 2 k 1 M × 1 matrix given as
C = c 1 , 0 , c 1 , 1 , , c 1 , M 1 , c 2 , 0 , c 2 , 1 , , c 2 , M 1 , , c 2 k 1 , 0 , , c 2 k 1 , M 1 T = c 1 , c 2 , , c 2 k 1 T ,
and Ψ ( x ) appearing in Equation (11) is the Fibonacci wavelet matrix of the order 2 k 1 M × 1 given by
Ψ ( x ) = ψ 1 , 0 , ψ 1 , 1 , , ψ 1 , M 1 , ψ 2 , 0 , ψ 2 , 1 , , ψ 2 , M 1 , , ψ 2 k 1 , 0 , , ψ 2 k 1 , M 1 T = ψ 1 ( x ) , ψ 2 ( x ) , , ψ 2 k 1 ( x ) T .
We consider the collocation points as follows:
x l = l 0.5 2 k 1 M , l = 1 , 2 , , 2 k 1 M .

4. Operational Matrix of Integration (OMI)

To obtain the operational matrix of integration associated with Fibonacci wavelets given by Equation (7), we proceed with
0 x Ψ ( s ) d s = P Ψ ( x ) + P ¯ .
where P denotes the Fibonacci wavelet operational matrix of order 2 k 1 M × 2 k 1 M and P ¯ denotes the 2 k 1 M × 1 column matrix. we construct the operational matrix P for different values of k and M.
For k = 1 and M = 6, we have
0 x ψ 1 , 0 ( s ) d s = 0 1 3 0 0 0 0 Ψ ( x ) ,
0 x ψ 1 , 1 ( s ) d s = 3 2 0 7 5 0 0 0 Ψ ( x ) ,
0 x ψ 1 , 2 ( s ) d s = 0 5 6 7 0 239 42 0 0 Ψ ( x ) ,
0 x ψ 1 , 3 ( s ) d s = 105 2 239 0 7 2 239 0 1943 4 717 0 Ψ ( x ) ,
0 x ψ 1 , 4 ( s ) d s = 0 0 0 717 5 1943 0 4 2582 5 21,373 Ψ ( x ) ,
0 x ψ 1 , 5 ( s ) d s = 385 4 2582 0 0 0 21,373 24 2582 0 Ψ ( x ) + 1,268,209 24 33,566 ψ 1 , 6 ( x ) .
Here,
P = 0 1 3 0 0 0 0 3 2 0 7 5 0 0 0 0 5 6 7 0 239 42 0 0 105 2 239 0 7 2 239 0 1943 4 717 0 0 0 0 717 5 1943 0 4 2582 5 21,373 385 4 2582 0 0 0 21,373 24 2582 0 a n d P ¯ = 0 0 0 0 0 1,268,209 24 33,566 ψ 1 , 6 ( x ) .
Similarly, for k = 1 and M = 10 , we obtain
0 x ψ 1 , 0 ( s ) d s = 0 1 3 0 0 0 0 0 0 0 0 Ψ ( x ) ,
0 x ψ 1 , 1 ( s ) d s = 3 2 0 7 5 0 0 0 0 0 0 0 Ψ ( x ) ,
0 x ψ 1 , 2 ( s ) d s = 0 5 6 7 0 239 42 0 0 0 0 0 0 Ψ ( x ) ,
0 x ψ 1 , 3 ( s ) d s = 105 2 239 0 7 2 239 0 1943 4 717 0 0 0 0 0 Ψ ( x ) ,
0 x ψ 1 , 4 ( s ) d s = 0 0 0 717 5 1943 0 4 2582 5 21,373 0 0 0 0 Ψ ( x ) ,
0 x ψ 1 , 5 ( s ) d s = 385 4 2582 0 0 0 21,373 24 2582 0 1,268,209 24 33,566 0 0 0 Ψ ( x ) ,
0 x ψ 1 , 6 ( s ) d s = 0 0 0 0 0 4 33,566 7 1,268,209 0 2,827,883 7 1,268,209 0 0 Ψ ( x ) ,
0 x ψ 1 , 7 ( s ) d s = 3 5005 4 2,827,883 0 0 0 0 0 1,268,209 8 2,827,883 0 28,195,421 4 48,074,011 0 Ψ ( x ) ,
0 x ψ 1 , 8 ( s ) d s = 0 0 0 0 0 0 0 48,074,011 18 28,195,421 0 5,016,284,989 18 535,712,999 Ψ ( x ) ,
0 x ψ 1 , 9 ( s ) d s = 3 323,323 25,081,424,945 0 0 0 0 0 0 0 535,712,999 5 5,016,284,989 0 Ψ ( x ) + 11,941,544,471 10 5,016,284,989 ψ 1 , 10 ( x ) .
In this case,
P = 0 1 3 0 0 0 0 0 0 0 0 3 2 0 7 5 0 0 0 0 0 0 0 0 5 6 7 0 239 42 0 0 0 0 0 0 105 2 239 0 7 2 239 0 1943 4 717 0 0 0 0 0 0 0 0 717 5 1943 0 4 2582 5 21,373 0 0 0 0 385 4 2582 0 0 0 21,373 24 2582 0 1,268,209 24 33,566 0 0 0 0 0 0 0 0 4 33,566 7 1,268,209 0 2,827,883 7 1,268,209 0 0 3 5005 4 2,827,883 0 0 0 0 0 1,268,209 8 2,827,883 0 28,195,421 4 48,074,011 0 0 0 0 0 0 0 0 48,074,011 18 28,195,421 0 5,016,284,989 18 535,712,999 3 323,323 25,081,424,945 0 0 0 0 0 0 0 535,712,999 5 5,016,284,989 0
and P ¯ = 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 11,941,544,471 10 5,016,284,989 ψ 1 , 10 ( x ) T .

5. Stability Analysis and Solution of Dengue Fever SIR Model by Fibonacci Wavelet Collocation Method

Consider the Dengue Fever SIR model given by [6,7]
d S d t = μ h ( 1 S ( t ) ) α S ( t ) R ( t )
d I d t = α S ( t ) R ( t ) β I ( t )
d R d t = γ ( 1 R ( t ) ) I ( t ) δ R ( t )
with the initial conditions S ( 0 ) = 5,070,822 5,071,126 , I ( 0 ) = 304 5,071,126 , R ( 0 ) = 0.01 , α = 0.006, β = 0.333333, γ = 0.375 , δ = 0.02941 , and μ h = 0.0045 .
Here, R ˜ = α · γ β · δ is called the threshold parameter and R 0 = R ˜ is the basic reproductive number. If R ˜ < 1 , then the SIR model is stable [7]. In our case, the value of R ˜ = 0.2295 < 1 . Hence, our model is stable, and a solution is bounded for the parameter values on which we have obtained the numerical solution.
For k = 1 and M = 6 , we have n = 1 and m = 0 , 1 , , 5 .
Let us express
d S d t = C 1 T Ψ
d I d t = C 2 T Ψ
d R d t = C 3 T Ψ
where
C i = c i , 0 , c i , 1 , c i , 2 , c i , 3 , c i , 4 , c i , 5 T , i = 1 , 2 , 3
and
Ψ = ψ 1 , 0 , ψ 1 , 1 , ψ 1 , 2 , ψ 1 , 3 , ψ 1 , 4 , ψ 1 , 5 T .
Now, integrating from Equations (36)–(38) within 0 to t, and using initials conditions, we obtain
S ( t ) = C 1 T ( P Ψ + P ¯ ) + S ( 0 ) d T Ψ
I ( t ) = C 2 T ( P Ψ + P ¯ ) + I ( 0 ) d T Ψ
R ( t ) = C 3 T ( P Ψ + P ¯ ) + R ( 0 ) d T Ψ
where P and P ¯ are given by (22) and d = 1 2 1 , 0 , 0 , 0 , 0 , 0 T such that 1 = d T Ψ .
Now let us put the values of Equations (36)–(41) in Equations (33)–(35), and we obtain
C 1 T Ψ ( t ) = μ h d T Ψ ( t ) C 1 T P Ψ ( t ) C 1 T P ¯ S ( 0 ) d T Ψ ( t ) α C 1 T P Ψ ( t ) + C 1 T P ¯ + S ( 0 ) d T Ψ ( t ) C 3 T P Ψ ( t ) + C 3 T P ¯ + R ( 0 ) d T Ψ ( t )
C 2 T Ψ ( t ) = α C 1 T P Ψ ( t ) + C 2 T P ¯ + S ( 0 ) d T Ψ ( t ) C 3 T P Ψ ( t ) + C 3 T P ¯ + R ( 0 ) d T Ψ ( t ) β C 2 T P Ψ ( t ) + C 2 T P ¯ + I ( 0 ) d T Ψ ( t )
C 3 T Ψ ( t ) = γ d T Ψ ( t ) C 3 T P Ψ ( t ) C 3 T P ¯ R ( 0 ) d T Ψ ( t ) C 2 T P Ψ ( t ) + C 2 T P ¯ + I ( 0 ) d T Ψ ( t ) δ C 3 T P Ψ ( t ) + C 3 T P ¯ + R ( 0 ) d T Ψ ( t )
Let us consider the collocation points defined in Equation (14); we have τ 1 = 1 12 , τ 2 = 3 12 , τ 3 = 5 12 , τ 4 = 7 12 , τ 5 = 9 12 , τ 6 = 11 12 .
Now let us put these collocation points τ i , i = 1 , 2 , , 6 in Equations (42)–(44), and we obtain
C 1 T Ψ ( τ i ) = μ h d T Ψ ( τ i ) C 1 T P Ψ ( τ i ) C 1 T P ¯ S ( 0 ) d T Ψ ( τ i ) α C 1 T P Ψ ( τ i ) + C 1 T P ¯ + S ( 0 ) d T Ψ ( τ i ) C 3 T P Ψ ( τ i ) + C 3 T P ¯ + R ( 0 ) d T Ψ ( τ i )
C 2 T Ψ ( τ i ) = α C 1 T P Ψ ( τ i ) + C 1 T P ¯ + S ( 0 ) d T Ψ ( τ i ) C 3 T P Ψ ( τ i ) + C 3 T P ¯ + R ( 0 ) d T Ψ ( τ i ) β C 2 T P Ψ ( τ i ) + C 2 T P ¯ + I ( 0 ) d T Ψ ( τ i )
C 3 T Ψ ( τ i ) = γ d T Ψ ( τ i ) C 3 T P Ψ ( τ i ) C 3 T P ¯ R ( 0 ) d T Ψ ( τ i ) C 2 T P Ψ ( τ i ) + C 2 T P ¯ + I ( 0 ) d T Ψ ( τ i ) δ C 3 T P Ψ ( τ i ) + C 3 T P ¯ + R ( 0 ) d T Ψ ( τ i ) ,
for i = 1 , 2 , , 6 . Here, we have 18 equations in 18 unknowns, and these unknowns can be found by using the Newton–Raphson method, which are as follows:
C 1 = 1.0 × 10 4 × ( 0.590091969653335 , 0.003694941812179 , 0.000994140897631 ,
0.015915277834494 , 0.004831804737754 , 0.010287372858342 ) ,
C 2 = 1.0 × 10 4 × ( 0.114626280088155 , 0.281209022967569 , 0.923381047043343 ,
0.867607040144165 , 0.517344226006418 , 0.614352083709031 ) ,
C 3 = 1.0 × 10 3 × ( 0.163726424204976 , 0.043819756365797 , 0.924298713138464 ,
0.299576632108260 , 0.617152147379955 , 0.360188732065404 ) .
By substituting the values of unknowns in Equations (39)–(41), we obtain the approximated solution of the Dengue Fever SIR model.

6. Error Analysis

The convergence of the proposed method is deduced from the following results [24,25]:
Theorem 1.
“Let f C M [ 0 , 1 ) and f * = C T Ψ ( x ) be a Fibonacci wavelet expansion of the real sufficiently smooth function f ( x ) [ 0 , 1 ) ; then the norm of the truncated error e ( x ) can be bounded as
e ( x ) = f f * M M ! 2 M + 1 ,
where M is the maximal level of resolution”.
Theorem 2.
“If f ( x ) L 2 ( R ) is a continuous function defined on [ 0 , 1 ) and | f ( x ) | M , then the Fibonacci wavelet expansion of f ( x ) defined in Equation (10) converges uniformly to the function f ( x ) , and also,
| c | M δ 2 k 1 2 r 2 m ,
where δ = 1 w m ”.

7. Comparison of Solutions Obtained from Different Numerical Methods (FWCM, BWCM, and RK4) and Their Error Analysis

On an 11th Gen Intel(R) Core(TM) i5-1135G7 @ 2.40GHz 1.38 GHz processor with 8 GB RAM, we used MATLAB R2023b for graphical representation and for numerical simulation.
We give our computational results here along with a discussion of them, taking into account the mathematical model given by Equations (33)–(35) and its initial conditions with the following specified parameters: S ( 0 ) = 5,070,822 5,071,126 , I ( 0 ) = 304 5,071,126 , R ( 0 ) = 0.01 , α = 0.006 , β = 0.333333 , γ = 0.375 , δ = 0.02941 , and μ h = 0.0045 . After obtaining the solution from the FWCM, we examined it with the Runge–Kutta(RK4) and BWCM [13] solutions.
Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6 illustrate that the FWCM’s solutions are accurate and produce better approximations than the BWCM. In these tables, we have also compared Absolute Error (AE) of the proposed method and the BWCM with the RK4 method. As shown in Table 4, Table 5 and Table 6, Absolute Error (AE) decreases rapidly as we increase the value of M (maximal degree of the polynomial). We have also compared graphical solutions and Absolute Error (AE) obtained from the FWCM, BWCM, and RK4. Graphical representations of solutions obtained from the FWCM, BWCM, and RK4 are drawn in Figure 1, Figure 2 and Figure 3. From these figures, it is evident that the proposed method yields solutions that are significantly closer to the RK4 solution. A graphical comparison of Absolute Error obtained from the FWCM with RK4 and the BWCM with RK4 is shown in Figure 4, Figure 5 and Figure 6. The graph’s straight line indicates that the developed technique (FWCM) provides a better approximation than the BWCM method and produces a consistent absolute error at all places within the defined interval.

8. Conclusions

In this article, we have discussed a novel numerical approach called Fibonacci wavelet collocation method (FWCM) for the Dengue Fever model. Using Fibonacci wavelets, we constructed an operational matrix of integration that is used to estimate a dengue model. The efficiency and accuracy of the proposed method are demonstrated by comparing the obtained results with other existing methods, like Runge–Kutta (RK4) and BWCM [13]. The tables and figures clearly show that the FWCM outperforms the BWCM (when compared with RK4) and produces more consistent and accurate approximation results as M (maximal degree of the polynomial) increases. The error decreases more rapidly when we increase the degree (m) of the Fibonacci polynomials.

Author Contributions

All authors contributed to the study conception and design. Material preparation, data collection, and analysis were collaboratively performed by each author. The first draft of the manuscript was written by A.K. (Amit Kumar) under the supervision of A.K. (Ayub Khan) and A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Inquiries about data availability should go directly to the authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Guzman, M.G.; Harris, E. Dengue. Lancet 2015, 385, 453–465. [Google Scholar] [CrossRef] [PubMed]
  2. Gould, E.; Pettersson, J.; Higgs, S.; Charrel, R.; De Lamballerie, X. Emerging arboviruses: Why today? One Health 2017, 4, 1–13. [Google Scholar] [CrossRef] [PubMed]
  3. Bhatt, S.; Gething, P.W.; Brady, O.J.; Messina, J.P.; Farlow, A.W.; Moyes, C.L.; Drake, J.M.; Brownstein, J.S.; Hoen, A.G.; Sankoh, O.; et al. The global distribution and burden of dengue. Nature 2013, 496, 504–507. [Google Scholar] [CrossRef] [PubMed]
  4. Buhler, C.; Winkler, V.; Runge-Ranzinger, S.; Boyce, R.; Horstick, O. Environmental methods for dengue vector control—A systematic review and meta-analysis. PLoS Negl. Trop. Dis. 2019, 13, e0007420. [Google Scholar] [CrossRef] [PubMed]
  5. Simmons, C.P.; Farrar, J.J.; van Vinh Chau, N.; Wills, B. Dengue. N. Engl. J. Med. 2012, 366, 1423–1432. [Google Scholar] [CrossRef] [PubMed]
  6. Nuraini, N.; Soewono, E.; Sidarto, K. Mathematical model of dengue disease transmission with severe dhf compartment. Bull. Malays. Math. Sci. Soc. 2007, 30, 143–157. [Google Scholar]
  7. Yaacob, Y. Analysis of a dengue disease transmission model without immunity. Mat. Univ. Teknol. Malays. 2007, 23, 75–81. [Google Scholar]
  8. Khalid, M.; Sultana, M.; Khan, F.S. Numerical Solution of SIR Model of Dengue Fever. Int. J. Comput. Appl. 2015, 118, 1–4. [Google Scholar] [CrossRef]
  9. Rangkuti, Y.M.; Side, S.; Noorani, M.S.M. Numerical Analytic Solution of SIR Model of Dengue Fever Disease in South Sulawesi using Homotopy Perturbation Method and Variational Iteration Method. J. Math. Fundam. Sci. 2014, 46, 91–105. [Google Scholar] [CrossRef]
  10. Mungkasi, S. Improved variational iteration solutions to the SIR model of dengue fever disease for the case of South Sulawesi. J. Math. Fundam. Sci. 2020, 52, 297–311. [Google Scholar] [CrossRef]
  11. Umar, M.; Sabir, Z.; Raja, M.A.Z.; Sánchez, Y.G. A stochastic numerical computing heuristic of SIR nonlinear model based on dengue fever. Results Phys. 2020, 19, 103585. [Google Scholar] [CrossRef]
  12. Lede, Y.K.; Mungkasi, S. Performance of the Runge-Kutta methods in solving a mathematical model for the spread of dengue fever disease. AIP Conf. Proc. 2019, 2202, 020044. [Google Scholar]
  13. Li, F.; Baskonus, H.M.; Kumbinarasaiah, S.; Manohara, G.; Gao, W.; Ilhan, E. An Efficient Numerical Scheme for Biological Models in the Frame of Bernoulli Wavelets. Comput. Model. Eng. Sci. 2023, 137, 2381–2408. [Google Scholar] [CrossRef]
  14. Bulut, F.; Oruç, Ö.; Esen, A. Higher order Haar wavelet method integrated with strang splitting for solving regularized long wave equation. Math. Comput. Simul. 2022, 197, 277–290. [Google Scholar] [CrossRef]
  15. Keshavarz, E.; Ordokhani, Y.; Razzaghi, M. The Taylor wavelets method for solving the initial and boundary value problems of Bratu-type equations. Appl. Numer. Math. 2018, 128, 205–216. [Google Scholar] [CrossRef]
  16. Rahimkhani, P.; Ordokhani, Y.; Babolian, E. Müntz-Legendre wavelet operational matrix of fractional-order integration and its applications for solving the fractional pantograph differential equations. Numer. Algorithms 2017, 77, 1283–1305. [Google Scholar] [CrossRef]
  17. Sahu, P.K.; Ray, S.S. A new Bernoulli wavelet method for accurate solutions of nonlinear fuzzy Hammerstein–Volterra delay integral equations. Fuzzy Sets Syst. 2017, 309, 131–144. [Google Scholar] [CrossRef]
  18. Yousefi, S.; Banifatemi, A. Numerical solution of Fredholm integral equations by using CAS wavelets. Appl. Math. Comput. 2006, 183, 458–463. [Google Scholar] [CrossRef]
  19. Nemati, S.; Lima, P.M.; Torres, D.F.M. Numerical solution of a class of third-kind Volterra integral equations using Jacobi wavelets. Numer. Algorithms. 2021, 86, 675–691. [Google Scholar] [CrossRef]
  20. Debnath, L.; Shah, F.A. Lecture Notes on Wavelet Transforms; Springer: Cham, Switzerland, 2017. [Google Scholar] [CrossRef]
  21. Ahmad, K. Wavelet Packets and Their Statistical Applications; Springer: Singapore, 2018. [Google Scholar] [CrossRef]
  22. Falcón, S.; Plaza, Á. On k–Fibonacci sequences and polynomials and their derivatives. Chaos Solitons Fractals 2009, 39, 1005–1019. [Google Scholar] [CrossRef]
  23. Sabermahani, S.; Ordokhani, Y.; Yousefi, S. Fibonacci wavelets and their applications for solving two classes of time-varying delay problems. Optim. Control Appl. Methods. 2020, 41, 395–416. [Google Scholar] [CrossRef]
  24. Irfan, M.; Shah, F.A.; Nisar, K.S. Fibonacci wavelet method for solving Pennes bioheat transfer equation. Int. J. Wavelets Multiresolut. Inf. Process. 2021, 19, 2150023. [Google Scholar] [CrossRef]
  25. Rafiq, M.; Abdullah, A. Numerical Investigation using Fibonacci Wavelet Collocation Method for Solving Modified Unstable Nonlinear Schrödinger Equation. Int. J. Appl. Comput. Math 2023, 9, 118. [Google Scholar] [CrossRef]
Figure 1. Comparison of solution S ( t ) vs. t.
Figure 1. Comparison of solution S ( t ) vs. t.
Mathematics 12 02565 g001
Figure 2. Comparison of solution I ( t ) vs. t.
Figure 2. Comparison of solution I ( t ) vs. t.
Mathematics 12 02565 g002
Figure 3. Comparison of solution R ( t ) vs. t.
Figure 3. Comparison of solution R ( t ) vs. t.
Mathematics 12 02565 g003
Figure 4. Comparison and error plots for S ( t ) vs. t.
Figure 4. Comparison and error plots for S ( t ) vs. t.
Mathematics 12 02565 g004
Figure 5. Comparison and error plots for I ( t ) vs. t.
Figure 5. Comparison and error plots for I ( t ) vs. t.
Mathematics 12 02565 g005
Figure 6. Comparison and error plots for R ( t ) vs. t.
Figure 6. Comparison and error plots for R ( t ) vs. t.
Mathematics 12 02565 g006
Table 1. Solution S ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 6 , and their error analysis.
Table 1. Solution S ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 6 , and their error analysis.
tFWCM
For M = 6
RK4BWCM [13]
for M = 6
AE of FWCM with RK4
1.0 × 10 5 ×
AE of BWCM with RK4
1.0 × 10 3 ×
00.9999385215689020.9999400527614580.99994005270.1531192550.000000061
0.10.9999326059794050.9999340895898390.99988016320.1483610430.053926389
0.20.9999267138835460.9999281453076370.99982055240.1431424090.107592907
0.30.9999208446807900.9999222197743980.99976104120.1375093600.161178574
0.40.9999149981909810.9999163128532280.99970181240.1314662240.214500453
0.50.9999091746076550.9999104244106770.99964276320.1249803020.267661210
0.60.9999033744513450.9999045543166230.99958384570.1179865270.320708616
0.70.9998975985228900.9998987024441660.99952524210.1103921270.373460344
0.80.9998918478567460.9998928686695240.99946672140.1020812770.426147269
0.90.9998861236742910.9998870528719290.99940847410.0929197630.478578771
1.00.9998804273371350.9998812549335270.99935034120.0827596390.530913733
Table 2. Solution I ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 6 , and their error analysis.
Table 2. Solution I ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 6 , and their error analysis.
tFWCM
for M = 6
1.0 × 10 3 ×
RK4BWCM [13]
for M = 6
AE of FWCM with RK4
1.0 × 10 4 ×
AE of BWCM with RK4
00.1513885958022870.5994723854228820.10000000000.9144135720.099940052
0.10.1546756889817610.6387462429032760.09970415740.9080106460.099640282
0.20.1574705989496610.6765731246118320.09942163250.8981328640.099353975
0.30.1597361823479670.7130017660040470.09913523650.8843600570.099063936
0.40.1613979100606830.7480792704099030.09885632510.8658998300.098781517
0.50.1623388679077790.7818511636709110.09857236520.8415375150.098494180
0.60.1623947573391470.8143614469482170.09829425730.8095861260.098212821
0.70.1613488961285390.8456526477630170.09802125470.7678363130.097936689
0.80.1589272190675190.8757658724313730.09774763250.7135063210.097660055
0.90.1547932786594090.9047408382299790.09747642510.6431919480.097385951
1.00.1485432458132390.9326159505101280.09720805240.5528165070.097114790
Table 3. Solution R ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 6 , and their error analysis.
Table 3. Solution R ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 6 , and their error analysis.
tFWCM
for M = 6
RK4BWCM [13]
for M = 6
AE of FWCM with RK4
1.0 × 10 3 ×
AE of BWCM with RK4
00.0094638880960410.0100000000000000.00005994720.5361119030.009940052
0.10.0094385694268230.0099729287940270.00011690120.5343593670.009856027
0.20.0094156179233930.0099460800472230.00017181390.5304621230.009774266
0.30.0093953621662630.0099194478400340.00022475370.5240856730.009694694
0.40.0093783274253050.0098930264474220.02019578690.5146990220.010302760
0.50.0093652952976590.0098668103323780.00032497760.5015150340.009541832
0.60.0093573633456300.0098407941396410.00037238770.4834307940.009468406
0.70.0093560047345950.0098149726896380.00041807700.4589679550.009396895
0.80.0093631278709050.0097893409726150.00046210330.4262131010.009327237
0.90.0093811360397870.0097638941429740.00050452260.3827581030.009259371
1.00.0094129870432460.0097386275137920.00054538890.3256404700.009193238
Table 4. Solution S ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 10 , and their error analysis.
Table 4. Solution S ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 10 , and their error analysis.
tFWCM
for M = 10
RK4BWCM [13]
for M = 10
AE of FWCM with RK4
1.0 × 10 13 ×
AE of BWCM with RK4
1.0 × 10 3 ×
00.9999400527613810.9999400527614580.99994005300.7682743330.000000238
0.10.9999340895897610.9999340895898390.99988018610.7782663400.053903489
0.20.9999281453075590.9999281453076370.99982052330.7849276780.107622007
0.30.9999222197743200.9999222197743980.99976106280.7827072320.161156974
0.40.9999163128531510.9999163128532280.99970180290.7771561170.214509953
0.50.9999104244106010.9999104244106770.99964274190.7571721020.267682510
0.60.9999045543165500.9999045543166230.99958387830.7227551890.320676016
0.70.9998987024441000.9998987024441660.99952521030.6572520300.373492144
0.80.9998928686694690.9998928686695240.99946673660.5495603970.426132069
0.90.9998870528718910.9998870528719290.99940845560.3730349360.478597271
1.00.9998812549335180.9998812549335270.99935036590.0877076180.530889033
Table 5. Solution I ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 10 , and their error analysis.
Table 5. Solution I ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 10 , and their error analysis.
tFWCM
for M = 10
1.0 × 10 4 ×
RK4
1.0 × 10 4 ×
BWCM [13]
for M = 10
AE of FWCM with RK4
1.0 × 10 12 ×
AE of BWCM with RK4
00.5994723844282160.5994723854228820.00005994720.0994666520.000000000
0.10.6387462424059690.6387462429032760.00011690120.0497306320.000053026
0.20.6765731245825760.6765731246118320.00017181390.0029256110.000104156
0.30.7130017664166210.7130017660040470.00022475370.0412574430.000153453
0.40.7480792712424990.7480792704099030.02019578690.0832595440.020120978
0.50.7818511649080850.7818511636709110.00032497760.1237174420.000246792
0.60.8143614485839020.8143614469482170.00037238770.1635684310.000290951
0.70.8456526498051320.8456526477630170.00041807700.2042115390.000333511
0.80.8757658718059000.8757658693284210.00046210330.2477479070.000374526
0.90.9047408412032900.9047408382299790.00050452260.2973311290.000414048
1.00.9326159540868190.9326159505101280.00054538890.3576691080.000452127
Table 6. Solution R ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 10 , and their error analysis.
Table 6. Solution R ( t ) obtained by RK, FWCM, and BWCM for k = 1 , M = 10 , and their error analysis.
tFWCM
for M = 10
RK4BWCM [13]
for M = 10
AE of FWCM with RK4
1.0 × 10 12 ×
AE of BWCM with RK4
00.0099999999998130.0100000000000000.10000000000.1871159470.090000000
0.10.0099729287937810.0099729287940270.09970931880.2459647060.089736390
0.20.0099460800469230.0099460800472230.09942137790.3001141000.089475297
0.30.0099194478396840.0099194478400340.09913610270.3494600440.089216654
0.40.0098930264470290.0098930264474220.09885342110.3936000830.088960394
0.50.0098668103319460.0098668103323780.09857326330.4317050180.088706452
0.60.0098407941391790.0098407941396410.09829556180.4623038060.088454767
0.70.0098149726891550.0098149726896380.09802025120.4829522190.088205278
0.80.0097893409721250.0097893409726150.09774726820.4897089670.087957927
0.90.0097638941424980.0097638941429740.09747655170.4764157810.087712657
1.00.0097386275133580.0097386275137920.09720804250.4336322960.087469414
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

Kumar, A.; Khan, A.; Abdullah, A. Fibonacci Wavelet Collocation Method for Solving Dengue Fever SIR Model. Mathematics 2024, 12, 2565. https://doi.org/10.3390/math12162565

AMA Style

Kumar A, Khan A, Abdullah A. Fibonacci Wavelet Collocation Method for Solving Dengue Fever SIR Model. Mathematics. 2024; 12(16):2565. https://doi.org/10.3390/math12162565

Chicago/Turabian Style

Kumar, Amit, Ayub Khan, and Abdullah Abdullah. 2024. "Fibonacci Wavelet Collocation Method for Solving Dengue Fever SIR Model" Mathematics 12, no. 16: 2565. https://doi.org/10.3390/math12162565

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