Computer Modeling in Engineering & Sciences

DOI: 10.32604/cmes.2022.019154

ARTICLE

Efficient Numerical Scheme for the Solution of HIV Infection CD4+ T-Cells Using Haar Wavelet Technique

Rohul Amin1, Şuayip Yüzbaş2,* and Shah Nazir3

1Department of Mathematics, University of Peshawar, Peshawar, 25120, Pakistan
2Department of Mathematics, Faculty of Science, Akdeniz University, Antalya, Tr-07058, Turkey
3Department of Computer Science, University of Swabi, Swabi, 23430, Pakistan
*Corresponding Author: Şuayip Yüzbaş. Email: suayipyuzbasi@gmail.com
Received: 06 September 2021; Accepted: 05 November 2021

1  Introduction

Many models have been developed by mathematicians in the last decade to explain the immunological response to Human Immunodeficiency Virus (HIV) infection. Due to the scarcity of CD4+ T-cells, HIV disease is considered to result in concealment of the immune system (referred to generally as T4-cells or T-helper cells), cells which play a focal part in the human immune system. A class of white blood cells called CD4+ T-cells is essentially infected by HIV and this selective depletion of CD4+ T-cells means that a focal part in the resistant direction fills in as a clinical marker to estimate HIV disease movement. Mathematical models play an important role in the dynamics of these infectious diseases [1–6]. The level of CD4+ T-cells in the fringe blood is controlled at a level between 800 and 1200 mm−3 in a normal human body. These cells are the most inexhaustible white blood cell of the human safe framework, which battle against diseases [7]. Many authors introduced different numerical methods for the investigation of CD4+ T-cells [8 –14].

In recent years, many researchers have studied on new analytical and numerical methods for model problems characterized by differential equations. Yüzbaşi et al. [15] presented the Pelllucas collocation method for the solutions of two population models and residual correction. Mallawi et al. [16] utilized Legendre collocation method to the spacetime variable fractional-order advectiondispersion equation. Yüzbaşi [17] used an operational method for solutions of Riccati type differential equations. Chu et al. [18] studied the generalized (2+1) dimensional shallow water equation. The solution for the fractional system of HIV-1 infection of CD4+ T-cells was studied by Khater et al. [19] by using AtanganaBaleanu fractional derivative. Khater et al. [20] developed a semi analytical and numerical scheme for a biological model. Khater et al. [21] used the trigonometric quintic B-spline technique for the solutions of complex nonlinear FokasLenells equations. In reference [22] the nonlinear phi-four equation is solved through two analytical and semi-analytical techniques. Khater et al. [23] used trigonometric Quintic B-spline method for the solution of conformable fractional nonlinear time-space telegraph equation. Khater et al. [24] investigated the analytical and numerical solutions of the modified BenjaminBonaMahony equation via the modified Bspline collection method. Khater et al. [25] found the solution of nonlinear KleinFockGordon equation by using generalized exponential function and generalized Riccati expansion methods. Khater et al. [26] investigated the analytical and semi-analytical solutions of the time-fractional CahnAllen equation by using the Adomian decomposition method. Khater et al. [27] found the analytical solutions of the nonlinear Schrodinger equation with the higher-order through Kudryashov method. They also found the solutions of quadratic cubic fractional nonlinear Schrodinger equation by Adomian decomposition process [28]. Khater et al. [29] used the trigonometric quintic and exponential cubic B-spline schemes for the solutions of the nonlinear Klein-Gordon-Zakharov model. Yue et al. [30] found a solution of the fractional nonlinear HirotaSatsumaShallow water wave equation by using a modified Kudryashov method. Khater et al. [31] found solutions of the Fisher-Kolmogorov-Petrovskii-Piskunov model by employing the modified Kudryashov and trigonometric-quantic B-spline methods. Li et al. [32] found wave solutions of the (2+1) dimensional Kadomtsev Petviashvili Benjamin Bona Mahony model. Smadi et al. [33] study the accuracy of solution for fractional order an SEIR epidemic model by using the homotopy analysis method. Freihet et al. [34] found a solution of a fractional stiff system using residual functions algorithm. Smadi et al. [35] developed analytical technique for coupled system of fractional partial differential equations, for solutions of nonlinear fractional Kundu-Eckhaus equations [36] and for solution of coupled fractional resonant Schrodinger equations [37]. The same authors analyzed and investigated the analytical solution of the seventh-order fractional Sawada Kotera Ito, Lax, and Kaup Kupershmidt equations [38]. Moreover, the authors also analyzed and studied fuzzy fractional differential equations in terms of Atangana-Baleanu Caputo differential operators equipped with uncertain constraints coefficients and initial conditions [39].

Here we develop an accurate scheme by using HWC technique for the solution of the HIV infection of CD4+ T-cells. The main contributions of this work as:

•   To develop efficient numerical scheme by utilizing HWC technique for HIV infection CD4+ T-cells

•   To design numerical scheme using HWC technique

•   To evaluate the efficacy of the established technique in some instances and compare the results with other techniques Runge-Kutta, LADM-Pade [1], Bessel collocation technique [2], PIA(1,1) [3], MVIM [4], DTM [7] and Galerkin technique [5] available in the literature

In this article we find the numerical solution of first order model of the form [5]

{dTdt=λ−αT+rT(1−T+ITmax)−k∗VT,dIdt=k∗VT−βI,dVdt=N∗βI−γV,(1)

with initial conditions:

T(0)=T0,I(0)=I0andV(0)=V0.(2)

Here T(t), I(t), V(t) are used for concentration of uninfected cells, infected cells and free virus particles of CD4+ T-cells by HIV in the blood. rT(1−T+ITmax) is logistic growth of the healthy cells, Tmax is the most extreme level of cells in the human body, k is the steady rate which the body produces cells, k* V T is the frequency of HIV infection of healthy cells, k* > 0 is the rate of virus infection, N*β is the rate of virus production by contaminated cells, where N* is the average number of particles of infection produced by the infected T-cell, and γ is the rate of death of particles of the virus.

The paper is structured as: In Section 2, Haar functions are defined. HWC technique for solution of HIV infection is given in Section 3. In Section 4 error estimation about the model is given. In Section 5, one example is given. Conclusion is given in the last Section 6.

2  Haar Wavelet

Here we discuss Haar functions, integration of Haar functions and collocation points.

Definition 2.1. Scaling function on [α1, α2) is [40]

h1(t)={1fort∈[α1,α2),0elsewhere.(3)

Mother wavelet on [α1, α2) is

h2(t)={1fort∈[α1,α1+α22),−1fort∈[α1+α22,α2),0elsewhere.(4)

The other terms can be written as

hi(t)={1fort∈[η1,η2),−1fort∈[η2,η3),0elsewhere,(5)

where η1=α1+(α2−α1)ζd,η2=α1+(α2−α1)ζ+0.5d,η3=α1+(α2−α1)ζ+1d, where d = 2r, and r=0,1,…,d−1. If we take interval [0, 1], then values of η1, η2 and η3 are: η1=ζd, η2=1/2+ζd, η3=1+ζd. Any member u(t) in L2[0, 1), is written as u(t)=∑k=1∞λkhk(t), we truncate this series is u(t)≈∑k=1Nλkhk(t).

Using the notation

pi,1(t)=∫0thi(x)dx,(6)

and

pi,1(t)={t−ρ1fort∈[η1,η2),η3−tfort∈[η2,η3),0elsewhere.(7)

Generally,

pi,n(t)=∫0tpi,n−1(x)dx.(8)

Thus pi, n(t) is obtained as under [40]

pi,n(t)={0fort∈[0,ρ1),(t−ρ1)nn!fort∈[ρ1,ρ2),[(t−ρ1)n−2(ρ1−ρ2)n]n!fort∈[ρ2,ρ3),1n![(t−ρ1)n−2(ρ1−ρ2)n+(t−ρ3)n],fort∈[ρ3,1).(9)

Definition 2.2. The [β1, β2] interval for HWC scheme is discretized as [40]

ti=β1+(β2−β1)i−1/22Mi=1,2,3,4,…,2M=N.(10)

In the above Eq. (10), a collocation point (CP) are defined. Some of the recent work using HWC technique can be seen in [41–47].

3  Numerical Method

The implementation of the HWC method is discussed in this section in order to find the HIV model solution provided in Eq. (1). Using Haar functions, the derivative of the unknown function in the method is approximated and the expression for the unknown function is obtained by using initial condition and integration. By applying the Haar technique to Eq. (1) and putting the CPs, we get a system of algebraic equations. The Broyden technique is used to find solution of this system. At last, the approximate solution at CPs is obtained using these coefficients. We use the symbols Θ1=∑i=1N and Θ2=∑i=1M.

First, we consider that T′(t), I′(t) and V′(t) are square integrable function and hence can be written as

T′(t)=Θ1aihi(t),I′(t)=Θ1bihi(t)andV′(t)=Θ1cihi(t).(11)

Integrating with respect to t, we have

TN(t)=T0+Θ1aipi,1(t),IN(t)=I0+Θ1bipi,1(t),andVN(t)=V0+Θ1cipi,1(t).(12)

Putting Eq. (11) and Eq. (12) in Eq. (1), we have

Θ1aihi(t)=λ−α(T0+Θ1aipi,1(t))−k∗(v0+Θ1cipi,1(t))(T0+Θ1aipi,1(t))+r(T0+Θ1aipi,1(t))[1−(T0+Θ1aipi,1(t))(I0+Θ1bipi,1(t))Tmax],Θ1bihi(t)=k∗(v0+Θ1cipi,1(t))(T0+Θ1aipi,1(t))−β(I0+Θ1bipi,1(t)),Θ1cihi(t)=N∗β(I0+Θ1bipi,1(t))−γ(V0+Θ1cipi,1(t)),

by simplification we have

Θ1aihi(t)−λ+α(T0+Θ1aipi,1(t))+k∗(v0+Θ1cipi,1(t))(T0+Θ1aipi,1(t))−r(T0+Θ1aipi,1(t))[1−(T0+Θ1aipi,1(t))(I0+Θ1bipi,1(t))Tmax]=0,Θ1bihi(t)−k∗(v0+Θ1cipi,1(t))(T0+Θ1aipi,1(t))+β(I0+∑i=1Nbipi,1(t))=0,Θ1cihi(t)−N∗β(I0+Θ1bipi,1(t))+γ(V0+Θ1cipi,1(t))=0,

putting the CPs (10), we have

F1,j=Θ1aihi(t)−λ+α(T0+Θ1aipi,1(t))+k∗(v0+Θ1cipi,1(t))(T0+Θ1aipi,1(t))−r(T0+Θ1aipi,1(t))[1−(T0+Θ1aipi,1(t))(I0+Θ1bipi,1(t))Tmax],F2,j=Θ1bihi(t)−k∗(v0+Θ1cipi,1(t))(T0+Θ1aipi,1(t))+β(I0+Θ1bipi,1(t))=0,F3,j=Θ1cihi(t)−N∗β(I0+Θ1bipi,1(t))+γ(V0+Θ1cipi,1(t))=0.

This is solved by Broyden’s method. Jacobian is

J=[Jjk]3N×3N,(13)

where

Jjk={∂F1,j∂ak=hk(tj)+αpk,1(tj)−rpk,1(tj)−k∗[V0pk,1(tj)+Θ1cipi,1(tj)pk,1(tj)+rTmax[(T0+Θ1aipi,1(tj))2pk,1(tj)+I0pk,1(tj)+pk,1(tj)Θ1cipi,1(tj)],∂F1,j∂bk=T0pk,1(tj)+Θ1aipi,1(tj)pk,1(tj),∂F1,j∂ck=−k∗T0pk,1(tj)+Θ1aipi,1(tj)pk,1(tj),∂F2,j∂ak=−k∗V0pk,1(tj)(1+Θ1pk,1(tj)),∂F2,j∂bk=hk(tj)+βpk,1(tj),∂F2,j∂ck=−k∗I0pk,1(tj)(1+Θ1pk,1(tj)),∂F3,j∂ak=0,∂F3,j∂bk=−N∗βpk,1(tj),∂F3,j∂ck=hk(tj)+γpk,1(tj).

The solution of this gives the values of unknown coefficients ai’s, bi’s and ci’s. The required solution TN(t), IN(t) and VN(t) at CPs is calculated by putting ai, bi ci’s in Eq. (12).

4  Error Estimation

Here, we study the residual error estimation [5] for HIV model (1) utilizing HWC technique. The residual functions R1, N(t), R2, N(t) and R3, N(t) are

{R1,N(t)=TN′−λ+αTN−rTN(1−TN+INTmax)+k∗VNTN,R2,N(t)=IN′−k∗VNTN+βIN,R3,N(t)=VN′−N∗βIN+γVN,(14)

Define the error function as

{e1,N(t)=T(t)−TN(t),e2,N(t)=I(t)−IN(t),e2,N(t)=V(t)−VN(t),(15)

where T(t), I(t) and V(t) are exact solutions. So, we have

{T(t)=e1,N(t)+TN(t),I(t)=e2,N(t)+IN(t),V(t)=e3,N(t)+VN(t).(16)

Also, we can write

{T′(t)−TN′(t)=[T(t)−TN(t)]′=[e1,N(t)]′I′(t)−IN′(t)=[I(t)−IN(t)]′=[e2,N(t)]′V′(t)−VN′(t)=[V(t)−VN(t)]′=[e3,N(t)]′.(17)

By subtracting system (14) from system (1), we have

{T′(t)−TN′(t)=−α[T(t)−TN(t)]+r[T(t)−TN(t)]−rTmax[T2(t)−TN2(t)+T(t)I(t)−TN(t)IN(t)]−k∗[V(t)T(t)−VN(t)TN(t)]−R1,N(t),I′(t)−IN′(t)=k∗[V(t)T(t)−VN(t)TN(t)]−β[I(t)−IN(t)]−R2,N(t),V′(t)−VN′(t)=N∗β[I(t)−IN(t)]−γ[V(t)−VN(t)]−R3,N(t).(18)

By using systems (5)–(6) in Eq. (7) and simplifying, we obtain

{[e1,N(t)]′=−αe1,N(t)+re1,N(t)−rTmax[e1,N2(t)+e1,N(t)e2,N(t)+e1,N(t)IN+e2,N(t)TN]−k∗[e1,N(t)e3,N(t)+e3,N(t)TN+e1,N(t)VN]−R1,N(t),[e2,N(t)]′=k∗[e1,N(t)e3,N(t)+e3,N(t)TN+e1,N(t)VN]−βe2,N(t)−R2,N(t),[e3,N(t)]′=N∗βe2,N(t)−γe3,N(t)−R3,N(t).(19)

where e1, N(t), e2, N(t) and e3, N(t) are unknowns functions. The initial conditions for approximate solution TN(t), IN(t) and VN(t) are

TN(0)=T0,IN(0)=I0andVN(0)=V0,(20)

so initial conditions for system (19) are e1, N(0) = 0, e2, N(0) = 0 and e3, N(0) = 0, where e1, N(t), e2, N(t) and e3, N(t) are estimated by e1, N, M(t), e2, N, M(t) and e3, N, M(t) the Haar wavelet technique.

Let [e1,N,M(t)]′, [e2,N,M(t)]′ and [e3,N,M(t)]′ are in L2[0, 1), so

[e1,N,M(t)]′=Θ2ξihi(t),[e2,N,M(t)]′=Θ2μihi(t)and[e3,N,M(t)]′=Θ2νihi(t).(21)

Integrating the above system (21), with respect to t we obtain the following expression:

e1,N,M(t)=Θ2ξipi,1(t),e2,N,M(t)=Θ2μipi,1(t),ande3,N,M(t)=Θ2νipi,1(t).(22)

Applying Haar approximations, we have

Θ2ξihi(t)=−αΘ2ξipi,1(t)−k∗[Θ2ξipi,1(t)Θ2νipi,1(t)+TNΘ2νipi,1(t)+VNΘ2ξipi,1(t)]−rTmax[(Θ2ξipi,1(t))2+Θ2ξipi,1(t)Θ2μipi,1(t)+INΘ2ξipi,1(t)+TNΘ2μipi,1(t)]−R1,N(t),Θ2μihi(t)=k∗[Θ2ξipi,1(t)Θ2νipi,1(t)+TNΘ2νipi,1(t)+VNΘ2ξipi,1(t)]−βΘ2μipi,1(t)−R2,N(t),Θ2νihi(t)=N∗βΘ2μipi,1(t)−γΘ2νipi,1(t)−R3,N(t),

After simplification, we have

F1,j=Θ2ξihi(t)+αΘ2ξipi,1(t)+k∗[Θ2ξipi,1(t)Θ2νipi,1(t)+TNΘ2νipi,1(t)+VNΘ2ξipi,1(t)]+rTmax[(Θ2ξipi,1(t))2+Θ2ξipi,1(t)Θ2μipi,1(t)+INΘ2ξipi,1(t)+TNΘ2μipi,1(t)]+R1,N(t)=0,F2,j=Θ2μihi(t)−k∗[Θ2ξipi,1(t)Θ2νipi,1(t)+TNΘ2νipi,1(t)+VNΘ2ξipi,1(t)]+βΘ2μipi,1(t)+R2,N(t)=0,F3,j=Θ2νihi(t)−N∗βΘ2μipi,1(t)+γΘ2νipi,1(t)+R3,N(t),

putting the discrete CPs (10), we obtain

F1,j=Θ2ξihi(tj)+αΘ2ξipi,1(tj)+k∗[Θ2ξipi,1(tj)Θ2νipi,1(tj)+TNΘ2νipi,1(tj)+VNΘ2ξipi,1(tj)]+rTmax[(Θ2ξipi,1(tj))2+Θ2ξipi,1(tj)Θ2μipi,1(tj)+INΘ2ξipi,1(tj)+TNΘ2μipi,1(tj)]+R1,N(tj),F2,j=Θ2μihi(tj)−k∗[Θ2ξipi,1(tj)Θ2νipi,1(tj)+TNΘ2νipi,1(tj)+VNΘ2ξipi,1(tj)]+βΘ2μipi,1(tj)+R2,N(tj),F3,j=Θ2νihi(tj)−N∗βΘ2μipi,1(tj)+γΘ2νipi,1(tj)+R3,N(tj).

Broyden method is used for solution of above system. Jacobian is

J=[Jjk]3M×3M,(23)

where

Jjk={∂F1,j∂ξk=hk(tj)+αpk,1(tj)+k∗[pk,1(tj)Θ2νipi,1(tj)+VNpk,1(tj)]+rTmax[2Θ2ξipi,1(tj)pk,1(tj)+pk,1(tj)Θ2pi,1(tj)+INpk,1(tj)],∂F1,j∂μk=rTmax[Θ2ξipi,1(tj)pk,1(tj)+TNpk,1(tj)],∂F1,j∂νk=−k∗[Θ2xiipi,1(tj)pk,1(tj)+TNpk,1(tj)],∂F2,j∂ξk=−k∗[pk,1(tj)Θ2νpk,1(tj)+VNpk,1(tj)],∂F2,j∂μk=hk(tj)+βpk,1(tj),∂F2,j∂νk=−k∗[pk,1(tj)Θ2ξipi,1(tj)+TNpk,1(tj)],∂F3,j∂ξk=0,∂F3,j∂μk=−N∗βpk,1(tj),∂F3,j∂νk=hk(tj)+γpk,1(tj).

The unknown coefficients ξi, μi and νi are obtained from the solution of this system. The approximate solution at discrete CPs is obtained by plugging ξi, νi νii=1,2,…M in Eq. (22). Substituting the values of e1, N, M(t), e2, N, M(t) and e3, N, M(t) in system (16), we get the required solution.

5  Numerical Applications

The performance of the HWC technique is tested on example in this section. The numerical results are compared with Runge-Kutta technique, Bessel collocation technique, LADM-Pade and Galerkin technique available in the literature.

Problem 1. Consider the following system [5]

dTdt=0.1−0.02+T(1−T+I1500)−0.0027VTdIdt=−0.0027VT−0e.3I,dVdt=10(0.3)I−2.4V,(24)

with T(0) = 0.1, I(0) = 0 and V(0) = 0.1. The interval of study is 0≤ t ≤ 1.

The first order derivatives dTdt, dIdt and dVdt in above system (24) are approximated by Haar functions. Let dTdt, dIdt and dVdt are square integrable functions then

dT(t)dt=Θ1aihi(t),dI(t)dt=Θ1bihi(t)anddV(t)dt=Θ1cihi(t).(25)

By using initial conditions and integration, we obtain the approximate solution of this system in terms of Haar functions

T(t)=0.1+Θ1aipi,1(t),I(t)=Θ1bipi,1(t),andV(t)=0.1+Θ1cipi,1(t).(26)

Putting these approximations and CPs in above system (24), we obtain a system of nonlinear algebraic equations which is then solved by the method of Broyden’s. Error estimation is also calculated in a similar way, which aims to reduce the error of numerical solution. The residual functions R1, N(t), R2, N(t) and R3, N(t) are calculated as discussed in Eq. (14). The errors functions e1, N(t), e2, N(t), and e3, N(t) are obtained as discussed in Eq. (15).

The error functions for distinct CPs and distinct values of time are given in Tables 1–3. The error functions obtained for variable N = 3 at different values of t are compared with Laplace Adomian decomposition technique, Runge-Kutta technique, modified variational iteration technique, Pade approximation, the perturbation-iteration algorithm, Bessel collocation technique, differential transform technique and exponential Galerkin technique available in literature. From the tables we see that as N increases, the values of our results near to those of the other techniques. This show that our approximate solution become accurate as the number of discrete CPs N increases. Even a batter accuracy is obtained by taking more discrete CPs. The graph of numerical solution is also given in Fig. 1. The estimated results of error functions for distinct number of CPs are shown in Fig. 2. Due to the simplicity of the Haar wavelet it is effective for solution of the first-order of HIV infection CD4+ T-Cells model. However, HWC scheme has disadvantages too. This method use constant box functions and due to this we need a large number of collocation points in order to achieve better accuracy. This disadvantage can be overcome if Haar wavelet is replaced with some other wavelets having better approximating properties.

Figure 1: Numerical solution of problem 1

Figure 2: Estimated error functions for problem 1

6  Conclusion

Haar collocation scheme is developed for the solution of the HIV CD4+ T-cells model. Also, we discussed a procedure known as residual error estimation, whose aim is to get better arrangements utilizing the obtained solution. The technique is tested on one example, and the results are compared with other methods available in the literature. The comparison of the present HWC technique with Runge-Kutta technique, Bessel collocation technique, LADM-Pade and Galerkin technique is given. The results demonstrate that Haar technique is effective and precise for distinct numbers of CPs. The results show that taking the large value of variable N, the HWC technique give the accurate results of the numerical solution. The error functions using various numbers of CPs are also calculated and reported in the table. From results, it is seen that proposed technique gives good results for this model. The proposed technique is easily implemented in any software packages. MATLAB software is used to obtain the numerical results.

Acknowledgement: The authors are very grateful to the anonymous referees for their constructive comments which have improved the paper.

Funding Statement: No funding is available to support this study.

Conflicts of Interest: On behalf of all authors, the corresponding author states that there is no conflict of interest.

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.