Research article

Modeling the variable transmission rate and various discharges on the spread of Malaria


  • Received: 15 August 2022 Revised: 10 October 2022 Accepted: 17 October 2022 Published: 01 November 2022
  • Natural and household discharges are the natural breeding grounds of various mosquito species, including female Anopheles mosquitoes, which transmit the Plasmodium parasite, causing the spread of the life-threatening disease malaria. Apart from that, population migrations also have a substantial impact on malaria transmission, claiming about half a million lives every year around the world. To assess the effects of the cumulative density of households and other natural discharges, and emigration-dependent interaction rates on the dissemination of the vector-borne infectious disease malaria, we propose and analyze a non-linear mathematical model. The model comprises five dependent variables, namely, the density of the susceptible human population, the density of the infective human population, the density of the susceptible female Anopheles mosquito population, the density of the infective mosquito population and cumulative density of household and other natural discharges. In the model, the density of the mosquito population is supposed to follow logistic growth, whose intrinsic growth rate is a linear function of the cumulative density of household and other natural discharges. The nonlinear model is analyzed by using the stability theory of differential equations, numerical simulations and sensitivity analysis. The analysis shows that an increase in non-emigrating population causes increased incidence of malaria. It is also found that an increase in household and other natural discharges accelerates the occurrence of malaria. A basic differential sensitivity analysis is carried out to assess the sensitivity of model solutions with respect to key parameters. The model's numerical simulations demonstrate the analytical findings.

    Citation: Jitendra Singh, Maninder Singh Arora, Sunil Sharma, Jang B. Shukla. Modeling the variable transmission rate and various discharges on the spread of Malaria[J]. Electronic Research Archive, 2023, 31(1): 319-341. doi: 10.3934/era.2023016

    Related Papers:

    [1] Yi Liu, Lei Gao, Tianxu Zhao . Partially doubly strictly diagonally dominant matrices with applications. Electronic Research Archive, 2023, 31(5): 2994-3013. doi: 10.3934/era.2023151
    [2] Natália Bebiano, João da Providência, Wei-Ru Xu . Approximations for the von Neumann and Rényi entropies of graphs with circulant type Laplacians. Electronic Research Archive, 2022, 30(5): 1864-1880. doi: 10.3934/era.2022094
    [3] Yan Li, Yaqiang Wang . Some new results for $ B_1 $-matrices. Electronic Research Archive, 2023, 31(8): 4773-4787. doi: 10.3934/era.2023244
    [4] Yongge Tian . Characterizations of matrix equalities involving the sums and products of multiple matrices and their generalized inverse. Electronic Research Archive, 2023, 31(9): 5866-5893. doi: 10.3934/era.2023298
    [5] Zoran Pucanović, Marko Pešović . Analyzing Chebyshev polynomial-based geometric circulant matrices. Electronic Research Archive, 2024, 32(9): 5478-5495. doi: 10.3934/era.2024254
    [6] Jiachen Mu, Duanzhi Zhang . Multiplicity of symmetric brake orbits of asymptotically linear symmetric reversible Hamiltonian systems. Electronic Research Archive, 2022, 30(7): 2417-2427. doi: 10.3934/era.2022123
    [7] Xiuyun Guo, Xue Zhang . Determinants and invertibility of circulant matrices. Electronic Research Archive, 2024, 32(7): 4741-4752. doi: 10.3934/era.2024216
    [8] Yifan Luo, Kaisheng Lei, Qingzhong Ji . On the sumsets of units in a ring of matrices over $ \mathbb{Z}/m\mathbb{Z} $. Electronic Research Archive, 2025, 33(3): 1323-1332. doi: 10.3934/era.2025059
    [9] Jun-e Feng, Rong Zhao, Yanjun Cui . Simplification of logical functions with application to circuits. Electronic Research Archive, 2022, 30(9): 3320-3336. doi: 10.3934/era.2022168
    [10] Daochang Zhang, Dijana Mosić, Liangyun Chen . On the Drazin inverse of anti-triangular block matrices. Electronic Research Archive, 2022, 30(7): 2428-2445. doi: 10.3934/era.2022124
  • Natural and household discharges are the natural breeding grounds of various mosquito species, including female Anopheles mosquitoes, which transmit the Plasmodium parasite, causing the spread of the life-threatening disease malaria. Apart from that, population migrations also have a substantial impact on malaria transmission, claiming about half a million lives every year around the world. To assess the effects of the cumulative density of households and other natural discharges, and emigration-dependent interaction rates on the dissemination of the vector-borne infectious disease malaria, we propose and analyze a non-linear mathematical model. The model comprises five dependent variables, namely, the density of the susceptible human population, the density of the infective human population, the density of the susceptible female Anopheles mosquito population, the density of the infective mosquito population and cumulative density of household and other natural discharges. In the model, the density of the mosquito population is supposed to follow logistic growth, whose intrinsic growth rate is a linear function of the cumulative density of household and other natural discharges. The nonlinear model is analyzed by using the stability theory of differential equations, numerical simulations and sensitivity analysis. The analysis shows that an increase in non-emigrating population causes increased incidence of malaria. It is also found that an increase in household and other natural discharges accelerates the occurrence of malaria. A basic differential sensitivity analysis is carried out to assess the sensitivity of model solutions with respect to key parameters. The model's numerical simulations demonstrate the analytical findings.



    In 1999, the concept of k-subdirect sums of square matrices was proposed by Fallat and Johnson [1], which is a generalization of the usual sum of matrices [2]. The subdirect sum of matrices plays an important role in many areas, such as matrix completion problems, global stiffness matrices in finite elements and overlapping subdomains in domain decomposition methods [1,2,3,4,5].

    An important question for subdirect sums is whether the k-subdirect sum of two square matrices in one class of matrices lies in the same class. This question has attracted widespread attention in different classes of matrices and produced a variety of results. In 2005, Bru et al. gave sufficient conditions ensuring that the subdirect sum of two nonsingular M-matrices was also a nonsingular M-matrix [3]. Then the following year, they further came to the conclusion of the k-subdirect sum of S-SDD matrices is also an S-SDD matrix [2]. In [6], Chen and Wang succeeded in producing some sufficient conditions that the k-subdirect sum of SDD1 matrices is an SDD1 matrix. In [7], Li et al. gave some sufficient conditions such that the k-subdirect sum of doubly strictly diagonally dominant (DSDD) matrices is in the class of DSDD matrices. In addition, the k-subdirect sum of other classes of matrices were mentioned, such as Nekrasov matrices [8,9,10], quasi-Nekrasov (QN) matrices [11], SDD(p) matrices [12], weakly chained diagonally dominant matrices [13], Ostrowski-Brauer Sparse (OBS) matrices [14], {i0}-Nekrasov matrices [15], {p1,p2}-Nekrasov matrices [16], Dashnic-Zusmanovich (DZ) matrices [17], and B-matrices [18,19].

    GSDD1 matrices as a new subclass of H-matrices was proposed by Dai et al. in 2023 [20]. In this paper, we focus on the subdirect sum of GSDD1 matrices, and some sufficient conditions such that the k-subdirect sum of GSDD1 matrices with SDD matrices belong to GSDD1 matrices are given. Numerical examples are presented to illustrate the corresponding results.

    Now, some definitions are listed as follows.

    Definition 1.1. ([2]) Let A and B be two square matrices of order n1 and n2, respectively, and k be an integer such that 1kmin{n1,n2}, and let A and B be partitioned into 2×2 blocks as follows:

    A=(A11A12A21A22),B=(B11B12B21B22), (1.1)

    where A22 and B11 are square matrices of order k. Following [1], we call the square matrix of order n=n1+n2k given by

    C=(A11A120A21A22+B11B120B21B22)

    the k-subdirect sum of A and B, denoted by C=AkB. We can use the elements in A and B to represent any element in C. Before that, let us define the following set of indices:

    S1={1,2,...,n1k},S2={n1k+1,n1k+2,...,n1},S3={n1+1,...,n}. (1.2)

    Obviously, S1S2S3=N:={1,2,...,n}. Denoting C=AkB=[cij], A=[aij] and B=[bij], then

    cij={aij,iS1,jS1S2,0,iS1,jS3,aij,iS2,jS1,aij+bin1+k,jn1+k,iS2,jS2,bin1+k,jn1+k,iS2,jS3,0,iS3,jS1,bin1+k,jn1+k,iS3,jS2S3.

    Definition 1.2. ([20]) Given a matrix A=[aij]Cn×n, where Cn×n is the set of complex matrices. Let

    ri(A)=jN,ji|aij|,iN.
    NA={i||aii|ri(A)},
    ¯NA={i||aii|>ri(A)}.

    It is easy to obtain that ¯NA is the complement of NA in N, i.e., ¯NA=NNA.

    Definition 1.3. ([6]) A matrix A=[aij]Cn×n is called a strictly diagonally dominant (SDD) matrix if

    |aii|>ri(A),iN.

    Definition 1.4. ([20]) A matrix A=[aij]Cn×n is called a GSDD1 matrix if

    {ri(A)>pi¯NA(A),i¯NA,(ri(A)pi¯NA(A))(|ajj|pjNA(A))>piNA(A)pj¯NA(A),i¯NA,jNA,

    where

    piNA(A):=jNA{i}|aij|,pi¯NA(A):=j¯NA{i}rj(A)|ajj||aij|,iN.

    Remark 1.1. From Definitions 1.3 and 1.4, it is easy to obtain that if a matrix A is an SDD matrix with ri(A)>0, then it is a GSDD1 matrix.

    First of all, a counterexample is given to show that the subdirect sum of two GSDD1 matrices may not necessarily be a GSDD1 matrix.

    Example 2.1. Consider the following GSDD1 matrices A and B, where

    A=(432143013.5),B=(2.5201212.31.84).

    and the 1-subdirect sum C=A1B is

    C=(43200143000162000121002.31.84).

    However, C is not a GSDD1 matrix because

    (r3(C)p3¯NC(C))(|c11|p1NC(C))=(30)(43)=3=3×1=p3NC(C)p1¯NC(C).

    Example 2.1 shows that the subdirect sum of GSDD1 matrices is not a GSDD1 matrix. Then, a meaningful discussion is concerned with: under what conditions will the subdirect sum of GSDD1 matrices is in the class of GSDD1 matrices?

    In order to obtain the main results, several lemmas are introduced that will be used in the sequel.

    Lemma 2.1. If matrix A=[aij]Cn×n is a GSDD1 matrix, then |ajj|pjNA(A)>0 holds for all jNA.

    Proof. According to the definition of GSDD1 matrices, we get

    (ri(A)pi¯NA(A))(|ajj|pjNA(A))>piNA(A)pj¯NA(A).

    Since ri(A)pi¯NA(A)>0,piNA(A), and pj¯NA(A) are all nonnegative, |ajj|pjNA(A)>0 is obtained.

    Lemma 2.2. ([20]) If A=[aij]Cn×n is a GSDD1 matrix, then there is at least one entry aij0, ij, i¯NA, jN.

    Lemma 2.3. ([20]) If A=[aij]Cn×n is a GSDD1 matrix with NA=, then A is an SDD matrix, and there is at least one entry aij0, ij, i¯NA, j¯NA.

    Now, we consider the 1-subdirect sum of GSDD1 matrices.

    Theorem 2.1. Let A=[aij] and B=[bij] be square matrices of order n1 and n2 partitioned as in (1.1), respectively. And let k=1, S1={1,2,,n11}, S2={n1}, and S3={n1+1,n1+2,,n1+n21}. We assume that A is a GSDD1 matrix, and B is an SDD matrix with ri(B)>0 for all i¯NB. If all diagonal entries of A22 and B11 are positive (or all negative), n1¯NA and

    rn1(A)|an1,n1|rn1(A)+r1(B)|an1,n1+b11|,

    then the 1-subdirect sum C=A1B is a GSDD1 matrix.

    Proof. According to the 1-subdirect sum C=A1B, we have

    rn1(C)=rn1(A)+r1(B).

    From n1¯NA, we know |an1,n1|>rn1(A). Because all diagonal entries of A22 and B11 are positive (or negative), we have

    |cn1,n1|=|an1,n1+b11|=|an1,n1|+|b11|>rn1(A)+r1(B)=rn1(C).

    Since A is a GSDD1 matrix, B is an SDD matrix with ri(B)>0 for all i¯NB, C=A1B, and according to Lemmas 2.2 and 2.3, we know that ri(C)0 for all i¯NC. Therefore, for any i¯NC,

    ri(C)=jN{i}|cij|>j¯NC{i}rj(C)|cjj||cij|=pi¯NC(C).

    For any jNC, we easily get jNCS1=NAS1S1. For the three different selection ranges of i, that is, i¯NCS1=¯NAS1S1, i¯NCS2={n1}, and i¯NCS3S3, therefore, we divide the proof into three cases.

    Case 1. For i¯NCS1=¯NAS1S1, jNC, we have

    ri(C)=ri(A),
    pi¯NC(C)=j¯NC{i},jS1rj(C)|cjj||cij|+j¯NC{i},jS2rj(C)|cjj||cij|+j¯NC{i},jS3rj(C)|cjj||cij|=j¯NA{i},jS1rj(A)|ajj||aij|+rn1(A)+r1(B)|an1,n1+b11||ai,n1|+0j¯NA{i},jS1rj(A)|ajj||aij|+rn1(A)|an1,n1||ai,n1|=pi¯NA(A),
    |cjj|=|ajj|, (2.1)
    pjNC(C)=jNC{j}|cjj|=jNA{j}|ajj|=pjNA(A), (2.2)
    piNC(C)=jNC{i}|cij|=jNA{i}|aij|=piNA(A),
    pj¯NC(C)=j¯NC{j},jS1rj(C)|cjj||cjj|+j¯NC{j},jS2rj(C)|cjj||cjj|+j¯NC{j},jS3rj(C)|cjj||cjj|=j¯NA{j},jS1rj(A)|ajj||ajj|+rn1(A)+r1(B)|an1,n1+b11||aj,n1|+0j¯NA{j},jS1rj(A)|ajj||ajj|+rn1(A)|an1,n1||aj,n1|=pj¯NA(A). (2.3)

    Therefore, we obtain that

    (ri(C)pi¯NC(C))(|cjj|pjNC(C))(ri(A)pi¯NA(A))(|ajj|pjNA(A))>piNA(A)pj¯NA(A)piNC(C)pj¯NC(C).

    Case 2. For i¯NCS2={n1}, jNC,

    rn1(C)=rn1(A)+r1(B),
    pn1NC(C)=jNC{n1}|cn1,j|=jNA{n1}|an1,j|=pn1NA(A).
    pn1¯NC(C)=j¯NC{n1},jS1rj(C)|cjj||cn1,j|+j¯NC{n1},jS3rj(C)|cjj||cn1,j|=j¯NA{n1}rj(A)|ajj||an1,j|+j¯NB{1}rj(B)|bjj||b1j|=pn1¯NA(A)+p1¯NB(B).

    We know that the results of the |cjj|, pjNC(C), and pj¯NC(C) are the same as (2.1), (2.2), and (2.3). Because B is an SDD matrix with ri(B)>0 for all i¯NB, we clearly get

    r1(B)p1¯NB(B)>0.

    Hence,

    (rn1(C)pn1¯NC(C))(|cjj|pjNC(C))=(rn1(A)+r1(B)pn1¯NA(A)p1¯NB(B))(|ajj|pjNA(A))>(rn1(A)pn1¯NA(A))(|ajj|pjNA(A))>pn1NA(A)pj¯NA(A)pn1NC(C)pj¯NC(C).

    Case 3. For i¯NCS3S3, jNC, in particular, we obtain that

    piNC(C)=jNC{i}|cij|=0.

    So we easily come up with

    (rn1(C)pn1¯NC(C))(|cjj|pjNC(C))=(rn1(C)pn1¯NC(C))(|ajj|pjNA(A))>0=piNC(C)pj¯NC(C).

    From Cases 1–3, we have that for any i¯NC and jNC, the C matrix satisfies the definition of the GSDD1 matrix. The conclusion is as follows.

    Theorem 2.2. Let A=[aij] and B=[bij] be square matrices of order n1 and n2 partitioned as in (1.1), respectively. And let k, S1, S2, and S3 be as in Theorem 2.1. Likewise, we assume A is a GSDD1 matrix, and B is an SDD matrix with ri(B)>0 for all i¯NB. If all diagonal entries of A22 and B11 are positive (or all negative), n1NA, rn1(A)+r1(B)|an1,n1|+|b11| and

    min2ln2(rl(B)pl¯NB(B))maxm¯NA(rm(A)pm¯NA(A)),
    minm¯NApmNA(A)max2ln2|bl1|,

    then C=A1B is a GSDD1 matrix.

    Proof. Since A is a GSDD1 matrix, B is an SDD matrix with ri(B)>0 for all i¯NB, n1NA, and rn1(A)+r1(B)|an1,n1|+|b11|, we get n1NC and then NC=NA.

    For any i¯NC, by Lemmas 2.2 and 2.3, we have ri(C)0 and then

    ri(C)=jN{i}|cij|>j¯NC{i}rj(C)|cjj||cij|=pi¯NC(C).

    Since n1NC, i.e., i¯NCS2=, we prove it according to the two different selection ranges of i, namely i¯NCS1=¯NAS1S1 and i¯NCS3S3. For any jNC, that is, jNCS1=NAS1S1 and jNCS2=NAS2={n1}. Therefore, we prove it from the following cases.

    Case 1. For i¯NCS1=¯NAS1S1, jNCS1=NAS1S1, we obtain that

    ri(C)=ri(A), (2.4)
    pi¯NC(C)=j¯NC{i},jS1rj(C)|cjj||cij|+j¯NC{i},jS3rj(C)|cjj||cij|=j¯NA{i}rj(A)|ajj||aij|+0=pi¯NA(A), (2.5)
    |cjj|=|ajj|, (2.6)
    pjNC(C)=jNC{j}|cjj|=jNA{j}|ajj|=pjNA(A), (2.7)
    piNC(C)=jNC{i}|cij|=jNA{i}|aij|=piNA(A), (2.8)
    pj¯NC(C)=j¯NC{j},jS1rj(C)|cjj||cjj|+j¯NC{j},jS3rj(C)|cjj||cjj|=j¯NA{j}rj(A)|ajj||ajj|+0=pj¯NA(A). (2.9)

    Therefore,

    (ri(C)pi¯NC(C))(|cjj|pjNC(C))=(ri(A)pi¯NA(A))(|ajj|pjNA(A))>piNA(A)pj¯NA(A)=piNC(C)pj¯NC(C).

    Case 2. For i¯NCS1=¯NAS1S1, jNCS2=NAS2={n1}, we know that ri(C), piNC(C), and pi¯NC(C) have the same results as (2.4), (2.8), and (2.5). Moreover,

    |cn1,n1|=|an1,n1+b11|=|an1,n1|+|b11|, (2.10)
    pn1NC(C)=jNC{n1}|cn1,j|=jNA{n1}|an1,j|=pn1NA(A), (2.11)
    pn1¯NC(C)=j¯NC{n1},jS1rj(C)|cjj||cn1,j|+j¯NC{n1},jS3rj(C)|cjj||cn1,j|=j¯NA{n1}rj(A)|ajj||an1,j|+j¯NB{1}rj(B)|bjj||b1j|=pn1¯NA(A)+p1¯NB(B). (2.12)

    Hence, we obtain that

    (ri(C)pi¯NC(C))(|cn1,n1|pn1NC(C))=(ri(A)pi¯NA(A))(|an1,n1|+|b11|pn1NA(A))=(ri(A)pi¯NA(A))(|an1,n1|pn1NA(A))+(ri(A)pi¯NA(A))|b11|>piNA(A)pn1¯NA(A)+piNA(A)p1¯NB(B)=piNC(C)pn1¯NC(C).

    Case 3. For i¯NCS3S3, jNCS1=NAS1S1, we have

    ri(C)=rin1+1(B)=rl(B), (2.13)
    pi¯NC(C)=j¯NC{i},jS1rj(C)|cjj||cij|+j¯NC{i},jS3rj(C)|cjj||cij|=0+j¯NB{l},j{2,,n2}rj(B)|bjj||blj|pl¯NB(B), (2.14)
    piNC(C)=jNC{i},jS1|cij|+jNC{i},j=n1|ci,n1|=0+|bl1|=|bl1|, (2.15)

    where l=in1+1. We have the same values of |cjj|, pjNC(C), and pj¯NC(C) as (2.6), (2.7), and (2.9). Therefore,

    (ri(C)pi¯NC(C))(|cjj|pjNC(C))(rl(B)pl¯NB(B))(|ajj|pjNA(A))(rm(A)pm¯NA(A))(|ajj|pjNA(A))>pmNA(A)pj¯NA(A)|bl1|pj¯NA(A)=piNC(C)pj¯NC(C).

    Case 4. For i¯NCS3S3, jNCS2=NAS2={n1}, we get that the values of ri(C), piNC(C), and pi¯NC(C) are the same as (2.13), (2.15), and (2.14). Moreover, the results of |cjj|, pjNC(C), and pj¯NC(C) are the same as (2.10), (2.11), and (2.12). Hence, we obtain that

    (ri(C)pi¯NC(C))(|cn1,n1|pn1NC(C))(rl(B)pl¯NB(B))(|an1,n1|+|b11|pn1NA(A))=(rl(B)pl¯NB(B))|b11|+(rl(B)pl¯NB(B))(|an1,n1|pn1NA(A))(rm(A)pm¯NA(A))|b11|+(rm(A)pm¯NA(A))(|an1,n1|pn1NA(A))>pmNA(A)p1¯NB(B)+pmNA(A)pn1¯NA(A)|bl1|(p1¯NB(B)+pn1¯NA(A))=piNC(C)pn1¯NC(C).

    From Cases 1–4, we definitively get that C is a GSDD1 matrix.

    The following Example 2.2 shows that Theorem 2.1 may not necessarily hold when k2.

    Example 2.2. Consider the following matrices:

    A=(311.71141122410113),B=(311021103),

    where A is a GSDD1 matrix and B is an SDD matrix with ri(B)>0 for all i¯NB. It is easy to verify that A and B satisfy the conditions of Theorem 2.1 and A1B is a GSDD1 matrix. However, C=A2B is not a GSDD1 matrix. In fact,

    C=(311.71014110227210115100103).

    By computation,

    ¯NC={2,4,5},NC={1,3},
    (r5(C)p5¯NC(C))(|c11|p1NC(C))=(10)(31.7)=1.3<1.35=1×1.35=p5NC(C)p1¯NC(C).

    Therefore, C=A2B is not a GSDD1 matrix.

    The following Example 2.3 shows that Theorem 2.2 may not necessarily hold when k2.

    Example 2.3. Consider the following matrices:

    A=(5221040100321112),B=(3201154.30.95.117),

    where A is a GSDD1 matrix and B is an SDD matrix with ri(B)>0 for all i¯NB. It is easy to verify that A and B satisfy the conditions of Theorem 2.2 and A1B is a GSDD1 matrix. However, C=A2B is not a GSDD1. In fact,

    C=(522100401000600112174.3000.95.117).

    By computation, r3(C)p3¯NC(C)=0, therefore, C=A2B is not a GSDD1 matrix.

    Those are sufficient conditions to ensure that the 1-subdirect sum of GSDD1 matrices with SDD matrices is a GSDD1 matrix. In fact, as the value of k increases, the situation becomes more complicated, so that the adequate conditions we give will also be more complicated.

    Next, some sufficient conditions ensuring that the k-subdirect (k2) sum of GSDD1 matrices with SDD matrices is a GSDD1 matrix are given.

    Theorem 2.3. Let A=[aij] and B=[bij] be square matrices of order n1 and n2 partitioned as in (1.1), respectively. And let 2kmin{n1,n2}, S1, S2, and S3 be as in (1.2). We assume A is a GSDD1 matrix and B is an SDD matrix with ri(B)>0 for all i¯NB. If all diagonal entries of A22 and B11 are positive (or all negative), i¯NA for any iS2 and

    j¯NA{i},jS2λj|ajj+bjn1+k,jn1+k||aij|j¯NA{i},jS2rj(A)|ajj||aij|,(iS1S2)
    j¯NB{in1+k}j{1,,k}λj+n1k|aj+n1k,j+n1k+bjj||bin1+k,j|j¯NB{in1+k}j{1,,k}rj(B)|bjj||bin1+k,j|,(iS2)
    λiri(A)+pin1+k¯NB(B),(iS2)

    where λi=ri(A)+rin1+k(B)+n1j=n1k+1ji|aij+bin1+k,jn1+k|n1j=n1k+1ji(|aij|+|bin1+k,jn1+k|), then the k-subdirect sum C=AkB is a GSDD1 matrix.

    Proof. Since A is a GSDD1 matrix with i¯NA for any iS2, we get |aii|>ri(A). According to the k-subdirect sum C=AkB, we have ri(C)=λiri(A)+rin1+k(B). Because all diagonal entries of A22 and B11 are positive (or negative), we get |cii|=|aii|+|bin1+k,in1+k|. Therefore, we obtain that |cii|>ri(C), that is, for any iS2, i¯NC. Since A is a GSDD1 matrix, B is an SDD matrix with ri(B)>0 for all i¯NB, and C=AkB, by Lemmas 2.2 and 2.3 we know that ri(C)0 for i¯NCS1S3=¯NAS1S3. For iS2, by sufficient conditions, we have λiri(A)+pin1+k¯NB(B), which means that λi>0. Therefore, for any i¯NC, we obtain that

    ri(C)=jN{i}|cij|>j¯NC{i}rj(C)|cjj||cij|=pi¯NC(C).

    Moreover, for any jNC, we get jNCS1=NAS1S1. For any i¯NC, similarly, we prove it from the following three cases, which are i¯NCS1=¯NAS1S1, i¯NCS2=¯NAS2S2, and i¯NCS3S3.

    Case 1. For i¯NCS1=¯NAS1S1, jNC, we have

    ri(C)=ri(A),
    pi¯NC(C)=j¯NC{i},jS1rj(C)|cjj||cij|+j¯NC{i},jS2rj(C)|cjj||cij|+j¯NC{i},jS3rj(C)|cjj||cij|=j¯NA{i},jS1rj(A)|ajj||aij|+j¯NA{i},jS2λj|ajj+bjn1+k,jn1+k||aij|+0j¯NA{i},jS1rj(A)|ajj||aij|+j¯NA{i},jS2rj(A)|ajj||aij|=pi¯NA(A),
    |cjj|=|ajj|, (2.16)
    pjNC(C)=jNC{j}|cjj|=jNA{j}|ajj|=pjNA(A), (2.17)
    piNC(C)=jNC{i}|cij|=jNA{i}|aij|=piNA(A),
    pj¯NC(C)=j¯NC{j},jS1rj(C)|cjj||cjj|+j¯NC{j},jS2rj(C)|cjj||cjj|+j¯NC{j},jS3rj(C)|cjj||cjj|=j¯NA{j},jS1rj(A)|ajj||ajj|+j¯NA{j},jS2λj|ajj+bjn1+k,jn1+k||ajj|+0j¯NA{j},jS1rj(A)|ajj||ajj|+j¯NA{j},jS2rj(A)|ajj||ajj|=pj¯NA(A). (2.18)

    Therefore,

    (ri(C)pi¯NC(C))(|cjj|pjNC(C))(ri(A)pi¯NA(A))(|ajj|pjNA(A))>piNA(A)pj¯NA(A)piNC(C)pj¯NC(C).

    Case 2. For i¯NCS2=¯NAS2S2, jNC, we obtain that

    ri(C)=λi,
    pi¯NC(C)=j¯NC{i},jS1rj(C)|cjj||cij|+j¯NC{i},jS2rj(C)|cjj||cij|+j¯NC{i},jS3rj(C)|cjj||cij|=j¯NA{i},jS1rj(A)|ajj||aij|+j¯NA{i},jS2λj|ajj+bjn1+k,jn1+k||aij+bin1+k,jn1+k|+j¯NB{in1+k}j{k+1,,n2}rj(B)|bjj||bin1+k,j|j¯NA{i},jS1rj(A)|ajj||aij|+j¯NA{i},jS2λj|ajj+bjn1+k,jn1+k||aij|+j¯NA{i},jS2λj|ajj+bjn1+k,jn1+k||bin1+k,jn1+k|+j¯NB{in1+k}j{k+1,,n2}rj(B)|bjj||bin1+k,j|j¯NA{i},jS1rj(A)|ajj||aij|+j¯NA{i},jS2rj(A)|ajj||aij|+j¯NB{in1+k}j{1,,k}λj+n1k|aj+n1k,j+n1k+bjj||bin1+k,j|+j¯NB{in1+k}j{k+1,,n2}rj(B)|bjj||bin1+k,j|pi¯NA(A)+j¯NB{in1+k}j{1,,k}rj(B)|bjj||bin1+k,j|+j¯NB{in1+k}j{k+1,,n2}rj(B)|bjj||bin1+k,j|=pi¯NA(A)+pin1+k¯NB(B),
    piNC(C)=jNC{i}|cij|=jNA{i}|aij|=piNA(A).

    We know that |cjj|, pjNC(C), and pj¯NC(C) are the same as (2.16), (2.17), and (2.18). Therefore,

    (ri(C)pi¯NC(C))(|cjj|pjNC(C))(λipi¯NA(A)pin1+k¯NB(B))(|ajj|pjNA(A))ri(A)+pin1+k¯NB(B)pi¯NA(A)pin1+k¯NB(B)×(|ajj|pjNA(A))=(ri(A)pi¯NA(A))(|ajj|pjNA(A))>piNA(A)pj¯NA(A)piNC(C)pj¯NC(C).

    Case 3. For i¯NCS3S3, jNC, specifically, we obtain that

    pNCi(C)=jNC/{i}|cij|=0.

    Hence,

    (ri(C)pi¯NC(C))(|cjj|pjNC(C))=(ri(C)pi¯NC(C))(|ajj|pjNA(A))>0=piNC(C)pj¯NC(C).

    Therefore, we get that ri(C)pi¯NC(C)>0 and (ri(C)pi¯NC(C))(|cjj|pjNC(C))>piNC(C)pj¯NC(C) for any i¯NC, jNC.

    Corollary 2.1. Let A=[aij] and B=[bij] be square matrices of order n1 and n2 partitioned as in (1.1), respectively. And let 2kmin{n1,n2}, S1, S2, and S3 be as in (1.2). We assume A is a GSDD1 matrix and B is an SDD matrix with ri(B)>0 for all i¯NB. If all diagonal entries of A22 and B11 are positive (or all negative), i¯NA for any iS2 and

    λj|ajj+bjn1+k,jn1+k|min{rj(A)|ajj|,rjn1+k(B)|bjn1+k,jn1+k|},(jS2)
    λiri(A)+pin1+k¯NB(B),

    where λi is the same as λi of Theorem 2.3 and iS2, then the k-subdirect sum C=AkB is a GSDD1 matrix.

    Proof. For the inequality

    λi|aii+bin1+k,in1+k|ri(A)|aii|,

    multiplying both sides of this inequality by |aij| (iS1S2,ji) and summing for every j¯NA{i} (jS2), we have

    j¯NA{i},jS2λj|ajj+bjn1+k,jn1+k||aij|j¯NA{i},jS2rj(A)|ajj||aij|.

    Similarly, for iS2, we obtain that

    j¯NB{in1+k}j{1,,k}λj+n1k|aj+n1k,j+n1k+bjj||bin1+k,j|j¯NB{in1+k}j{1,,k}rj(B)|bjj||bin1+k,j|.

    By Theorem 2.3, we obtain that the k-subdirect sum C=AkB is a GSDD1 matrix.

    Theorem 2.4. Let A=[aij] and B=[bij] be square matrices of order n1 and n2 partitioned as in (1.1), respectively. And let 2kmin{n1,n2}, S1, S2, and S3 be as in (1.2). We assume A is a GSDD1 matrix and B is an SDD matrix with ri(B)>0 for all i¯NB. If all diagonal entries of A22 and B11 are positive (or all negative), iNA for any iS2, |aii|+|bin1+k,in1+k|λi and

    |bjn1+k,jn1+k|j¯NB{jn1+k}j{1,,k}|aj,j+n1k+bjn1+k,j|pjn1+k¯NB(B),(jS2)
    mink+1ln2(rl(B)pl¯NB(B))maxm¯NA(rm(A)pm¯NA(A)),
    minm¯NApmNA(A)maxk+1ln2j{1,,k}|blj|,

    where λi is the same as λi of Theorem 2.3, then the k-subdirect sum C=AkB is a GSDD1 matrix.

    Proof. Since A is a GSDD1 matrix with iNA for any iS2 and |aii|+|bin1+k,in1+k|λi, we have |aii|ri(A) and |cii|=|aii|+|bin1+k,in1+k|λi=ri(C), that is, for any iS2, we have iNC. Moreover, we know that iNCS1=NAS1S1, which means that NC=NA. Combining Lemmas 2.2 and 2.3, we get that ri(C)0 for i¯NCS1S3=¯NAS1S3. Therefore, for any i¯NC, we obtain that

    ri(C)=jN{i}|cij|>j¯NC{i}rj(C)|cjj||cij|=pi¯NC(C).

    Since i¯NCS2=, we prove it from the following two aspects, which are i¯NCS1=¯NAS1S1 and i¯NCS3S3. For any jNC, that is, jNCS1=NAS1S1 and jNCS2=NAS2S2. Therefore, we prove it from the following cases.

    Case 1. For i¯NCS1=¯NAS1S1, jNCS1=NAS1S1, we get

    ri(C)=ri(A), (2.19)
    pi¯NC(C)=j¯NC{i},jS1rj(C)|cjj||cij|+j¯NC{i},jS3rj(C)|cjj||cij|=j¯NA{i}rj(A)|ajj||aij|+0=pi¯NA(A), (2.20)
    |cjj|=|ajj|, (2.21)
    pjNC(C)=jNC{j}|cjj|=jNA{j}|ajj|=pjNA(A), (2.22)
    piNC(C)=jNC{i}|cij|=jNA{i}|aij|=piNA(A), (2.23)
    pj¯NC(C)=j¯NC{j},jS1rj(C)|cjj||cjj|+j¯NC{j},jS3rj(C)|cjj||cjj|=j¯NA{j}rj(A)|ajj||ajj|+0=pj¯NA(A). (2.24)

    Therefore, we obtain that

    (ri(C)pi¯NC(C))(|cjj|pjNC(C))=(ri(A)pi¯NA(A))(|ajj|pjNA(A))>piNA(A)pj¯NA(A)=piNC(C)pj¯NC(C).

    Case 2. For i¯NCS1=¯NAS1S1, jNCS2=NAS2S2, we know that ri(C) and (2.19)} are equal, pi¯NC(C) and (2.20) are equal, and piNC(C) and (2.23) are equal. Moreover,

    |cjj|=|ajj+bjn1+k,jn1+k|=|ajj|+|bjn1+k,jn1+k|, (2.25)
    pjNC(C)=jNC{j},jS1|cjj|+jNC{j},jS2|cjj|=jNA{j},jS1|ajj|+j¯NB{jn1+k}j{1,,k}|aj,j+n1k+bjn1+k,j|, (2.26)
    pj¯NC(C)=j¯NC{j},jS1rj(C)|cjj||cjj|+j¯NC{j},jS3rj(C)|cjj||cjj|=j¯NA{j}rj(A)|ajj||ajj|+j{k+1,,n2}rj(B)|bjj||bjn1+k,j|pj¯NA(A)+pjn1+k¯NB(B). (2.27)

    Hence,

    (ri(C)pi¯NC(C))(|cjj|pjNC(C))=(ri(A)pi¯NA(A))×(|ajj|+|bjn1+k,jn1+k|jNA{j},jS1|ajj|j¯NB{jn1+k}j{1,,k}|aj,j+n1k+bjn1+k,j|)=(ri(A)pi¯NA(A))(|ajj|jNA{j},jS1|ajj|)+(ri(A)pi¯NA(A))(|bjn1+k,jn1+k|j¯NB{jn1+k}j{1,,k}|aj,j+n1k+bjn1+k,j|)(ri(A)pi¯NA(A))(|ajj|pjNA(A))+piNA(A)pjn1+k¯NB(B)>piNA(A)pj¯NA(A)+piNA(A)pjn1+k¯NB(B)piNC(C)pj¯NC(C).

    Case 3. For i¯NCS3S3, jNCS1=NAS1S1, we obtain that

    ri(C)=rin1+k(B)=rl(B), (2.28)
    pi¯NC(C)=j¯NC{i},jS1rj(C)|cjj||cij|+j¯NC{i},jS3rj(C)|cjj||cij|=0+j¯NB{l}j{k+1,,n2}rj(B)|bjj||blj|pl¯NB(B), (2.29)
    piNC(C)=jNC{i},jS1|cij|+jNC{i},jS2|cij|=0+j{1,,k}|bin1+k,j|=j{1,,k}|blj|, (2.30)

    where l=in1+k. We know that |cjj|, pjNC(C), and pj¯NC(C) are the same as (2.21), (2.22), and (2.24). Therefore,

    (ri(C)pi¯NC(C))(|cjj|pjNC(C))(rl(B)pl¯NB(B))(|ajj|pjNA(A))(rm(A)pm¯NA(A))(|ajj|pjNA(A))>pmNA(A)pj¯NA(A)j{1,,k}|blj|pj¯NA(A)=piNC(C)pj¯NC(C).

    Case 4. For i¯NCS3S3, jNCS2=NAS2S2, we obtain that the values of ri(C), piNC(C), and pi¯NC(C) are equal to (2.28), (2.30), and (2.29). Moreover, the results of |cjj|, pjNC(C), and pj¯NC(C) are the same as (2.25), (2.26), and (2.27). Hence, we arrive at

    (ri(C)pi¯NC(C))(|cjj|pjNC(C))(rl(B)pl¯NB(B))×(|ajj|+|bjn1+k,jn1+k|(jNA{j},jS1|ajj|+j¯NB{jn1+k}j{1,,k}|aj,j+n1k+bjn1+k,j|))=(rl(B)pl¯NB(B))(|ajj|jNA{j},jS1|ajj|)+(rl(B)pl¯NB(B))(|bjn1+k,jn1+k|j¯NB{jn1+k}j{1,,k}|aj,j+n1k+bjn1+k,j|)(rm(A)pm¯NA(A))(|ajj|pjNA(A))+(rm(A)pm¯NA(A))pjn1+k¯NB(B)>pmNA(A)pj¯NA(A)+pmNA(A)pjn1+k¯NB(B)=pmNA(A)(pj¯NA(A)+pjn1+k¯NB(B))j{1,,k}|blj|(pj¯NA(A)+pjn1+k¯NB(B))piNC(C)pj¯NC(C).

    In conclusion, for any i¯NC, jNC, we successfully derive that ri(C)pi¯NC(C)>0 and} (ri(C)pi¯NC(C))(|cjj|pjNC(C))>piNC(C)pj¯NC(C). Therefore, C=AkB is a GSDD1 matrix.

    Example 2.4. Consider the following matrices:

    A=(7.512212.5170.3120.21.11.3510.810.410.26.51.20.90.310.20.96.61.40.70.90.11.218),B=(651.5211.51.2662.31.60.91.42671.31.233.42660.60.42.111.877),

    where A is a GSDD1 matrix with i¯NA for all iS2, and B is an SDD matrix with ri(B)>0 for all i¯NB. By computation, we derive NA={1,3},¯NA={2,4,5,6}. Moreover,

    λ4|a44+b11|=5.571.50.077<0.5693.76.5=r4(A)|a44|,λ5|a55+b22|=5.272.60.072<0.5763.86.6=r5(A)|a55|,
    λ6|a66+b33|=5.475=0.072<0.4883.98=r6(A)|a66|,

    we get that j¯NA{i},jS2λj|ajj+bjn1+k,jn1+k||aij|j¯NA{i},jS2rj(A)|ajj||aij| is true for iS1S2.

    λ4|a44+b11|0.077<0.092665=r1(B)|b11|,λ5|a55+b22|0.072<0.091666=r2(B)|b22|,
    λ6|a66+b33|=0.072<0.0885.967=r3(B)|b33|,

    we have that the second sufficient condition in Theorem 2.3 is true.

    λ4=5.5>4.252=3.7+0.552r4(A)+p1¯NB(B),λ5=5.2>4.393=3.8+0.593r5(A)+p2¯NB(B),
    λ6=5.4>4.471=3.9+0.571r6(A)+p3¯NB(B),

    we get that the third sufficient condition in Theorem 2.3 is met. Therefore, by Theorem 2.3, C=A3B is a GSDD1 matrix. In fact,

    C=(7.512212.500170.3120.2001.11.3510.81000.410.271.50.31.111.50.310.20.372.60.91.60.90.70.90.10.21751.31.200033.42660.60000.42.111.877),

    where NC={1,3},¯NC={2,4,5,6,7,8}. By computation,

    r2(C)=4.5,p2¯NC(C)0.235,p2NC(C)=1.3;r4(C)=5.5,p4¯NC(C)0.983,p4NC(C)=0.6;
    r5(C)=5.2,p5¯NC(C)1.011,p5NC(C)=0.5;r6(C)=5.4,p6¯NC(C)0.925,p6NC(C)=0.8;
    r7(C)=9,p7¯NC(C)0.66,p7NC(C)=0;r8(C)=5.3,p8¯NC(C)0.499,p8NC(C)=0;
    |c11|=7.5,p1NC(C)=2,p1¯NC(C)1.048;|c33|=5,p3NC(C)=1.1,p3¯NC(C)1.042.

    It is not difficult to find that ri(C)pi¯NC(C)>piNC(C) and |cjj|pjNC(C)>pj¯NC(C) when i¯NC, jNC. So we deduce that ri(C)>pi¯NC(C) and (ri(C)pi¯NC(C))(|cjj|pjNC(C))>piNC(C)pj¯NC(C) are true when i¯NC, jNC. Thus, C=A3B is a GSDD1 matrix.

    Example 2.5. Consider the following matrices:

    A=(620.5110.81.20.180.70.311.30.80.50.87.71.11.20.30.12.11.50.981.80.61.70.30.71.418.42.52.81.62.5211.79.21.10.81.21.62.41.81.59),B=(40210.71520.81.24532.51919.22.52.1541.12126.91.82.40.9612530.80.511.30.265101.40.40.70.59.968),

    where A is a GSDD1 matrix and B is an SDD matrix with ri(B)>0 for all i¯NB. By computation, NA={1,4,5,6,7},S2={4,5,6,7},

    |a44|+|b11|=48<48.1=λ4,|a55|+|b22|=53.4<53.6=λ5,
    |a66|+|b33|=63.2<63.5=λ6,|a77|+|b44|=70<70.2=λ7.

    Moreover,

    |b11|j¯NB{1},j{1,...,4}|a4,j+3+b1j|=32.2>10.633p1¯NB(B),
    |b22|j¯NB{2},j{1,...,4}|a5,j+3+b2j|=32>14.101p2¯NB(B),
    |b33|j¯NB{3},j{1,...,4}|a6,j+3+b3j|=44.5>14.965p3¯NB(B),
    |b44|j¯NB{4},j{1,...,4}|a7,j+3+b4j|=50.2>15.908p4¯NB(B).
    mink+1ln2(rl(B)pl¯NB(B))=r6(B)p6¯NB(B)7.944>3.836r2(A)p2¯NA(A)=maxm¯NA(rm(A)pm¯NA(A)),
    minm¯NApmNA(A)=p3NA(A)=3.2>3=j{1,...,4}|b5j|=j{1,...,4}|b6j|=maxk+1ln2j{1,,k}|blj|.

    Hence, the conditions in Theorem 2.4 are met. By Theorem 2.4, C=A4B is a GSDD1 matrix. In fact,

    C=(620.5110.81.2000.180.70.311.30.8000.50.87.71.11.20.30.1002.11.50.9483.81.62.41520.80.30.71.42.253.45.55.31919.21.62.523.53.863.22.22126.90.81.21.64.24.22.4702530.80000.511.30.265100001.40.40.70.59.968).

    By computation, NC={1,4,5,6,7}, ¯NC={2,3,8,9}. Moreover,

    r2(C)=4.2,p2¯NC(C)0.364,p2NC(C)=3.5;r3(C)=4,p3¯NC(C)=0.42,p3NC(C)=3.2;
    r8(C)=13,p8¯NC(C)1.897,p8NC(C)=3;r9(C)=12.9,p9¯NC(C)=1.98,p9NC(C)=3.
    |c11|=6,p1NC(C)=4,p1¯NC(C)1.31;|c44|=48,p4NC(C)=9.9,p4¯NC(C)8.201;
    |c55|=53.4,p5NC(C)=13.3,p5¯NC(C)8.537;|c66|=63.2,p6NC(C)=11.1,p6¯NC(C)11.655;
    |c77|=70,p7NC(C)=11.6,p7¯NC(C)12.304.

    We see that ri(C)pi¯NC(C)>piNC(C) and |cjj|pjNC(C)>pj¯NC(C) when i¯NC, jNC. Therefore, we obtain that ri(C)>pi¯NC(C) and (ri(C)pi¯NC(C))(|cjj|pjNC(C))>piNC(C)pj¯NC(C) are true when i¯NC, jNC. Therefore, C=A4B is a GSDD1 matrix.

    Remark 2.1. Since the subdirect sum of matrices does not satisfy the commutative law, if we change "A is a GSDD1 matrix, and B is an SDD matrix" to "A is an SDD matrix, and B is a GSDD1 matrix", then we will obtain new sufficient conditions by using similar proofs in this paper.

    In this paper, some sufficient conditions are given to show that the subdirect sum of GSDD1 matrices with SDD matrices is in the class of GSDD1 matrices, and these conditions are only dependent on the elements of the given matrices. Furthermore, some numerical examples are also presented to illustrate the corresponding theoretical results.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This work was partly supported by the National Natural Science Foundation of China (31600299), Natural Science Basic Research Program of Shaanxi, China (2020JM-622), and the Postgraduate Innovative Research Project of Baoji University of Arts and Sciences (YJSCX23YB33).

    The authors declare there is no conflicts of interest.



    [1] WHO Malaria: Fact sheets, 2022. Available from: https://www.who.int/news-room/fact-sheets/detail/malaria.
    [2] S. Sato, Plasmodium-a brief introduction to the parasites causing human malaria and their basic biology, J. Physiol. Anthropol., 40 (2021), 1–13. https://doi.org/10.1186/s40101-020-00251-9 doi: 10.1186/s40101-020-00251-9
    [3] E. Pondeville, N. Puchot, M. Lang, F. Cherrier, F. Schaffner, C. Dauphin-Villemant, et al., Evolution of sexually-transferred steroids and mating-induced phenotypes in Anopheles mosquitoes, Sci. Rep., 9 (2019), 4669. https://doi.org/10.1038/s41598-019-41094-4 doi: 10.1038/s41598-019-41094-4
    [4] N. T. J. Bailey, The Biomathematics of Malaria, Charles Griffin & Company Ltd., 1982.
    [5] L. Cai, Z. Li, C. Yang, J. Wang, Global analysis of an environmental disease transmission model linking within-host and between-host dynamics, Appl. Math. Model., 86 (2020), 404–423. https://doi.org/10.1016/j.apm.2020.05.022 doi: 10.1016/j.apm.2020.05.022
    [6] F. Brauer, C. Castillo-Chavez, Z. Feng, Mathematical Models in Epidemiology, Springer, New York, 2019.
    [7] M. Ghosh, P. Chandra, P. Sinha, J. B. Shukla, Modelling the carrier dependent infectious diseases with environmental effects, Appl. Math. Comput., 157 (2004), 385–402. https://doi.org/10.1016/S0096-3003(03)00564-2 doi: 10.1016/S0096-3003(03)00564-2
    [8] D. J. Rogers, S. E. Randolph, The global spread of malaria in a future warmer world, Science, 289 (2000), 1763–1766. https://doi.org/10.1126/science.289.5485.176 doi: 10.1126/science.289.5485.176
    [9] C. L. Mitchell, M. M. Janko, M. K. Mwandagalirwa, A. K. Tshefu, J. K. Edwards, B. W. Pence, et al., Impact of extractive industries on malaria prevalence in the Democratic Republic of the Congo: a population-based cross-sectional study, Sci. Rep., 12 (2022), 1–9.
    [10] F. J. Colón-González, M. O. Sewe, A. M. Tompkins, H. Sjödin, A. Casallas, J. Rocklöv, et al., Projecting the risk of mosquito-borne diseases in a warmer and more populated world: a multi-model, multi-scenario intercomparison modelling study, Lancet Planet. Health, 5 (2021), 404–414. https://doi.org/10.1016/S2542-5196(21)00132-7 doi: 10.1016/S2542-5196(21)00132-7
    [11] R. M. Anderson, R. M. May, Infectious Diseases of Humans: Dynamics and Control, Oxford University Press, London, U.K., 1991.
    [12] M. S. Arora, S. Singh, A. Omar, A. Malviya, J. B. Shukla, Effect of global warming temperature on the spread of bacteria dependent infectious diseases, Int. J. Clim. Change: Impacts Responses, 13 (2021), 1–19. https://doi.org/10.18848/1835-7156/CGP/v13i02/1-19 doi: 10.18848/1835-7156/CGP/v13i02/1-19
    [13] J. B. Shukla, V. Singh, A. K. Mishra, Modeling the spread of an infectious disease with bacteria and carriers in the environment, Nonlinear Anal.: Real World Appl., 12 (2011), 2541–2551. https://doi.org/10.1016/j.nonrwa.2011.03.003 doi: 10.1016/j.nonrwa.2011.03.003
    [14] S. Singh, Modeling the effect of global warming on the spread of carrier dependent infectious diseases, Model. Earth. Syst. Environ., 3 (2017), 3–39. https://doi.org/10.1007/s40808-017-0272-5 doi: 10.1007/s40808-017-0272-5
    [15] S. Singh, A. Omar, A. Modelling and analysis the effect of global warming on the spread of carrier dependent infectious diseases, Jnanabha, 50 (2020), 189–206.
    [16] A. ul Rehman, R. Singh, T. Abdeljawad, E. Okyere, L. Guran, Modeling, analysis and numerical solution to malaria fractional model with temporary immunity and relapse, Adv. Differ. Equations, 390 (2021). https://doi.org/10.1186/s13662-021-03532-4
    [17] S. Singh, J. B. Shukla, P. Chandra, Modelling and analysis of the spread of malaria: Environmental and ecological effects, J. Biol. Syst., 13 (2005), 1–11. https://doi.org/10.1142/S0218339005001367 doi: 10.1142/S0218339005001367
    [18] M. Ghosh, A. A. Lashari, X. Li, Biological control of malaria: A mathematical model, Appl. Math. Comput., 219 (2013), 7923–7939. https://doi.org/10.1016/j.amc.2013.02.053 doi: 10.1016/j.amc.2013.02.053
    [19] E. Ndamuzi, P. Gahungu, Mathematical modeling of Malaria transmission dynamics: Case of burundi, J. Appl. Math. Phys., 9 (2021), 2447–2460. https://doi.org/10.4236/jamp.2021.910156 doi: 10.4236/jamp.2021.910156
    [20] H. Wu, Z. Hu, Malaria transmission model with transmission-blocking drugs and a time delay, Math. Probl. Eng., (2021), 1–17. https://doi.org/10.1155/2021/1339086
    [21] B. Traoré, O. Koutou, B. Sangaré, A global mathematical model of malaria transmission dynamics with structured mosquito population and temperature variations, Nonlinear Anal.: Real World Appl., 53 (2020). https://doi.org/10.1016/j.nonrwa.2019.103081
    [22] A. Elaiw, A. A. Agha, Global analysis of a reaction-diffusion within-host Malaria infection model with adaptive immune response, Mathematics, 8 (2020), 563. https://doi.org/10.3390/math8040563 doi: 10.3390/math8040563
    [23] S. Noeiaghdam, A. Dreglea, H. Işık, M. Suleman, A comparative study between discrete stochastic arithmetic and floating-point arithmetic to validate the results of fractional order model of Malaria infection, Mathematics, 9 (2021), 1435. https://doi.org/10.3390/math9121435 doi: 10.3390/math9121435
    [24] S. Noeiaghdam, S. Micula, Dynamical strategy to control the accuracy of the nonlinear bio-mathematical model of Malaria infection, Mathematics, 9 (2021), 1031. https://doi.org/10.3390/math9091031 doi: 10.3390/math9091031
    [25] S. Singh, J. Singh, J. B. Shukla, Modelling and analysis of the effects of density dependent contact rates on the spread of carrier dependent infectious diseases with environmental discharges, Model. Earth Syst. Environ., 5 (2019), 21–32. https://doi.org/10.1007/s40808-018-0515-0 doi: 10.1007/s40808-018-0515-0
    [26] S. Tang, L. Chen, Density-dependent birth rate, birth pulses and their population dynamic consequences, J. Math. Biol., 44 (2002), 185–199. https://doi.org/10.1007/s002850100121 doi: 10.1007/s002850100121
    [27] H. W. Hethcote, Q. Linda, Disease transmission models with density-dependent demographics, Math. Biosci., 128 (1992), 93–130. https://doi.org/10.3817/0992093130 doi: 10.3817/0992093130
    [28] S. Singh, J. Singh, S. K. Sharma, J. B. Shukla, The effect of density dependent emigration on spread of infectious diseases: a modelling study, Rev. Latinoam. Hipertens., 14 (2019), 83–88.
    [29] J. Zhou, H. W. Hethcote, Population size dependent incidence in models for diseases without immunity, J. Math. Biol., 32 (1994), 809–834. https://doi.org/10.1007/BF00168799 doi: 10.1007/BF00168799
    [30] J. B. Shukla, S. Singh, J. Singh, S. K. Sharma, Modelling and analysis of bacteria dependent infectious diseases with variable contact rates, Comput. Materals Contin., 68 (2021), 1859–1875. https://doi.org/10.32604/cmc.2021.012095 doi: 10.32604/cmc.2021.012095
    [31] H. W. Hethcote, Qualitative analyses of communicable disease models, Math. Biosci., 28 (1976), 335–356. https://doi.org/10.1016/0025-5564(76)90132-2 doi: 10.1016/0025-5564(76)90132-2
    [32] D. Greenhalgh, Some threshold and stability results for epidemic models with a density dependent death rate, Theor. Popul. Biol., 42 (1990), 130–151. https://doi.org/10.1016/0040-5809(92)90009-I doi: 10.1016/0040-5809(92)90009-I
    [33] D. M. Bortz, P. W. Nelson, Sensitivity analysis of a nonlinear lumped parameter model of HIV infection dynamics, Bull. Math. Biol., 66 (2004), 1009–1026. https://doi.org/10.1016/j.bulm.2003.10.011 doi: 10.1016/j.bulm.2003.10.011
  • This article has been cited by:

    1. Muhammad Aslam Noor, Khalida Inayat Noor, Absolute Value Variational Inclusions, 2021, 2581-8147, 121, 10.34198/ejms.8122.121153
    2. Muhammad Aslam Noor, Khalida Inayat Noor, General Variational Inclusions and Nonexpansive Mappings, 2022, 2581-8147, 145, 10.34198/ejms.9222.145164
    3. Muhammad Aslam Noor, Khalida Inayat Noor, Some Novel Aspects of Quasi Variational Inequalities, 2022, 2581-8147, 1, 10.34198/ejms.10122.166
    4. Muhammad Aslam Noor, Khalida Inayat Noor, Savin Treanţă, Kamsing Nonlaopon, On three-step iterative schemes associated with general quasi-variational inclusions, 2022, 61, 11100168, 12051, 10.1016/j.aej.2022.05.031
    5. Muhammad Aslam Noor, Khalida Inayat Noor, Iterative Methods and Sensitivity Analysis for Exponential General Variational Inclusions, 2023, 2581-8147, 53, 10.34198/ejms.12123.53107
    6. Mudasir A. Malik, Mohd Iqbal Bhat, Bisma Zahoor, Solvability of a class of set-valued implicit quasi-variational inequalities: A Wiener–Hopf equation method, 2022, 9, 26667207, 100169, 10.1016/j.rico.2022.100169
    7. Saudia Jabeen, Jorge E. Macías-Díaz, Muhammad Aslam Noor, Muhammad Bilal Khan, Khalida Inayat Noor, Design and convergence analysis of some implicit inertial methods for quasi-variational inequalities via the Wiener–Hopf equations, 2022, 182, 01689274, 76, 10.1016/j.apnum.2022.08.001
    8. Muhammad Noor, Khalida Noor, New inertial approximation schemes for general quasi variational inclusions, 2022, 36, 0354-5180, 6071, 10.2298/FIL2218071N
    9. Mohd Iqbal Bhat, Mudasir Malik, Bisma Zahoor, Solvability of a Class of Setvalued Implicit Quasi-Variational Inequalities: A Wiener-Hopf Equation Method, 2022, 1556-5068, 10.2139/ssrn.4134292
    10. Muhammad Aslam Noor, Khalida Inayat Noor, Michael Th. Rassias, 2023, Chapter 13, 978-3-031-46486-7, 237, 10.1007/978-3-031-46487-4_13
    11. Muhammad Aslam Noor, Khalida Inayat Noor, Mixed Variational Inequalities and Nonconvex Analysis, 2024, 2581-8147, 873, 10.34198/ejms.14524.8731029
    12. Muhammad Aslam Noor, Khalida Inayat Noor, Some Computational Methods for Solving Extended General Bivariational Inclusions, 2023, 2581-8147, 133, 10.34198/ejms.13123.133163
    13. Muhammad Aslam Noor, Khalida Inayat Noor, New Iterative Methods and Sensitivity Analysis for Inverse Quasi Variational Inequalities, 2025, 2581-8147, 495, 10.34198/ejms.15425.495539
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1811) PDF downloads(86) Cited by(1)

Figures and Tables

Figures(12)  /  Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog