Loading [MathJax]/extensions/TeX/boldsymbol.js
Research article

Insecticidal Efficacy of Syzygium aromaticum, Tephrosia vogelii and Croton dichogamus Extracts against Plutella xylostella and Trichoplusia ni on Brassica oleracea crop in Northern Tanzania

  • The insecticidal efficacy of 10%, 5% and 1% w/v of Tephrosia vogelii, Croton dichogamus and Syzygium aromaticum aqueous plant extracts were assessed against larvae of Plutella xylostella and Trichoplusia ni on Brassica oleracea var. capitata crop field. Synthetic organophosphate pesticide (Chlorpyrifos) was used as a positive control and the negative controls were water and water plus soap. It was revealed that, aqueous plant extracts significantly (P ≤ 0.05) controlled the number of P. xylostella and T. ni larvae compared to negative controls. The 10% concentration of aqueous plant extracts showed significant higher efficacy in terms of reducing the insect population and their damage than 1% and 5 % concentration. The population of P. xylostella larvae per B. oleracea in five weeks of treatment applications at 10% w/v of T. vogelii, C. dichogamus and S. aromaticum aqueous plant extracts were 0.08, 0.15, 0.13, 0.05 and 0.08; 0.08, 0.20, 0.15, 0.13 and 0.18; 0.03, 0.05, 0.15, 0.18 and 0.13, respectively which was significantly (P ≤ 0.05) lower than in water (1.13, 1.68, 2.28, 2.20 and 3.28) and water plus soap (0.75, 1.60, 2.58, 1.83 and 3.30) negative controls respectively. The number of T. ni larvae per B. oleracea in five weeks of treatment applications at 10% w/v of T. vogelii, C. dichogamus and S. aromaticum of aqueous plant extracts were 0.00, 0.03, 0.05, 0.03 and 0.00; 0.03, 0.08, 0.05, 0.08 and 0.08; 0.05, 0.03, 0.00, 0.05 and 0.03, respectively which was significantly (P ≤ 0.05) lower than in water (0.50, 0.63, 0.60, 0.48 and 0.78) and water plus soap (0.30, 0.48, 68, 0.65 and 0.80) negative controls. The percentage damage of B. oleracea in five weeks of treatment applications at 10% w/v of T. vogelii, C. dichogamus and S. aromaticum aqueous plant extracts were 10.0, 6.3, 7.5, 7.5 and 5.6; 11.3, 11.3, 11.3, 11.3 and 12.5; 8.8, 8.1, 6.9, 8.1 and 9.4, respectively compared with water (33.8, 33.1, 38.8, 45.0 and 70.6) and water plus soap (30.0, 31.9, 41.3, 41.3 and 56.3). These pesticidal plants can be recommended for smallholder farmers to significantly control P. xylostella and T. ni larvae in B. oleracea crop.

    Citation: Nelson Mpumi, Revocatus L. Machunda, Kelvin M. Mtei, Patrick A. Ndakidemi. Insecticidal Efficacy of Syzygium aromaticum, Tephrosia vogelii and Croton dichogamus Extracts against Plutella xylostella and Trichoplusia ni on Brassica oleracea crop in Northern Tanzania[J]. AIMS Agriculture and Food, 2021, 6(1): 185-202. doi: 10.3934/agrfood.2021012

    Related Papers:

    [1] Li-Bin Liu, Yige Liao, Guangqing Long . Error estimate of BDF2 scheme on a Bakhvalov-type mesh for a singularly perturbed Volterra integro-differential equation. Networks and Heterogeneous Media, 2023, 18(2): 547-561. doi: 10.3934/nhm.2023023
    [2] Xiongfa Mai, Ciwen Zhu, Libin Liu . An adaptive grid method for a singularly perturbed convection-diffusion equation with a discontinuous convection coefficient. Networks and Heterogeneous Media, 2023, 18(4): 1528-1538. doi: 10.3934/nhm.2023067
    [3] Li-Bin Liu, Limin Ye, Xiaobing Bao, Yong Zhang . A second order numerical method for a Volterra integro-differential equation with a weakly singular kernel. Networks and Heterogeneous Media, 2024, 19(2): 740-752. doi: 10.3934/nhm.2024033
    [4] Dilip Sarkar, Shridhar Kumar, Pratibhamoy Das, Higinio Ramos . Higher-order convergence analysis for interior and boundary layers in a semi-linear reaction-diffusion system networked by a k-star graph with non-smooth source terms. Networks and Heterogeneous Media, 2024, 19(3): 1085-1115. doi: 10.3934/nhm.2024048
    [5] Chaoqun Huang, Nung Kwan Yip . Singular perturbation and bifurcation of diffuse transition layers in inhomogeneous media, part II. Networks and Heterogeneous Media, 2015, 10(4): 897-948. doi: 10.3934/nhm.2015.10.897
    [6] Yongqiang Zhao, Yanbin Tang . Approximation of solutions to integro-differential time fractional wave equations in Lpspace. Networks and Heterogeneous Media, 2023, 18(3): 1024-1058. doi: 10.3934/nhm.2023045
    [7] Chaoqun Huang, Nung Kwan Yip . Singular perturbation and bifurcation of diffuse transition layers in inhomogeneous media, part I. Networks and Heterogeneous Media, 2013, 8(4): 1009-1034. doi: 10.3934/nhm.2013.8.1009
    [8] Gianni Dal Maso, Francesco Solombrino . Quasistatic evolution for Cam-Clay plasticity: The spatially homogeneous case. Networks and Heterogeneous Media, 2010, 5(1): 97-132. doi: 10.3934/nhm.2010.5.97
    [9] Ciro D’Apice, Umberto De Maio, T. A. Mel'nyk . Asymptotic analysis of a perturbed parabolic problem in a thick junction of type 3:2:2. Networks and Heterogeneous Media, 2007, 2(2): 255-277. doi: 10.3934/nhm.2007.2.255
    [10] Rémi Goudey . A periodic homogenization problem with defects rare at infinity. Networks and Heterogeneous Media, 2022, 17(4): 547-592. doi: 10.3934/nhm.2022014
  • The insecticidal efficacy of 10%, 5% and 1% w/v of Tephrosia vogelii, Croton dichogamus and Syzygium aromaticum aqueous plant extracts were assessed against larvae of Plutella xylostella and Trichoplusia ni on Brassica oleracea var. capitata crop field. Synthetic organophosphate pesticide (Chlorpyrifos) was used as a positive control and the negative controls were water and water plus soap. It was revealed that, aqueous plant extracts significantly (P ≤ 0.05) controlled the number of P. xylostella and T. ni larvae compared to negative controls. The 10% concentration of aqueous plant extracts showed significant higher efficacy in terms of reducing the insect population and their damage than 1% and 5 % concentration. The population of P. xylostella larvae per B. oleracea in five weeks of treatment applications at 10% w/v of T. vogelii, C. dichogamus and S. aromaticum aqueous plant extracts were 0.08, 0.15, 0.13, 0.05 and 0.08; 0.08, 0.20, 0.15, 0.13 and 0.18; 0.03, 0.05, 0.15, 0.18 and 0.13, respectively which was significantly (P ≤ 0.05) lower than in water (1.13, 1.68, 2.28, 2.20 and 3.28) and water plus soap (0.75, 1.60, 2.58, 1.83 and 3.30) negative controls respectively. The number of T. ni larvae per B. oleracea in five weeks of treatment applications at 10% w/v of T. vogelii, C. dichogamus and S. aromaticum of aqueous plant extracts were 0.00, 0.03, 0.05, 0.03 and 0.00; 0.03, 0.08, 0.05, 0.08 and 0.08; 0.05, 0.03, 0.00, 0.05 and 0.03, respectively which was significantly (P ≤ 0.05) lower than in water (0.50, 0.63, 0.60, 0.48 and 0.78) and water plus soap (0.30, 0.48, 68, 0.65 and 0.80) negative controls. The percentage damage of B. oleracea in five weeks of treatment applications at 10% w/v of T. vogelii, C. dichogamus and S. aromaticum aqueous plant extracts were 10.0, 6.3, 7.5, 7.5 and 5.6; 11.3, 11.3, 11.3, 11.3 and 12.5; 8.8, 8.1, 6.9, 8.1 and 9.4, respectively compared with water (33.8, 33.1, 38.8, 45.0 and 70.6) and water plus soap (30.0, 31.9, 41.3, 41.3 and 56.3). These pesticidal plants can be recommended for smallholder farmers to significantly control P. xylostella and T. ni larvae in B. oleracea crop.


    Consider the following singularly perturbed Fredholm integro-differential equation (SPFIDE) in the interval ˉI=[0,T]:

    {Lu(t):=εu(t)+f(t,u(t))+λT0K(t,s,u(s))ds=0, tI=(0,T],u(0)=A, (1.1)

    where 0<ε1 is a perturbation parameter, A is a given constant and λ is a real parameter. We assume that f(t,u)C1(ˉI×R), K(t,s,u)C1(ˉI×ˉI×R) and there exist constants α, β such that 0<α|f/u|, |K/u|β. Under these assumptions, the problem (1.1) has a unique solution (see [1]).

    It is well known that the SPFIDEs arise widely in scientific fields such as mathematical biology [2], material mechanics [3], hydrodynamic [4] and so on.These problems depend on such a small positive parameter ε that the solution varies rapidly when ε0. Due to the presence of this perturbation parameter ε, classical numerical methods on a uniform mesh fail to give accurate results. Therefore, it is necessary to develop suitable numerical methods that are ε-uniformly convergent for solving these problems.

    Over the past few years, there has been a growing interest in the numerical methods for Volterra integro-differential equations (see, e.g., [5,6,7,8,9,10]) and Fredholm integro-differential equations (see, e.g., [11,12]). When the differential term of these integro-differential equations contains a small positive perturbation parameter ε, these problems are called singularly perturbed integro-differential equations. Recently, some robust numerical methods are proposed to solve singularly perturbed Volterra integro-differential equations[13,14,15,16]. Meanwhile, the authors in [17,18] developed fitted finite difference schemes on a uniform mesh for second-order SPFIDEs and gave some convergence results based on the prior information of the exact solution. Durmaz, et.al., [19] proposed a second-order uniformly convergent finite difference scheme on a Shishkin mesh for a singularly perturbed Fredholm integro-differential equation with the reduced second type Fredholm equation. In [20], the authors presented a fitted finite difference approach on a Shishkin mesh for a first-order singularly perturbed Fredholm integro-differential initial value problem with integral condition. Kumar et al. [21] proposed a non-standard finite difference scheme with Haar wavelet basis functions for a singularly perturbed partial integro-differential equation. Recently, Cakir et al., [22] solved first-order nonlinear SPFIDEs on a Sinshkin mesh with first-order convergence rate.

    From the literatures we have mentioned above, the existing numerical methods of SPFIDEs given in [19,20,21] are layer-adapted mesh approaches, which require a priori information about the location and width of the boundary layer. Thus, adaptive grid methods by equidistributing monitor functions are widely used to solve some singularly perturbed problems, see [25,26,27,28,29] for example. The advantages of these adaptive grid methods is to cluster automatically the grid points within the boundary layer. To the best of our knowledge, there is no report about this adaptive grid method for problem (1.1). Therefore, the aim of this paper is to solve problem (1.1) numerically by a finite difference scheme on an adaptive grid obtained by equidistributing a positive monitor function. The discrete scheme is constructed by using the backward Euler formula and right rectangle formula to approximate derivative term and nonlinear integral term, respectively. It is proved under some additional conditions that the proposed adaptive grid method is first-order uniformly convergent, with respect to ε.

    The rest of this paper is as follows. In Section 2, preliminary results of the exact solution are laid out. A discretization scheme is established in Section 3, in which a prior error analysis and a posterior error estimate are carried out successively. Numerical results obtained by the adaptive gird algorithm are given in Section 4 to support the theoretical analyses. The paper ends with a summary of the main conclusions in Section 5.

    Notation Throughout this paper, C which is not necessarily the same at each occurrence, indicates a positive constant independent of the mesh division parameter N or the singular perturbation parameter ε. To simplify the notation, we set vk=v(tk) for any function v(t). In our estimates, the maximum norm of a continuous function v(t) with the domain [0,T] is defined as v(t)=esssupt[0,T]|v(t)| as well as the maximum norm of a discrete vector \mathit{\boldsymbol{x}} = \left \{ x_i \right \}_{i = 0}^{N} with N+1 elements is defined as \left \| \mathit{\boldsymbol{x}} \right \|_{\infty} = \underset{i = 0, 1, 2, \cdots, N}{\max }\left | x_{i} \right | .

    In this section, we list the bounds for the exact solution u\left(t\right) and its first-order derivative.

    Lemma 2.1. [22, Lemma 1] Assume the constant \lambda satisfies

    \begin{equation} \left | \lambda \right | < \frac{\alpha }{\underset{0\leq t\leq T}{\max}\int_{0}^{T}\left | G\left ( t, s \right ) \right |ds}. \end{equation} (2.1)

    Then we have

    \begin{eqnarray} &&\left \| u \right \| _{\infty } \le C_{0}, \end{eqnarray} (2.2)
    \begin{eqnarray} && \left | u'\left(t\right) \right |\le C\left ( 1+\frac{1}{\varepsilon}e^{-\frac{\alpha t}{\varepsilon}} \right ) , \ 0\le t\le T, \end{eqnarray} (2.3)

    where

    \begin{eqnarray*} &&C_{0} = \frac{\left | A \right |+\frac{1}{\alpha}\left \|q\right \|_{\infty}}{1-\frac{1}{\alpha}\left | \lambda \right |\underset{0\leq t\leq T}{\max}\int_{0}^{T}\left | G\left ( t, s \right ) \right |ds}, \\ &&G\left(t, s\right) = \frac{\partial }{\partial u} K\left(t, s, \gamma u\right), \ 0 < \gamma < 1, \\ &&q\left(t\right) = -f\left(t, 0\right)-\gamma \int_{0}^{T} K\left(t, s, 0\right)ds. \end{eqnarray*}

    Corollary 2.1. For any two functions v\left(t\right) and w\left(t\right) satisfying

    \begin{equation} v\left(0\right) = w\left(0\right) = A, \ \end{equation} (2.4)

    and

    \begin{equation} Lv\left(t\right)-Lw\left(t\right) = \tilde{F}\left(t\right), \ t\in I, \ \end{equation} (2.5)

    where \tilde{F}(t) is a bounded piece-wise continuous function, we have

    \begin{equation} \left \| v\left(t\right)-w\left(t\right) \right \|_{\infty}\le C\left \| Lv\left(t\right)-Lw\left(t\right) \right \|_{\infty}.\ \end{equation} (2.6)

    Proof. The proof is similar to [15, Corollary 2.1].

    Let \bar{\Omega } _{N}: = \left\{0 = t_{0} < t_{1} < \dots < t_{N} = T \right\} be an arbitrary non-uniform mesh and h_{i} = t_{i}-t_{i-1} be the local mesh size for i = 1, 2, \cdots, N . For a given mesh function \left \{ v_{i} \right \} _{i = 0}^{N} , define the backward finite difference operator as follows:

    \begin{equation} D^-v_i = \frac{v_i-v_{i-1}}{h_i}, \ i = 1, 2, \cdots, N. \end{equation} (3.1)

    Then to construct the discretization scheme for problem (1.1), we integrate Eq (1.1) over \left(t_{i-1}, t_i \right) and use the right rectangle rule to approximate the integral part, which yield,

    \begin{equation} \left\{ \begin{aligned} &\varepsilon D^{-}u_{i}+f\left(t_i, u_i\right)+\lambda \sum\limits_{j = 1}^{N}h_jK\left(t_i, t_j, u_j\right)+R_i = 0, \ i = 1, 2, \cdots, N, \\ &u_0 = A, \end{aligned} \right. \end{equation} (3.2)

    where

    \begin{equation} R_i: = R_i^{\left(1\right)}+R_i^{\left(2\right)}+R_i^{\left(3\right)} \end{equation} (3.3)

    and

    \begin{eqnarray*} && R_i^{\left(1\right)} = -h_i^{-1}\int_{t_{i-1}}^{t_i}(t-t_{i-1})\frac{d }{d t} f\left(t, u\left(t\right)\right)dt, \\ &&R_i^{\left(2\right)} = -\lambda h_i^{-1}\int_{t_{i-1}}^{t_i}\left(t-t_{i-1}\right)\int_0^T\frac{\partial }{\partial t} K\left(t, s, u\left(s\right)\right)dsdt, \\ && R_i^{\left(3\right)} = -\lambda \sum\limits_{j = 1}^{N}\int_{t_{j-1}}^{t_j}\left(s-t_{j-1}\right)\frac{d }{d s} K\left(t_i, s, u\left(s\right)\right)ds.\ \end{eqnarray*}

    Neglecting the truncation error R_i in Eq (3.2), we obtain the discretization scheme of problem (1.1)

    \begin{equation} \left\{ \begin{aligned}&L^Nu_i^N: = \varepsilon D^-u_i^N+f\left(t_i, u_i^N\right)+\lambda\sum\limits_{j = 1}^{N}h_jK\left(t_i, t_j, u_j^N\right) = 0, \ i = 1, 2, \cdots, N, \\ &u_0 = A, \end{aligned}\right. \end{equation} (3.4)

    where u_i^N is the approximation of u(t) at point t = t_i .

    Let e_i^N: = u_i^N-u_i , i = 0, 1, \cdots, N , be the absolute error at t_i of the numerical solution. Then we can obtain the following error equations

    \begin{equation} \left\{ \begin{aligned} &L^Nu_i^N-L^Nu_i = R_i, \ i = 1, 2, \cdots, N, \\ &e_0^N = 0, \end{aligned}\right. \end{equation} (3.5)

    where R_i is the local truncation error defined in Eq (3.3) at t_i .

    Lemma 3.1. For i = 1, 2, \cdots, N , the truncation error R_i defined in Eq (3.3) satisfies

    \begin{equation} \left | R_i \right | \le C\max\limits_{1\le i\le N} \int_{t_{i-1}}^{t_i}\left(1+\left | {u}'\left(t\right) \right |\right)dt. \ \end{equation} (3.6)

    Proof. At first, based on the conditions f\left(t, u\right)\in C^1(\bar{I}\times R) and 0 < \alpha\le\left | \partial f/\partial u \right | , we obtain

    \begin{equation} \begin{aligned} \left | R_i^{\left(1\right)} \right |& \le h_{i}^{-1}\int_{t_{i-1}}^{t_i}h_i \left(\left | \frac{\partial f\left(t, u\right)}{\partial t} \right |+\left | \frac{\partial f\left(t, u\right)}{\partial u} \right | \left | u'\left(t\right) \right | \right)dt\\ &\le C\int_{ t_{i-1}}^{t_i}\left(1+\left | u'\left(t\right) \right |\right)dt.\ \end{aligned} \end{equation} (3.7)

    Then, since K\left(t, s, u\right)\in C^1\left(\bar{I}\times\bar{I}\times R\right) and \left | \partial K /\partial u \right | \le \beta , we have

    \begin{equation} \begin{aligned} \left | R_i^{\left(2\right)} \right | &\le\left | \lambda \right | \int_{t_{i-1}}^{t_i} \int_{0}^{T}\left | \frac{\partial }{\partial t}K\left(t, s, u\left(s\right)\right) \right |dsdt\\ & \le Ch_i.\ \end{aligned} \end{equation} (3.8)

    and

    \begin{equation} \begin{aligned} \left | R_i^{\left(3\right)} \right |&\le \left | \lambda \right | \sum\limits_{j = 1}^{N}T \int_{t_{j-1}}^{t_j}\left(\left | \frac{\partial K\left(t_i, s, u\right)}{\partial s} \right | +\left | \frac{\partial K\left(t_i, s, u\right)}{\partial u} \right |\left | u'\left(s\right) \right | \right)ds \\ &\le C\max\limits_{1\le j\le N}\int_{t_{j-1}}^{t_j}\left(1+\left | u'\left(t\right) \right |\right) dt.\ \end{aligned} \end{equation} (3.9)

    Finally, the desired result of this lemma can be followed by Eqs (3.7), (3.8) and (3.9).

    Lemma 3.2. Under the assumption

    \begin{equation} \left | \lambda \right | < \frac{1}{\alpha}\max\limits_{1\le i\le N}\sum\limits_{j = 1}^{N}h_j\left | G_{ij} \right |, \ \end{equation} (3.10)

    we have

    \begin{equation} \left \| \mathit{\boldsymbol{e}}^N \right \| _{\infty}\le \frac{1}{\alpha}\left(1-\frac{1}{\alpha}\left | \lambda \right |\max\limits_{1\le i\le N}\sum\limits_{j = 1}^{N}h_i\left | G_{ij} \right | \right)^{-1}\left \|\mathit{\boldsymbol{R}} \right \|_{\infty}, \ \end{equation} (3.11)

    where \mathit{\boldsymbol{e}}^N = \left \{ e_i^N \right \} _{i = 0}^{N} , \mathit{\boldsymbol{R}} = \left \{ R_i \right \} _{i = 0}^{N} and G_{ij} = \frac{\partial }{\partial u}K\left(t_i, s_j, u_j+\zeta e_j^N\right), \ 0 < \zeta < 1 .

    Proof. Applying the mean value theorem to Eq (3.5), we get

    \begin{equation} \varepsilon D^-e_i^N+a_ie_i^N+\lambda \sum\limits_{j = 1}^{N}h_jG_{ij}e_j^N = R_i, \ i = 1, 2, \cdots, N, \end{equation} (3.12)

    where

    \begin{eqnarray} &&a_i = \frac{\partial }{\partial u}f\left(t_i, u_i+\xi e_i^N\right), \ 0 < \xi < 1, \end{eqnarray} (3.13)
    \begin{eqnarray} &&G_{ij} = \frac{\partial }{\partial u}K\left(t_i, s_j, u_j+\zeta e_j^N\right), \ 0 < \zeta < 1 . \end{eqnarray} (3.14)

    According to maximum principle for the operator \varepsilon D^-e_i^N+a_ie_i^N , we have

    \begin{equation} \left \| \mathit{\boldsymbol{e}}^N \right \| _\infty \le \frac{1}{\alpha}\left \| \mathit{\boldsymbol{R}} \right \| _\infty +\frac{1}{\alpha}\left | \lambda \right | \left \| \mathit{\boldsymbol{e}}^N \right \| _\infty \max\limits_{1\le i\le N}\sum\limits_{j = 1}^{N}h_j\left | G_{ij} \right |, \end{equation} (3.15)

    which immediately leads to the desired result with the assumption (3.10).

    Based on the above Lemmas 3.1–3.2, we get the following convergence results.

    Theorem 3.1. Let u(t) be the solution of problem (1.1) and u_i^N be the solution of discrete scheme (3.4). Then

    \begin{equation} \max\limits_{1\le i\le N}\left | u_i^N-u_i \right | \le C \max\limits_{1\le i\le N}\int_{t_{i-1}}^{t_i}\left(1+\left | {u}'\left(t\right) \right |\right)dt. \end{equation} (3.16)

    Corollary 3.1. Under the conditions of Theorem 3.1, there exists an adaptive grid \left \{ t_i \right \} _{i = 0}^{N} such that

    \begin{equation} \max\limits_{1\le i\le N}\left | u_i^N-u_i \right |\le CN^{-1}. \ \end{equation} (3.17)

    Proof. Based on the mesh equidistribution principle presented in [29], the mesh \left \{ t_i \right \}_{i = 0}^{N} given by our adaptive grid algorithm satisfies

    \begin{equation} \int_{t_{i-1}}^{t_i}M\left(t\right)dt = \frac{1}{N}\int_{0}^{T}M\left(t\right)dt, \ i = 1, 2, \cdots, N, \end{equation} (3.18)

    where M(t) is called the monitor function, which can be chosen as

    \begin{equation} M\left(t\right) = 1+\left | {u}'\left(t\right) \right |.\ \end{equation} (3.19)

    Therefore, it follows from Lemma 2.1 that

    \begin{equation} \begin{aligned} \max\limits_{1\le i\le N}\left | u_i^N-u_i \right | &\le C \max\limits_{1\le i\le N}\int_{t_{i-1}}^{t_i}\left(1+\left | {u}'\left(t\right) \right |\right)dt\\ & = \frac{C}{N}\int_{0}^{T}\left(1+\left | {u}'\left(t\right) \right | \right)dt\\ &\le \frac{C}{N}\int_{0}^{T}\left(1+\frac{1}{\varepsilon }\exp\left(-\frac{\alpha t}{\varepsilon } \right)\right)dt\\ &\le \frac{C}{N}\left(T+\frac{1}{\alpha}\left(1-\exp\left(-\frac{\alpha T}{\varepsilon } \right)\right)\right)\\ &\le\frac{C}{N}. \end{aligned} \end{equation} (3.20)

    In this section, we shall derive an a posteriori error estimation for the numerical solution \left \{ u_i^N \right \} _{i = 0}^{N} . Recall that \tilde{u}^N\left(t\right) is a piece-wise linear interpolation function through knots \left(t_i, u_i^N\right), \ i = 0, 1, \cdots, N . Then, for any t\in J_i: = \left [t_{i-1}, t_i \right] , we obtain

    \begin{equation} \tilde{u}^N\left(t\right) = u_i^N+D^-u_i^N\left(t-t_i\right), \ i = 1, 2, \cdots, N.\ \end{equation} (3.21)

    Theorem 3.2. Let u(t) be the exact solution of problem (1.1), \left \{ u_i^N \right \} _{i = 0}^{N} be the discrete solution of problem (3.4) and \tilde{u}^N(t) be its piece-wise linear interpolation function defined in Eq (3.21). Then we have

    \begin{equation} \left \| \tilde{u}^N\left(t\right)-u\left(t\right) \right \|_{\infty}\le C \max\limits_{1\le i\le N}\left(h_i+h_i\left | D^-u_i^N \right | \right).\ \end{equation} (3.22)

    Proof. For any t\in\left (t_{i-1}, t_i \right] , it follows from Eq (1.1) and Eq (3.4) that

    \begin{equation} \begin{aligned} L\tilde{u}^N\left(t\right)-Lu\left(t\right) & = \varepsilon D^-u_i^N+f\left(t, \tilde{u}^N\left(t\right)\right)+\lambda \int_0^TK\left(t, s, \tilde{u}^N\left(s\right)\right)ds\\ & = -f\left(t_i, u_i^N\right)-\lambda \sum\limits_{j = i}^{N}h_jK\left(t_i, t_j, u_j^N\right) \\ &\quad +f\left(t, \tilde{u}^N\left(t\right)\right)+\lambda \int_0^TK\left(t, s, \tilde{u}^N\left(s\right)\right)ds\\ & = P\left(t\right)+Q\left(t\right), \ \end{aligned} \end{equation} (3.23)

    where

    \begin{eqnarray} &&P\left(t\right) = f\left(t, \tilde{u}^N\left(t\right)\right)-f\left(t_i, u_i^N\right), \end{eqnarray} (3.24)
    \begin{eqnarray} && Q\left(t\right) = \lambda \int_{0}^{T}K(t, s, \tilde{u}^N\left(s)\right)ds-\lambda \sum\limits_{j = 1}^{N}h_jK\left(t_i, t_j, u_j^N\right) . \end{eqnarray} (3.25)

    With the assumptions of functions f\left(t, u\right) , K\left(t, s, u\right) and the definition of \tilde{u}^N\left(t\right) , we have

    \begin{equation} \begin{aligned} \left | P\left(t\right) \right | & = \left | f\left(t_i, u_i^N\right)+\int_{t_i}^{t}\frac{df\left(\tau, \tilde{u}^N\left(\tau\right)\right)}{d\tau}d\tau-f\left(t_i, u_i^N\right) \right | \\ &\le \int_{t}^{t_i}\left(\left | \frac{\partial f\left(\tau, \tilde{u}^N\right)}{\partial \tau} \right |+ \left | \frac{\partial f\left(\tau, \tilde{u}^N\right)}{\partial u} \right |\left | D^-u_i^N \right |\right) d\tau \\ &\le C h_i\left(1+\left | D^-u_i^N \right |\right) \end{aligned} \end{equation} (3.26)

    and

    \begin{equation} \begin{aligned} \left | Q\left(t\right) \right |& = \left | \lambda \sum\limits_{j = 1}^{N}\int_{t_{j-1}}^{t_j}\left(K\left(t, s, \tilde{u}^N\left(s\right)\right)-K\left(t_i, t_j, u_j^N\right)\right)ds \right | \\ &\le\left |\lambda\right | \sum\limits_{j = 1}^{N} \int_{t_{j-1}}^{t_j}\left(\left |\frac{\partial }{\partial t} K\left(\xi_1t+\left(1-\xi_1\right)t_i, s, \tilde{u}^N\left(s\right)\right)\left(t-t_{i}\right)\right |\right.\\ &\quad+\left |\frac{\partial }{\partial s} K\left(t_i, \xi_2s+\left(1-\xi_2\right)t_{j}, \tilde{u}^N\left(s\right)\right)\left(s-t_{j}\right)\right |\\ &\left. \quad+\left | \frac{\partial }{\partial u}K\left(t_i, t_{j}, \xi_3\tilde{u}^N\left(s\right)+\left(1-\xi_3\right)u_{j}^N\right)\left(\tilde{u}^N\left(s\right)-u_{j }^N\right)\right |\right) ds\\ &\le C\left(h_i+\max\limits_{1\le j\le N} h_j \left(1+\left | D^-u_j^N \right |\right)\right), \ \end{aligned} \end{equation} (3.27)

    where 0 < \xi_1 < 1 , 0 < \xi_2 < 1 , and 0 < \xi_3 < 1 . The result can be derived from Eqs (3.23), (3.26), (3.27) and Corollary 2.1.

    From Corollary 3.1, it is easy to conclude that there exists a mesh \left \{ t_i \right \}_{i = 0}^{N} and a monitor function M\left(t\right) given in Eq (3.19) such that the inequality (3.17) holds true. However, {u}'\left(t\right) is not available. Therefore, based on the a posterior error estimation (3.22), we choose the discrete analogue of M\left(t\right) as

    \begin{equation} {\tilde{M}_i = 1+\left | D^-u_i^N \right |, i = 1, 2, \cdots, N.\ } \end{equation} (3.28)

    Therefore, the idea is to adaptively design a mesh in which the values of monitor function tilde{M}_i are the same on each mesh interval. This is equivalent to find \left \{ \left(t_i, u_i^N\right) \right \}_{i = 0}^{N} , such that

    \begin{equation} h_i\tilde{M}_i = \frac{1}{N}\sum\limits_{j = 1}^{N}h_j\tilde{M}_j, \ i = 1, 2, \cdots, N.\ \end{equation} (3.29)

    Furthermore, to obtain this equidistributed mesh \left\{t_i\right\}_{i = 0}^N and the corresponding numerical solution u_i^N , we give the following iteration algorithm:

    Algorithm 1 Steps of adaptive grid algorithm
    1: Step 1: For a given N , let \{t_i^{(0)}\} _{i = 0}^{N} be an initial mesh with mesh step \frac{1}{N} . Choose a constant \mu^{*} > 1 that controls when the algorithm terminates.
    2: Step 2: For a given mesh \{t_i^{(k)}\}_{i = 0}^{N} and numerical solution \{u_i^{N, (k)}\}_{i = 0}^{N} , compute \tilde{M}_i^{(k)}, i = 1, 2, \cdots, N in Eq (3.28) and set \tilde{M}_0^{(k)} = 0 .
    3: Step 3: Set h_i^{(k)} = t_i^{(k)}-t_{i-1}^{(k)} for each i and set L_0^{(k)} = 0 and L_i^{(k)} = \sum_{j = 1}^{i}h_j^{(k)}\tilde{M}_j^{(k)} for i = 1, 2, \cdots, N . Define
                                            \begin{array}{*{20}{c}} {\mu^{(k)}: = \frac{N}{L_N^{(k)}}\max\limits_{i = 0, 1, \cdots, N}h_i^{(k)}\tilde{M}_i^{(k)}.} & {} & {} & {} & {(3.30)} \end{array}
    4: Step 4: Set Y_i^{(k)} = iL_N^{(k)}/N for i = 0, 1, \cdots, N . Interpolate (see [30, Remark 5.1]) to the points (L_i^{(k)}, t_i^{(k)}) . Generate the new mesh \{t_i^{(k+1)}\}_{i = 0}^{N} by evaluating this interpolant at the Y_i^{(k)} for i = 0, 1, \cdots, N .
    5: Step 5: If \mu^{(k)}\leq\mu^* , then take \{t_i^{(k+1)}\}_{i = 0}^{N} as the final mesh and compute \{u_i^{N, (k+1)}\}_{i = 0}^{N} . Otherwise return to Step 2.

    In Section 4.1, We first present the iterative scheme. Then numerical experiments are given in Section 4.2 to validate the theoretical result of this paper. All experiments were performed on a Windows 10 (64 bit) PC-Intel(R) Core(TM) i5-4200H CPU 2.80 GHz, 8 GB of RAM using MATLAB R2021a.

    In order to avoid solving the nonlinear equations (3.4), we apply the quasilinearization technique that performs a first-order Taylor expansion on the last iteration values and obtain

    \begin{equation} \left\{ \begin{aligned}&u_{i}^{N, \left(k\right)} = \frac{\varepsilon /h_i u_{i-1}^{N, \left(k\right)}+B_iu_i^{N, \left(k-1\right)}+C_i }{\varepsilon/h_i +B_i }, i = 1, 2, \cdots, N, \\ &u_0^{N, \left(k\right)} = A, \end{aligned} \right. \end{equation} (4.1)

    where

    \begin{eqnarray} &&B_i = \frac{\partial }{\partial u} f\left(t_i, u_i^{N, \left(k-1\right)}\right), \end{eqnarray} (4.2)
    \begin{eqnarray} &&C_i = -f\left(t_i, u_{i}^{N, \left(k-1\right)}\right)-\lambda \sum\limits_{j = 1}^{N}h_jK\left(t_i, t_j, u_j^{N, \left(k-1\right)}\right). \end{eqnarray} (4.3)

    For all the numerical experiments below, we choose \mu^* = 1.1 , which is defined in Step 4 in Algorithm 1.

    Example 4.1. We consider a SPFIDE in the form [22]

    \begin{equation} \left\{ \begin{aligned} &\varepsilon {u}'\left(t\right)+2u\left(t\right)+\tanh \left(u\left(t\right)\right)-e^t+\frac{1}{4}\int_{0}^{1}t^2\sin \left(u\left(s\right)\right)ds = 0, \ t\in\left(0, 1\right], \\ &u(0) = 1. \end{aligned} \right. \end{equation} (4.4)

    Since the analytic solution of this problem is not available, we use the following formulas to calculate the errors and the corresponding convergence rates:

    \begin{equation} e_{\varepsilon }^N = \left \| \hat{\mathit{\boldsymbol{u}}}^{2N}-\mathit{\boldsymbol{u}}^N \right \| _{\infty}, \ p_{\varepsilon }^N = \log_{2}{\left(\frac{e_{\varepsilon }^{N} }{e_{\varepsilon }^{2N} } \right)} , \ \end{equation} (4.5)

    where \hat{\mathit{\boldsymbol{u}}}^{2N} is the numerical solution obtained on the fine mesh \bar{\Omega } _{2N} = \bar{\Omega } _{N}\bigcup \left\{\frac{t_{i}+t_{i+1}}{2}\right\}_{i = 0}^{N-1} . The maximum errors and \varepsilon -uniform rates of the convergence are respectively defined as

    \begin{equation} e^N = \max\limits_\varepsilon e^N_\varepsilon, \ p^N = \log_{2}{\left(\frac{e^{N}}{e^{2N}}\right)}. \end{equation} (4.6)

    In the numerical experiments, we apply the presented adaptive grid algorithm to solve this problem. The resulting errors e_{\varepsilon }^N and the orders of the convergence p_{\varepsilon }^N , for particular values of \varepsilon and N are listed in Table 1. In addition, to compare the performance of the presented adaptive mesh with the Shishkin mesh [22] and the Bakhvalov mesh [23], some numerical results are given in Table 2.

    Table 1.  The errors and corresponding convergence rates for Example 4.1.
    \varepsilon N=64 N=128 N=256 N=512 N=1024
    2^{-2} 0.004663 0.002384 0.001222 0.000619 0.000312
    0.9680 0.9637 0.9811 0.9903 -
    2^{-4} 0.00595 0.003152 0.00163 0.000801 0.000412
    0.9167 0.9515 1.0244 0.9610 -
    2^{-6} 0.006539 0.003482 0.001818 0.000939 0.000477
    0.9089 0.9373 0.9532 0.9769 -
    2^{-8} 0.006711 0.003633 0.001893 0.000979 0.000501
    0.8853 0.9404 0.9514 0.9667 -
    2^{-10} 0.006799 0.003635 0.001913 0.000996 0.000509
    0.9035 0.9258 0.9424 0.9667 -
    2^{-12} 0.006881 0.003649 0.001926 0.000996 0.000512
    0.9151 0.9221 0.9505 0.9606 -
    2^{-14} 0.006854 0.003683 0.001924 0.001 0.000513
    0.8963 0.9365 0.9447 0.9633 -
    2^{-16} 0.00672 0.003752 0.001934 0.000999 0.000515
    0.8411 0.9562 0.9526 0.9560 -
    e^N 0.006881 0.003752 0.001934 0.001 0.000515
    p^N 0.8751 0.9562 0.9517 0.9569 -

     | Show Table
    DownLoad: CSV
    Table 2.  Comparisons of errors and corresponding convergence rates for Example 4.1.
    N \varepsilon=2^{-6} \varepsilon=2^{-8}
    Adaptive Bakhvalov Shishkin Adaptive Bakhvalov Shishkin
    3\times2^5 4.53E-03 2.81E-03 7.18E-03 4.66E-03 2.83E-03 7.16E-03
    0.93 0.98 0.75 0.91 0.98 0.75
    3\times2^6 2.39E-03 1.42E-03 4.26E-03 2.48E-03 1.44E-03 4.25E-03
    0.94 0.99 0.79 0.95 0.99 0.79
    3\times2^7 1.24E-03 7.16E-04 2.46E-03 1.29E-03 7.23E-04 2.45E-03
    0.97 0.99 0.83 0.96 0.99 0.83
    3\times2^8 6.33E-04 3.59E-04 1.39E-03 6.62E-04 3.63E-04 1.38E-03
    0.98 1.00 0.85 0.97 1.00 0.85
    3\times2^9 3.20E-04 1.80E-04 7.71E-04 3.37E-04 1.82E-04 7.69E-04
    1.05 1.00 0.86 0.98 1.00 0.86
    3\times2^{10} 1.54E-04 9.01E-05 4.23E-04 1.71E-04 9.09E-05 4.22E-04

     | Show Table
    DownLoad: CSV

    According to the results in Table 1, for fixed N , the error increases with a diminishing speed and the convergence rate goes away from 1 as \varepsilon decreases, while for fixed \varepsilon , the error deceases to half as N doubles and the convergence rate is getting closer to 1. According to the results in Table 2, our adaptive mesh demonstrates both less error and more accurate first-order convergence rate than Shishkin mesh. In general, the performance of our adaptive grid algorithm is better when \mu^* which is used to distribute the grid uniformly is closer to 1. However, the enhance becomes very limited after certain threshold and it can not work out successfully when \mu^* is so close to 1 due to the enormous amount of calculations. Here the limitation of the adaptive mesh lies.

    The behaviors of the numerical solution are presented in Figures 1, 2. Obviously, it can be seen that in these two figures the solution of the test problem decreases successively near to 0 at first and increases progressively close to 1 then and has a boundary layer at t = 0 . Figure 3 shows the \varepsilon -uniform convergence of the method with different \varepsilon . To be more physical, the first-order uniform convergence stands for our method in spite of violent changes of the numerical solution in the bound layer at t = 0 . From Figure 3, no matter how close \varepsilon tends to zero, the maximum point-wise errors are bounded by O\left(N^{-1}\right) which validates our theoretical analyses. For \varepsilon = 2^{-8} , Figure 4 displays how a mesh with N = 64 evolves through successive iterations of the algorithm by using monitor function (3.28).

    Figure 1.  Numerical results of Example 4.1 for N = 64 and \varepsilon = 2^{-4} .
    Figure 2.  Numerical results of Example 4.1 for N = 256 and various \varepsilon .
    Figure 3.  Maximum point-wise errors of log-log plot for Example 4.1.
    Figure 4.  Evolution of the mesh for Example 4.1.

    Example 4.2. We consider a first-order nonlinear singularly perturbed mixed type integro-differential equation in [24] in the form:

    \begin{equation} \left\{ \begin{aligned} &\varepsilon u'\left(t\right)+f\left(t, u\left(t\right)\right)+\frac{1}{2}\int_{0}^{t}u^2(s)ds+\frac{1}{2}\int_{0}^{1}u^3\left(s\right)ds = 0, \ t\in\left(0, 1\right], \\ &u(0) = 1. \end{aligned} \right. \end{equation} (4.7)

    where

    \begin{equation} f\left(t, u\left(t\right)\right) = \frac{\varepsilon}{4}u^2\left(t\right)+u\left(t\right)+\frac{\varepsilon }{6}e^{-\frac{3}{\varepsilon }}-\frac{5}{12}\varepsilon. \end{equation} (4.8)

    It is testified that it satisfies the assumption (3.10). The analytic solution of this problem is u\left(t\right) = e^{-t/\varepsilon} . We use the following formulas to calculate the errors and the corresponding convergence rates:

    \begin{equation} e_{\varepsilon }^N = \max\limits_{0\le i\le N} \left | u_i^{N}-u_i \right | , \ p_{\varepsilon }^N = \log_{2}{\left(\frac{e_{\varepsilon }^{N} }{e_{\varepsilon }^{2N} } \right)} .\ \end{equation} (4.9)

    The maximum errors and \varepsilon -uniform convergence rates are defined in Eq (4.6).

    The resulting errors e_{\varepsilon }^N and the orders of the convergence p_{\varepsilon }^N , for particular values of \varepsilon and N are listed in Table 3. In addition, to compare the performance of the presented adaptive mesh with the Shishkin mesh [22] and the Bakhvalov mesh [23], some numerical results are given in Table 4.

    Table 3.  The errors and corresponding convergence rates for Example 4.2.
    \varepsilon N=64 N=128 N=256 N=512 N=1024
    2^{-2} 0.008698 0.004386 0.002202 0.001103 0.000552
    0.9880 0.9939 0.9969 0.9985 -
    2^{-4} 0.011073 0.00562 0.002799 0.001413 0.00071
    0.9783 1.0057 0.9863 0.9929 -
    2^{-6} 0.012473 0.00639 0.003247 0.001638 0.0008
    0.9648 0.9767 0.9876 1.0335 -
    2^{-8} 0.013236 0.006774 0.003463 0.00176 0.000887
    0.9664 0.9681 0.9763 0.9881 -
    2^{-10} 0.013376 0.006909 0.003559 0.00181 0.000917
    0.9530 0.9570 0.9755 0.9814 -
    2^{-12} 0.013461 0.006954 0.003575 0.001835 0.000929
    0.9528 0.9601 0.9619 0.9822 -
    2^{-14} 0.013563 0.007002 0.003592 0.001836 0.000937
    0.9540 0.9631 0.9680 0.9704 -
    2^{-16} 0.015521 0.007146 0.003596 0.001842 0.000935
    1.1190 0.9908 0.9647 0.9784 -
    e^N 0.015521 0.007146 0.003596 0.001842 0.000937
    p^N 1.1190 0.9908 0.9647 0.9755 -

     | Show Table
    DownLoad: CSV
    Table 4.  Comparisons of errors and corresponding convergence rates for Example 4.2.
    N \varepsilon=2^{-6} \varepsilon=2^{-8}
    Adaptive Bakhvalov Shishkin Adaptive Bakhvalov Shishkin
    3\times2^5 8.43E-03 2.15E-02 1.91E-02 8.95E-03 2.69E-02 1.67E-02
    0.97 0.90 1.28 0.97 0.65 0.11
    3\times2^6 4.31E-03 1.15E-02 7.85E-03 4.58E-03 1.72E-02 1.54E-02
    0.98 0.94 1.29 0.97 0.67 0.80
    3\times2^7 2.18E-03 6.00E-03 3.21E-03 2.33E-03 1.08E-02 8.88E-03
    1.04 0.97 1.00 0.98 0.90 1.28
    3\times2^8 1.06E-03 3.06E-03 1.61E-03 1.18E-03 5.78E-03 3.65E-03
    0.97 0.98 0.85 0.99 0.93 1.32
    3\times2^9 5.39E-04 1.55E-03 8.88E-04 5.93E-04 3.02E-03 1.46E-03
    0.98 0.99 0.87 1.00 0.97 1.41
    3\times2^{10} 2.72E-04 7.78E-04 4.87E-04 2.98E-04 1.54E-03 5.51E-04

     | Show Table
    DownLoad: CSV

    The behaviors of the numerical solution are presented in Figures 5, 6. Obviously, it can be seen that in these two figures the solution of the test problem gradually decreases to zero with a decreasing velocity rate as well as has a boundary layer at t = 0 . Figure 7 shows the \varepsilon -uniform convergence of the method with different \varepsilon . From Figure 7, no matter how close \varepsilon tends to zero, the maximum point-wise errors are bounded by O\left(N^{-1}\right) which validates our theoretical analyses. For \varepsilon = 2^{-8} , Figure 8 displays how a mesh with N = 64 evolves through successive iterations of the algorithm by using monitor function (3.28).

    Figure 5.  Numerical results of Example 4.2 for N = 64 and \varepsilon = 2^{-4} .
    Figure 6.  Numerical results of Example 4.2 for N = 256 and various \varepsilon .
    Figure 7.  Maximum point-wise errors of log-log plot for Example 4.2.
    Figure 8.  Evolution of the mesh for Example 4.2.

    The unique but significant difference between these two examples can be found in Table 4. At the situation \varepsilon = 2^{-8} , adaptive mesh demonstrates better than Bakhvalov mesh. This is because the construction of the Bakhvalov mesh is based on prior information of the exact solution. Therefore the effect of this method is tightly relevant to the form of equation. The existing Bakhvalov mesh may display wonderful at suitable equations but poor when it is adverse. However, the adaptive mesh based on posterior error estimate requires none of prior information of the exact solution and displays steadily at various equations.

    This paper considers a nonlinear singularly perturbed Fredholm integro-differential equation for a first-order initial value problem. A first-order \varepsilon -uniformly convergent numerical method for solving this problem is presented, which comprises an adaptive grid based on the posterior error estimate and the mesh equidistribution principle. The difference scheme including the weights and the remainders is established using the backward Euler formula for the derivative term, together with the composite right rectangle quadrature rule for the integral term. A theoretical analysis based on prior error bound is conducted to prove the first-order convergence of the proposed method. Two examples show that our adaptive mesh demonstrates better than Shishkin mesh and performs as well as Bakhvalov mesh. In the future, we will try to extend the proposed adaptive grid method for solving other related integro-differential equations which can be found in [5,6,31].

    The authors declare there is no conflict of interest.



    [1] Baidoo P, Mochiah M (2016) Comparing the effectiveness of garlic (Allium sativum L.) and hot pepper (Capsicum frutescens L.) in the management of the major pests of cabbage Brassica oleracea (L.). Sustainable Agric Res 5: 83-91.
    [2] Mpumi N, Machunda RS, Mtei KM, et al. (2020) Selected insect pests of economic importance to Brassica oleracea, their control strategies and the potential threat to environmental pollution in Africa. Sustainability 12: 1-22.
    [3] Mudzingwa S, Muzemu S, Chitamba J (2013) Pesticidal efficacy of crude aqueous extracts of Tephrosia vogelii L., Allium sativum L. and Solanum incanum L. in controlling aphids (Brevicoryne brassicae L.) in rape (Brassica napus L.). J Agric Res 2: 157-163.
    [4] Norman C, Shealy M (2007) Illustrated Encyclopedia of Healing Remedies. Toledo, OH, USA: Elements Book Inc.
    [5] Gibson A (2012) Colewart and the cole crops. Los Angeles: University of California Los Angeles.
    [6] Baidoo P, Adam J (2012) The Effects of Extracts of Lantana camara (L.) and Azadirachta indica (A. Juss) on the Population Dynamics of Plutella xylostella, Brevicoryne brassicae and Hellula undalis on Cabbage. Sustainable Agric Res 1: 229-234.
    [7] Furlong MJ, Wright DJ, Dosdall LM (2013) Diamondback moth ecology and management: problems, progress, and prospects. Annu Rev Entomol 58: 517-541.
    [8] Labou B, Brévault T, Sylla S, et al. (2017) Spatial and temporal incidence of insect pests in farmers' cabbage fields in Senegal. Int J Trop Insect Sci 37: 225-233.
    [9] Begna F, Damtew T (2015) Evaluation of four botanical insecticides against Diamondback Moth, Plutella Xylostella L.(Lepidoptera: Plutellidae) on head cabbage in the central rift valley of Ethiopia. Sky J Agric Res 4: 97-105.
    [10] Timbilla J, Nyarko K (2004) A survey of cabbage production and constraints in Ghana. Ghana J Agric Sci 37: 93-101.
    [11] Fening K, Amoabeng B, Adama I, et al. (2013) Sustainable management of two key pests of cabbage, Brassica oleracea var. capitata L.(Brassicaceae), using homemade extracts from garlic and hot pepper. Org Agric 3: 163-173.
    [12] Capinera J (2001) Handbook of vegetable pests. Eds, editor: Cambridge, MA, USA: Academic Press.
    [13] Mpumi N, Mtei K, Machunda R, et al. (2016) The toxicity, persistence and mode of actions of selected botanical pesticides in Africa against insect pests in common beans, P. vulgaris: a review. Am J Plant Sci 7: 138-151.
    [14] Labou B, Brévault T, Bordat D, et al. (2016) Determinants of parasitoid assemblages of the diamondback moth, Plutella xylostella, in cabbage farmer fields in Senegal. Crop Prot 89: 6-11.
    [15] Mkenda P, Mwanauta R, Stevenson PC, et al. (2015) Extracts from field margin weeds provide economically viable and environmentally benign pest control compared to synthetic pesticides. PLoS ONE 10: 1-14.
    [16] Mkindi A, Mpumi N, Tembo Y, et al. (2017) Invasive weeds with pesticidal properties as potential new crops. Ind Crops Prod 110: 113-122.
    [17] Ondieki JJ (1996) The current state of pesticide management in Sub-Saharan Africa. Sci Total Environ 188: 30-34.
    [18] Weinberger K, Srinivasan R (2009) Farmers' management of cabbage and cauliflower pests in India and their approaches to crop protection. J Asia-Pac Entomol 12: 253-259.
    [19] Amoabeng BW, Gurr GM, Gitau CW, et al. (2014) Cost: Benefit analysis of botanical insecticide use in cabbage: Implications for smallholder farmers in developing countries. Crop Prot 57: 71-76.
    [20] de Cássia Seffrin R, Shikano I, Akhtar Y, et al. (2010) Effects of crude seed extracts of Annona atemoya and Annona squamosa L. against the cabbage looper, Trichoplusia ni in the laboratory and greenhouse. Crop Prot 29: 20-24.
    [21] Mochiah M, Banful B, Fening K, et al. (2011) Botanicals for the management of insect pests in organic vegetable production. J Entomol Nematol 3: 85-97.
    [22] Isman MB (2006) Botanical insecticides, deterrents, and repellents in modern agriculture and an increasingly regulated world. Annu Rev Entomol 51: 45-66.
    [23] Orwa C, Mutua A, Kindt R, et al. (2009) Agroforestree database: a tree species reference and selection guide version 4.0. World Agroforestry Centre ICRAF, Nairobi, KE.
    [24] Munthali MG, Gachene CK, Sileshi GW, et al. (2014) Amendment of Tephrosia Improved Fallows with Inorganic Fertilizers Improves Soil Chemical Properties, N Uptake, and Maize Yield in Malawi. Int J Agron 2014: 1-9.
    [25] Kamatou GP, Vermaak I, Viljoen AM (2012) Eugenol from the remote Maluku Islands to the international market place: A review of a remarkable and versatile molecule. Molecules 17: 6953-6981.
    [26] Tian BL, Liu QZ, Liu ZL, et al. (2015) Insecticidal potential of clove essential oil and its constituents on Cacopsylla chinensis (Hemiptera: Psyllidae) in laboratory and field. J Econ Entomol 108: 957-961.
    [27] Araujo AFdO, Ribeiro-Paes JT, Deus JTd, et al. (2016) Larvicidal activity of Syzygium aromaticum (L.) Merr and Citrus sinensis (L.) Osbeck essential oils and their antagonistic effects with temephos in resistant populations of Aedes aegypti. Memórias do Instituto Oswaldo Cruz 111: 443-449.
    [28] Silva CG, Oliveira J, da Camara CA (2018) Insecticidal activity of the ethanolic extract from Croton species against Plutella xylostella L.(Lepidoptera: Plutellidae). Revista Facultad Nacional de Agronomía Medellín 71: 8543-8551.
    [29] Aldhaher A, Langat M, Ndunda B, et al. (2017) Diterpenoids from the roots of Croton dichogamus Pax. Phytochemistry 144: 1-8.
    [30] Belmain SR, Amoah BA, Nyirenda SP, et al. (2012) Highly variable insect control efficacy of Tephrosia vogelii chemotypes. J Agric Food Chem 60: 10055-10063.
    [31] Patra S, Dhote V, Alam S, et al. (2013) Population dynamics of major insect pests and their natural enemies on cabbage under new alluvial zone of West Bengal. J Plant Prot Sci 5: 42-49.
    [32] Tanyi CB, Ngosong C, Ntonifor NN (2018) Effects of climate variability on insect pests of cabbage: adapting alternative planting dates and cropping pattern as control measures. Chem Biol Technol Agric 5: 1-11.
    [33] Ayalew G, Ogol C (2006) Occurrence of the diamondback moth (Plutella xylostella L.) and its parasitoids in Ethiopia: influence of geographical region and agronomic traits. J Appl Entomol 130: 343-348.
    [34] Grzywacz D, Rossbach A, Rauf A, et al. (2010) Current control methods for diamondback moth and other brassica insect pests and the prospects for improved management with lepidopteran-resistant Bt vegetable brassicas in Asia and Africa. Crop Prot 29: 68-79.
    [35] Belmain S, Haggar J, Holt J, et al. (2013) Managing legume pests in sub-Saharan Africa: Challenges and prospects for improving food security and nutrition through agro-ecological intensification. Chatham Maritime (United Kingdom): Natural Resources Institute, University of Greenwich, 34.
    [36] Belmain S, Stevenson P (2001) Ethnobotanicals in Ghana: Reviving and modernising age-old farmer practice. Pestic Outlook 12: 233-238.
    [37] Do Ngoc Dai TDT, Ogunmoye A, Eresanya OI, et al. (2015) Chemical constituents of essential oils from the leaves of Tithonia diversifolia, Houttuynia cordata and Asarum glabrum grown in Vietnam. Am J Essent Oils Nat Prod 2: 17-21.
    [38] Grzywacz D, Stevenson P, Mushobozi WL, et al. (2014) The use of indigenous ecological resources for pest control in Africa. Food Secur 6: 71-86.
    [39] Ileke K, Oni M (2011) Toxicity of some plant powders to maize weevil, Sitophilus zeamais (motschulsky)[Coleoptera: Curculiondae] on stored wheat grains (Triticum aestivum). Afr J Biotechnol 6: 3043-3048.
    [40] Kamanula J, Sileshi GW, Belmain SR, et al. (2010) Farmers' insect pest management practices and pesticidal plant use in the protection of stored maize and beans in Southern Africa. Int J Pest Manag 57: 41-49.
    [41] Mkenda PA, Mtei K, Ndakidemi PA (2014) Pesticidal efficacy of Tephrosia vogelii and Tithonia diversifolia against field insect pests of common beans[Phaseolus vulgaris L.]. within African farming communities. Afr J Appl Agric Sci Technol 2: 9-26.
    [42] Mwanauta RW, Mtei K, Ndakidemi PA (2014) Prospective Bioactive Compounds from Vernonia amygdalina, Lippia javanica, Dysphania ambrosioides and Tithonia diversifolia in Controlling Legume Insect Pests. Agric Sci 5: 1129-1139.
    [43] Khater HF (2012) Prospects of botanical biopesticides in insect pest management. Pharmacologia 3: 641-656.
    [44] Marchioro C, Foerster L (2011) Development and survival of the diamondback moth, Plutella xylostella (L.)(Lepidoptera: Yponomeutidae) as a function of temperature: effect on the number of generations in tropical and subtropical regions. Neotrop Entomol 40: 533-541.
    [45] Kobori Y, Amano H (2003) Effect of rainfall on a population of the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae). Appl Entomol Zool 38: 249-253.
  • This article has been cited by:

    1. Ajay Singh Rathore, Vembu Shanthi, A numerical solution of singularly perturbed Fredholm integro-differential equation with discontinuous source term, 2024, 446, 03770427, 115858, 10.1016/j.cam.2024.115858
    2. Abhilipsa Panda, Jugal Mohapatra, Ilhame Amirali, Muhammet Enes Durmaz, Gabil M. Amiraliyev, A numerical technique for solving nonlinear singularly perturbed Fredholm integro-differential equations, 2024, 220, 03784754, 618, 10.1016/j.matcom.2024.02.011
  • Reader Comments
  • © 2021 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(3543) PDF downloads(309) Cited by(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog