The levels of bank capital, risk and efficiency in the Eurozone and the U.S. in the aftermath of the financial crisis

Dimitra Loukia Kolia; Simeon Papadopoulos; Dimitra Loukia Kolia; Simeon Papadopoulos

doi:10.3934/QFE.2020004

Quantitative Finance and Economics

2020, Volume 4, Issue 1: 66-90. doi: 10.3934/QFE.2020004

Previous Article Next Article

Research article

The levels of bank capital, risk and efficiency in the Eurozone and the U.S. in the aftermath of the financial crisis

Department of Accounting and Finance, University of Macedonia, Thessaloniki, Greece

Received: 14 December 2019 Accepted: 06 February 2020 Published: 10 February 2020
JEL Codes: G21

This study investigates the development of the levels of capital, risk and efficiency of the Eurozone and the U.S. banking institutions after the financial crisis. Concerning the methodology, we estimate bank efficiency by applying Data Envelopment Analysis. We estimate bank capital by employing the ratio of the value of total equity to total assets and the Z-score is used as an indicator of bank risk. The findings convey that the efficiency level of the Eurozone banks is considerably lower than that of the U.S. banks. Moreover, the efficiency levels on average increase during the reported period while they reach their peak in the year 2014. Secondly, concerning capitalization, our findings indicate that the capital ratios of the banks of the same sector and different country unions have striking differences with each other. As for the risk ratio, U.S. banks record higher levels of risk than Eurozone banks. This study builds on the existing literature by thoroughly examining bank capital, risk and efficiency with a contemporaneous data set, as the research with data from the period 2013 and onwards is very limited. Additionally, our study is the first to focus on the comparison of U.S. and Eurozone bank samples. The comparison is of utmost importance as the country unions have different characteristics and a different speed of recovery from the financial crisis. We also separately investigate the results per bank type (investment, retail and commercial banks).

Keywords:

Citation: Dimitra Loukia Kolia, Simeon Papadopoulos. The levels of bank capital, risk and efficiency in the Eurozone and the U.S. in the aftermath of the financial crisis[J]. Quantitative Finance and Economics, 2020, 4(1): 66-90. doi: 10.3934/QFE.2020004

Related Papers:

[1]	Hui Xu, Jun Kong, Mengyao Liang, Hui Sun, Miao Qi . Video behavior recognition based on actional-structural graph convolution and temporal extension module. Electronic Research Archive, 2022, 30(11): 4157-4177. doi: 10.3934/era.2022210
[2]	Xite Yang, Ankang Zou, Jidi Cao, Yongzeng Lai, Jilin Zhang . Systemic risk prediction based on Savitzky-Golay smoothing and temporal convolutional networks. Electronic Research Archive, 2023, 31(5): 2667-2688. doi: 10.3934/era.2023135
[3]	Min Li, Ke Chen, Yunqing Bai, Jihong Pei . Skeleton action recognition via graph convolutional network with self-attention module. Electronic Research Archive, 2024, 32(4): 2848-2864. doi: 10.3934/era.2024129
[4]	Jiangtao Zhai, Haoxiang Sun, Chengcheng Xu, Wenqian Sun . ODTC: An online darknet traffic classification model based on multimodal self-attention chaotic mapping features. Electronic Research Archive, 2023, 31(8): 5056-5082. doi: 10.3934/era.2023259
[5]	Bingsheng Li, Na Li, Jianmin Ren, Xupeng Guo, Chao Liu, Hao Wang, Qingwu Li . Enhanced spectral attention and adaptive spatial learning guided network for hyperspectral and LiDAR classification. Electronic Research Archive, 2024, 32(7): 4218-4236. doi: 10.3934/era.2024190
[6]	Huimin Qu, Haiyan Xie, Qianying Wang . Multi-convolutional neural network brain image denoising study based on feature distillation learning and dense residual attention. Electronic Research Archive, 2025, 33(3): 1231-1266. doi: 10.3934/era.2025055
[7]	Chengyong Yang, Jie Wang, Shiwei Wei, Xiukang Yu . A feature fusion-based attention graph convolutional network for 3D classification and segmentation. Electronic Research Archive, 2023, 31(12): 7365-7384. doi: 10.3934/era.2023373
[8]	Zengyu Cai, Liusen Xu, Jianwei Zhang, Yuan Feng, Liang Zhu, Fangmei Liu . ViT-DualAtt: An efficient pornographic image classification method based on Vision Transformer with dual attention. Electronic Research Archive, 2024, 32(12): 6698-6716. doi: 10.3934/era.2024313
[9]	Shaohu Zhang, Jianxiao Ma, Boshuo Geng, Hanbin Wang . Traffic flow prediction with a multi-dimensional feature input: A new method based on attention mechanisms. Electronic Research Archive, 2024, 32(2): 979-1002. doi: 10.3934/era.2024048
[10]	Jiange Liu, Yu Chen, Xin Dai, Li Cao, Qingwu Li . MFCEN: A lightweight multi-scale feature cooperative enhancement network for single-image super-resolution. Electronic Research Archive, 2024, 32(10): 5783-5803. doi: 10.3934/era.2024267

Abstract

1. Introduction

Tungiasis is a parasitic skin condition brought on by female sand fleas and is sometimes referred to as jiggers or sand flea infestation. Tropical and subtropical regions, especially the regions of America, Asia, Africa and the Caribbean, are infected with this disease. The flea enters into the skin, usually in areas where the skin is thin, such as the toes, heels and the spaces between toes and fingers [1]. The female sand flea burrows into the skin and grows into a small, swollen, and painful lesion, often resembling a tiny white or black dot. The flea then lays eggs inside this lesion and after development to complete the life cycle. The primary symptom of tungiasis is intense itching, leading to scratching, which can cause secondary bacterial infections. The affected area may become red, swollen, and painful. Over time, the lesions may grow larger and more uncomfortable [2]. Tungiasis is usually diagnosed based on clinical signs and symptoms. The presence of a specific lesion in the toes or feet, especially in areas with a history of sand flea infestations, helps with diagnosis [3]. The most effective treatment for tungiasis is the removal of the embedded sand fleas using sterilized instruments or a needle. The lesion should be cleaned and disinfected to prevent infections. In severe cases, topical antibiotics may be prescribed if secondary infections have developed [4]. Preventive measures include avoiding walking barefoot in sandy or contaminated areas and wearing protective footwear.

Surgical extraction of embedded sand fleas under sterile medical conditions [5] can also be considered an efficacious therapeutic approach for its management. Tungiasis is a preventable condition, but it can cause significant discomfort and complications if left untreated [6]. If you suspect you have tungiasis or are experiencing symptoms, it is essential to seek medical attention for proper diagnosis and treatment. The infected areas often experience subpar living conditions and inadequate sanitation, making the inhabitants particularly susceptible to the infestation [7]. The afflicted populations face serious health and economical difficulties as a result of this neglected tropical illness. The mechanics of tungiasis transmission must be understood in order to adopt effective control measures and lessen the disease's effects. Jiggers, chigoe flea infection, and pulex penetrans infestation are a few of the historical names for tungiasis [8]. The female sand flea that transmits the disease burrows into the skin, concentrating on the feet, causing painful sores, inflammatory reactions, and secondary infections if left untreated.

Mathematical models have significantly contributed to enhancing comprehension of the fundamental mechanisms of biological processes for public health [9]–[15]. Over the years, researchers and public health experts have dedicated significant efforts to gain insight into the spread of various infections [16]. Previous research has concentrated on a variety of disease-related topics [17], [18], from clinical symptoms and epidemiology to control methods and therapies. The frequency and geographic distribution of tungiasis have been examined in studies, revealing insight on the disease's burden and its effects on afflicted people [19], [20]. Kahuru et al conducted a study focused on employing an optimal control approach within a mathematical model to analyze the dynamics of tungiasis within a community [21]. The researchers in [22] established the stability results of the transmission dynamics of the disease. After that, the theory of optimal control is utilized, all with the ultimate goal of reducing the numbers of infested humans, infested animals, and sand flea populations [23]. The processes of disease transmission and the potential effects of interventions like better sanitation, health education, and treatment campaigns have both benefited greatly from these modeling studies [24]. Mathematical modeling of infectious diseases enhances our understanding of the diseases, facilitates evidence-based decision-making, and helps implement targeted interventions to reduce the burden of the condition in affected communities. In this work, we will enhance our scientific comprehension of tungiasis dynamics to provide more valuable perspectives for efficient management approaches. Through the utilization of mathematical modeling and simulation methods, we will evaluate the prospective consequences of treatments and facilitate decision-making within the realm of public health to improve the effectiveness of interventions and control measures.

Fractional calculus is a branch of mathematics that deals with derivatives and integrals of non-integer orders. It extends the concepts of traditional calculus, which involves integer-order derivatives and integrals. Fractional calculus has found applications in various scientific fields [25]–[30], including epidemiology, where it can be used to model the spread of infectious diseases. It has been acknowledged that fractional calculus offers a more flexible framework for describing complex real-world problems [31], [32]. The Atangana-Baleanu fractional derivative is a type of fractional derivative that has been applied to various mathematical models, including epidemic models [33], [34]. Fractional calculus generalizes the concept of traditional derivatives and integrals to non-integer orders, allowing for a more flexible representation of complex phenomena, such as anomalous diffusion or power-law behaviors [35]. In the context of epidemic models, the Atangana-Baleanu fractional derivative has been utilized to introduce memory effects and long-range interactions into the system. In addition to this, these models can capture more realistic scenarios and provide insights into the long-term behavior of infectious diseases. Hence, we choose to represent the dynamics of tungiasis using a fractional framework to obtain more precise results.

This article is organized as follows: The fundamental concepts and results of the fractional theory related to the ABC operator is presented in Section 2. In Section 3, we formulated a mathematical model for tungiasis through the ABC operator in the fractional framework. In Section 4, the suggested model is subsequently examined. In Section 5, the solution of the recommended model is investigated with the help of fixed-point theory. Then, we introduced a numerical scheme to highlight the dynamical behaviour of the system in Section 6. Finally, ending remarks of the work are presented with future work in Section 7.

2. Preliminaries of fractional-calculus

The main relevant findings and definitions of classical Caputo[36] and Atangana-Baleanu fractional derivatives[37] will be presented here and will be useful in the following portion of the paper.

Definition 2.1. Let $f:[p,q]\rightarrow \mathbb{R}$ be a given function, then the fractional derivative of Caputo fractional derivative of order j is given by

$_{p}^{C} D_{t}^{j} (u (t)) = \frac{1}{Γ (n - j)} \int_{p}^{t} u^{n} (κ) (t - κ)^{n - j - 1} d κ,$

for $j \in (n - 1, n)$ , where $n \in Z$ .

Definition 2.2. The Atangana-Beleanu fractional derivative for a given function f in the Caputo form is introduced as

$_{p}^{A B C} D_{t}^{j} f (t) = \frac{B (j)}{1 - j} \int_{p}^{t} f^{'} (κ) E_{j} [- j \frac{{(t - κ)}^{j}}{1 - j}] d κ,$

where $f \in H^{1} (p, q)$ , q > p, and $j \in [0,1] .$ In addition to this, B(j) is the normalization function which satisfies the condition B(0) = B(1) = 1.

Definition 2.3. The integral of the ABC derivative is denoted by $_{p}^{A B C} I_{t}^{j} f (t)$ and is defined as:

$_{p}^{A B C} I_{t}^{j} f (t) = \frac{1 - j}{B (j)} f (t) + \frac{j}{B (j) Γ (j)} \int_{p}^{t} f (κ) (t - κ)^{j - 1} d κ .$

Theorem 2.1. Assume f in a manner that $f \in C [p, q]$ , then we have the equation below mentioned in [37]:

$\|{}^{{ABC}}_{p}D^{\jmath}_{t}(f(t))\| < \frac{B(\jmath)}{1-\jmath}\|f(t)\|, \;where\; \|f(t)\|=max_{p\leq t\leq q}|f(t)|.$

Furthermore, it has been established that ABC derivatives satisfy the Lipschitz condition [37].

$\|{}^{{ABC}}_{p}D^{\jmath}_{t}f_1(t)-{}^{{ABC}}_{p}D^{\jmath}_{t}f_2(t)\| < \mathcal{\vartheta}_1 \|f_1(t)-f_2(t)\|.$

3. Evaluation of the dynamics

Here we present a mathematical model of Tungiasis disease. In this approach, the total population is divided into four categories: group that practices good hygiene ( $\mathcal{P}$ ), susceptible ( $\mathcal{S}$ ), infected ( $\mathcal{I}$ ), and treated ( $\mathcal{T}$ ) groups. Then, the total population is

$\mathcal{N}(t)= \mathcal{P}(t)+\mathcal{S}(t)+\mathcal{I}(t)+\mathcal{T}(t).$

At birth, the population that is most susceptible is recruited at a rate of $(1 - γ) υ$ , whereas the group that practices good hygiene is recruited at γυ where υ represents the rate of recruitment at birth and γ represents the likelihood of being recruited into the category of good hygiene practices. The appropriate hygiene practice group ( $\mathcal{P}$ ) develops susceptibility ( $\mathcal{S}$ ) at a rate of φ. Individuals from the class ( $\mathcal{S}$ ) move to the class ( $\mathcal{I}$ ) with a rate Г. The infected individuals ( $\mathcal{I}$ ) transfer to the treated group ( $\mathcal{T}$ ) at a rate of ρ following therapy. We indicated the natural death rate by µ. The infection rate Г is defined as

$\Gamma=\dfrac{\sigma\varepsilon\mathcal{I}}{\mathcal{N}},$

where ϵ is the contact rate with infected people and σ is the likelihood of being infected after repeated exposure to infected people. Then our model with the above assumptions for tungiasis disease is as follows:

$\begin{equation} \left\{\begin{array}{rcl} \frac{d\mathcal{P}}{dt} &=& \gamma\upsilon-\left(\varphi+\mu\right)\mathcal{P}, \\ \frac{d\mathcal{S}}{dt} &=& \left(1-\gamma\right)\upsilon+\varphi\mathcal{P}-\left(\Gamma+\mu\right)\mathcal{S} , \\ \frac{d\mathcal{I}}{dt} &=& \Gamma\mathcal{S}-\left(\omega+\varrho+\mu\right)\mathcal{I} , \\ \frac{d\mathcal{T}}{dt} &=& \varrho\mathcal{I}-\mu\mathcal{T} , \end{array}\right. \end{equation}$

(3.1)

with appropriate initial condition

$\begin{equation} \begin{array}{rcl} \mathcal{P}(0) > 0,\ \mathcal{S}(0) > 0,\ \mathcal{I}(0) > 0,\ \mathcal{T}(0) > 0. \end{array} \end{equation}$

(3.2)

Table 1. Parameter values with interpretation used in the proposed system.

Parameter	Interpretation
γ	Probability of recruitment in $\mathcal{P}(t)$
υ	Birth or recruitment rate
φ	The losing rate of protection
µ	Death occurs naturally in each class
Г	Infection rate from susceptible to infected class
ω	Disease induced mortality rate
ρ	Treatment rate for the infection
σ	Rate of transmission of the infection
ϵ	Contact rate of infection

| Show Table

DownLoad: CSV

It is acknowledged that epidemic models through the Atangana-Baleanu fractional derivative can more accurately represents the dynamics of infectious diseases than conventional integer derivatives. This inclusion enables a more comprehensive exploration of epidemic spread in intricate settings characterized by long-range interactions and memory effects, fostering an improved comprehension of the underlying mechanisms. Consequently, the refined insights gained from these fractional derivative-based models hold the potential to enhance the precision of epidemic predictions and facilitate the development of more efficacious control strategies. Nonetheless, it is crucial to acknowledge that the introduction of fractional derivatives into mathematical models may introduce increased system complexity, necessitating the utilization of specialized mathematical and numerical techniques for rigorous analysis and simulation. Thus, the model (3.1) recommended above can be written in fractional form as

$\begin{equation} \left\{\begin{array}{rcl} ^{ABC}_0 D^{\jmath}_t \mathcal{P} &=& \gamma\upsilon-\left(\varphi+\mu\right)\mathcal{P} , \\ ^{ABC}_0 D^{\jmath}_t \mathcal{S} &=& \left(1-\gamma\right)\upsilon+\varphi\mathcal{P}-\left(\Gamma+\mu\right)\mathcal{S} ,\\ ^{ABC}_0 D^{\jmath}_t \mathcal{I} &=& \Gamma\mathcal{S}-\left(\omega+\varrho+\mu\right)\mathcal{I} , \\ ^{ABC}_0 D^{\jmath}_t \mathcal{T} &=& \varrho\mathcal{I}-\mu\mathcal{T} . \end{array}\right. \end{equation}$

(3.3)

All of the parameters in the recommended model (3.1) for tungiasis disease are considered to be positive. Furthermore, the term $_{0}^{A B C} D_{t}^{j}$ in the above model indicates the ABC derivative. The effects of memory which arise in the epidemiological process but are not incorporated in the classical operators, are the driving force behind fractional order models. In Table 1, we illustrated the parameters of the recommended model with description.

4. Analysis of the model

In this section, we deal with some analysis of the model, like steady-state analysis, and find the basic reproduction ratio $\mathcal{R}_0$ for the above model (3.3). As the system (3.3) models the population of humans, all state variable solutions with nonnegative beginning conditions are nonnegative $\forall t > 0$ and have a feasible region bound as follows:

$\Xi=\left(\left(\mathcal{P}, \mathcal{S}, \mathcal{I}, \mathcal{T}\right)\in\mathbb{R}_+^{4}; \mathcal{S}, \mathcal{P},\mathcal{I},\mathcal{T}\geq 0; \mathcal{N}\leq\frac{\upsilon}{\mu}\right).$

In the upcoming step, we will focus on the steady-states of the suggested system of the infection. First, we take the following:

$^{ABC}_0 D^{\jmath}_t \mathcal{P} = ^{ABC}_0 D^{\jmath}_t \mathcal{S} = ^{ABC}_0 D^{\jmath}_t \mathcal{I} = ^{ABC}_0 D^{\jmath}_t \mathcal{T} = 0,$

then, the model (3.3) becomes

$\begin{equation} \left\{\begin{array}{rcl} 0 &=& \gamma\upsilon-\left(\varphi+\mu\right)\mathcal{P} , \\ 0 &=& \left(1-\gamma\right)\upsilon+\varphi\mathcal{P}-\left(\Gamma+\mu\right)\mathcal{S} ,\\ 0 &=& \Gamma\mathcal{S}-\left(\omega+\varrho+\mu\right)\mathcal{I} , \\ 0 &=& \varrho\mathcal{I}-\mu\mathcal{T} . \end{array}\right. \end{equation}$

(4.1)

The disease-free equilibrium of an epidemic model refers to a stable state where the infectious disease has been eliminated from the population. Understanding the disease-free equilibrium is essential in assessing the effectiveness of control strategies and vaccination programs. It serves as a benchmark to measure the success of public health interventions in eliminating or controlling infectious diseases and preventing outbreaks. For the infection-free steady-state, we set

$\mathcal{P}=\mathcal{P}^o, \mathcal{S}=\mathcal{S}^o, \mathcal{I}=\mathcal{I}^o, \mathcal{T}=\mathcal{T}^o \ \ and \ \ \mathcal{N}=\mathcal{N}^o.$

Thus, we have

$\mathit{E}_o = \left(\mathcal{P}^o, \mathcal{S}^o, \mathcal{I}^o, \mathcal{T}^o \right) = \left(0, \frac{\left(1-\gamma\right)\upsilon}{\left(\Gamma+\mu\right)},0, 0\right).$

The endemic equilibrium (EE) of an epidemic model refers to a stable state where the prevalence of the infectious disease remains constant over time. Understanding the endemic equilibrium of an epidemic model is crucial in assessing the long-term behavior of the infection and evaluating the quality of control measures in maintaining the disease at a manageable level. It provides valuable insights into the stability and persistence of the infection in a population and is a fundamental concept in the study of infectious disease dynamics. We indicate the EE of the system by (E*) and substitute $\Gamma=\frac{\sigma\varepsilon\mathcal{I}}{\mathcal{N}}$ . Then from (4.1), we have

$\begin{equation} \left\{\begin{array}{rcl} 0 &=& \gamma\upsilon-\left(\varphi+\mu\right)\mathcal{P}^* , \\ 0 &=& \left(1-\gamma\right)\upsilon+\varphi\mathcal{P}^*-\left(\frac{\sigma\varepsilon\mathcal{I}^*}{\mathcal{N}^*}+\mu\right)\mathcal{S}^* ,\\ 0 &=& \frac{\sigma\varepsilon\mathcal{I}^*}{\mathcal{N}^*}\mathcal{S}^*-\left(\omega+\varrho+\mu\right)\mathcal{I}^* , \\ 0 &=& \varrho\mathcal{I}^*-\mu\mathcal{T}^*. \end{array}\right. \end{equation}$

(4.2)

Solving the above system (4.2), we have

$\mathit{E}^*\left( \begin{array}{cccc} \mathcal{P}^*\\ \mathcal{S}^*\\ \mathcal{I}^*\\ \mathcal{T}^* \end{array}\right) = \mathit{E}^*\left( \begin{array}{cccc} \frac{\gamma\upsilon}{\varphi+\mu}\\ \frac{\mathcal{N}^*\left(\omega+\varrho+\mu\right)}{\sigma\varepsilon}\\ \frac{1}{\left(\omega+\varrho+\mu\right)}\left[\frac{(\varphi-\gamma\mu+\mu)\upsilon}{(\varphi+\mu)} - \frac{\mathcal{N}^*\mu(\omega+\varrho+\mu)}{\sigma\varepsilon}\right]\\ \frac{\varrho}{\mu(\omega+\varrho+\mu)}\left[\frac{(\varphi-\gamma\mu+\mu)\upsilon}{(\varphi+\mu)} - \frac{\mathcal{N}^*\mu(\omega+\varrho+\mu)}{\sigma\varepsilon}\right] \end{array}\right).$

The basic reproduction number is a fundamental concept in epidemiology that quantifies the potential of an infection to spread. It represents the average number of secondary infections generated by a single infected individual in a fully susceptible environment. If the basic reproduction number is less than 1, it means that, on average, each infected individual will cause fewer than one new infection during their infectious period. In this case, the disease will not be able to sustain itself in the population, and it will eventually die out. If this parameter is equal to 1, this implies that, on average, each infected individual will cause exactly one new infection during their infectious period. The disease will reach a stable endemic equilibrium, where the number of new infections is balanced by the number of recoveries, resulting in a constant prevalence of the disease in the population. If this parameter is greater than 1, this implies that each infected individual, on average, will cause more than one new infection during their infectious period. In this case, the disease has the potential to spread, and if left uncontrolled, it may lead to an epidemic outbreak.

We symbolized this number as $\mathcal{R}_0$ and calculate it through the next-generation matrix method [38]. We take the infected class of the system (3.1), so we have

$\begin{equation} ^{ABC}_0 D^{\jmath}_t \mathcal{I} = \Gamma\mathcal{S}-\left(\omega+\varrho+\mu\right)\mathcal{I}. \end{equation}$

(4.3)

Let F be the matrix containing parameters that are entering into infected class and V be the matrix containing those parameters that leave infected classes while neglecting the negative signs, so we get

$F= \Gamma\mathcal{S} = \frac{\sigma\varepsilon\mathcal{I}\mathcal{S}}{\mathcal{N}},$

and

$V= \left(\omega+\varrho+\mu\right)\mathcal{I}.$

This implies that $\mathcal{R}_0= \dfrac{\sigma\varepsilon}{\omega+\varrho+\mu}$ , which is the required $\mathcal{R}_0$ of the recommended system.

Theorem 4.1. For j = 1, the infection-free steady-state E₀ of the system (3.3) is locally asymptotically stable if $\mathcal{R}_0<1$ , otherwise it is unstable.

Proof. Let j = 1 and take the Jacobian matrix of the model (3.3) at infection-free steady-state E₀ as

$\mathfrak{J}(E_o)=\begin{bmatrix} -(\varphi+\mu) &0 &0 &0 \\ \varphi &-\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}+\mu\right) &-\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &0 \\ 0 &\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &\frac{\sigma\varepsilon\mathcal{S}^o}{\mathcal{N}}- \left(\varrho+\omega+\mu\right) &0 \\ 0 &0 &\varrho &-\mu \end{bmatrix}.$

The above Jacobian matrix has the following characteristic equation:

$\left|\mathfrak{J}(E_o)-\lambda I\right| = 0,$

furthermore, we have

$\begin{vmatrix} -(\varphi+\mu)-\lambda &0 &0 &0 \\ \varphi &-\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}+\mu\right)-\lambda &-\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &0 \\ 0 &\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &\frac{\sigma\varepsilon\mathcal{S}^o}{\mathcal{N}}- \left(\varrho+\omega+\mu\right)-\lambda &0 \\ 0 &0 &\varrho &-\mu-\lambda \end{vmatrix}=0,$

and simplifying, we have

$\left(-(\varphi+\mu)-\lambda\right)\begin{vmatrix} -\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}+\mu\right)-\lambda &-\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &0 \\ \frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &\frac{\sigma\varepsilon\mathcal{S}^o}{\mathcal{N}}- \left(\varrho+\omega+\mu\right)-\lambda &0 \\ 0 &\varrho &-\mu-\lambda \end{vmatrix}=0,$

which implies that $λ_{1} = - (ϕ + µ)$ and

$\begin{vmatrix} -\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}+\mu\right)-\lambda &-\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &0 \\ \frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &\frac{\sigma\varepsilon\mathcal{S}^o}{\mathcal{N}}- \left(\varrho+\omega+\mu\right)-\lambda &0 \\ 0 &\varrho &-\mu-\lambda \end{vmatrix}=0,$

the simplification of which implies that

$\left(-\mu-\lambda\right)\begin{vmatrix} -\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}+\mu\right)-\lambda &-\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} \\ \frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &\frac{\sigma\varepsilon\mathcal{S}^o}{\mathcal{N}}- \left(\varrho+\omega+\mu\right)-\lambda \end{vmatrix}=0.$

From the above, we have $λ_{2} = - µ$ and

$\begin{vmatrix} -\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}+\mu\right)-\lambda &-\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} \\ \frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &\frac{\sigma\varepsilon\mathcal{S}^o}{\mathcal{N}}- \left(\varrho+\omega+\mu\right)-\lambda \end{vmatrix}=0$

Hence, all the eigen values are negative. Furthermore, we have $\mathcal{J}_1(E_0)$ given by

$\mathcal{J}_1(E_0)=\begin{bmatrix} -\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}+\mu\right) &-\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} \\ \frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &\frac{\sigma\varepsilon\mathcal{S}^o}{\mathcal{N}}- \left(\varrho+\omega+\mu\right) \end{bmatrix}=0.$

To show the required result, we have to show that Tace( $\mathcal{J}_1(E_o))<0$ and Det( $\mathcal{J}_1(E_o))>0$ , so

$Trace(\mathcal{J}_1(E_o))=-\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}+\mu\right)+\frac{\sigma\varepsilon\mathcal{S}^o}{\mathcal{N}}- \left(\varrho+\omega+\mu\right),$

$=-\mu+\sigma\varepsilon - \left(\varrho+\omega+\mu\right) ,$

$=-\mu+\sigma\varepsilon -\varrho-\omega-\mu ,$

$=-2\mu+\sigma\varepsilon -(\varrho+\omega) ,$

$Trace(\mathcal{J}_1(E_0))=-\left(2\mu-\sigma\varepsilon+(\varrho+\omega)\right) < 0,$

which implies that

$Trace(\mathcal{J}_1(E_o)) < 0.$

Also, we have

$Det\left(\mathcal{J}_1(E_0)\right)=\begin{vmatrix} -\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}+\mu\right) &-\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} \\ \frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}} &\frac{\sigma\varepsilon\mathcal{S}^o}{\mathcal{N}}- \left(\varrho+\omega+\mu\right) \end{vmatrix}$

$=\left[-\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}+\mu\right)\left(\frac{\sigma\varepsilon\mathcal{S}^o}{\mathcal{N}}- \left(\varrho+\omega+\mu\right)\right) - \left(-\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}\right)\left(\frac{\sigma\varepsilon\mathcal{I}^o}{\mathcal{N}}\right) \right].$

Because $\mathcal{R}_0<1$ is given, so the Det will be positive, that is, $Det(\mathcal{J}_1(E_0))>0.$ Hence, the infection-free steady-state of the recommended system is locally asymptotically stable (LAS) for $\mathcal{R}_0<1$ and unstable in other cases.

Theorem 4.2. For j = 1, the endemic steady-state (E*) of system (3.3) is LAS if $\mathcal{R}_0>1$ and is unstable in other circumstances.

Proof. Let j = 1 and take the Jacobian matrix of system (3.3) at endemic steady-state (E*) as

$\mathfrak{J}(E)^*=\begin{bmatrix} -(\varphi+\mu) &0 &0 &0 \\ \varphi &-\left(\frac{\sigma\varepsilon\mathcal{I}^*}{\mathcal{N}^*}+\mu\right) &-\frac{\sigma\varepsilon\mathcal{I}^*}{\mathcal{N}^*} &0 \\ 0 &\frac{\sigma\varepsilon\mathcal{I}^*}{\mathcal{N}^*} &\frac{\sigma\varepsilon\mathcal{S}^*}{\mathcal{N}^*}- \left(\varrho+\omega+\mu\right) &0 \\ 0 &0 &\varrho &-\mu, \end{bmatrix},$

furthermore, we get

$\mathfrak{J}(E^*)=\begin{bmatrix} -(\varphi+\mu) &0 &0 &0 \\ \varphi &-\left(\frac{\sigma\varepsilon W}{\mathcal{N}^*}+\mu\right) &-\frac{\sigma\varepsilon W}{\mathcal{N}^*}&0 \\ 0 &\frac{\sigma\varepsilon W}{\mathcal{N}^*} &-(\varrho+\omega+\mu) &0 \\ 0 &0 &\varrho &-\mu, \end{bmatrix},$

in which $W= \frac{1}{\left(\omega+\varrho+\mu\right)}\left[\frac{(\varphi-\gamma\mu+\mu)\upsilon}{(\varphi+\mu)} - \frac{\mathcal{N}^*\mu(\omega+\varrho+\mu)}{\sigma\varepsilon}\right].$ The characteristic equation is $(\varphi+\mu+\lambda)(\mu+\lambda)\left[\lambda^2+ \left(\frac{\sigma\varepsilon W}{\mathcal{N}^*}+\mu\right)\lambda+ \frac{\sigma\varepsilon(\omega+\varrho+\mu) W}{\mathcal{N}^*} \right]=0,$ which implies that $λ_{1} = - (ϕ + µ)$ and $λ_{2} = - µ$ with the following

$\begin{equation} \lambda^2+ \left(\frac{\sigma\varepsilon W}{\mathcal{N}^*}+\mu\right)\lambda+ \frac{\sigma\varepsilon(\omega+\varrho+\mu) W}{\mathcal{N}^*} = 0. \end{equation}$

(4.4)

Substituting the value of W into (4.4), we get

$\begin{equation} \lambda^2+m\lambda+n=0, \end{equation}$

(4.5)

in which $m=\frac{\mathcal{R}_o((\varphi-\gamma\mu+\mu)\upsilon)}{(\varphi+\mu)\mathcal{N}^*}$ and $n=\frac{\sigma\varepsilon(\varphi-\gamma\mu+\mu)\upsilon}{(\varphi+\mu)\mathcal{N}^*} - (\omega+\varrho+\mu)\mu .$ Hence it is clear that $m, n > 1$ , also $\mathcal{R}_0=\frac{\sigma\varepsilon}{\omega+\varrho+\mu} > 1.$ The eigenvalues provided by equation (4.5) are negative according to the Routh-Hurwitz criteria. Thus, whenever $\mathcal{R}_0 > 1$ , the EE of the proposed model of tungiasis is locally asymptotically stable.

5. Fractional order model solution

In this section of the paper, we utilize fixed point theory to demonstrate the existence and uniqueness of the fractional order model (3.3) solution. The system of equations can be represented in the following way:

$\begin{eqnarray} \left\{ \begin{array}{ll} ^{ABC}_0 D^{\jmath}_{t}s(t)=\mathcal{U}(t,v(t)),\\ v(0)=v_0, \;\; 0 < t < T < \infty. \end{array} \right. \end{eqnarray}$

(5.1)

In the above (5.1), $v(t)=(\mathcal{P}, \mathcal{S}, \mathcal{I}, \mathcal{T})$ is the vector form of the state variables and $\mathcal{U}$ is a continuous vector function. Moreover, $\mathcal{U}$ is given by

$\mathcal{U}=\left( \begin{array}{ccccc} \mathcal{U}_1\\ \mathcal{U}_2\\ \mathcal{U}_3\\ \mathcal{U}_4 \end{array} \right)=\left( \begin{array}{cccc} \gamma\upsilon-\left(\varphi+\mu\right)\mathcal{P} , \\ \left(1-\gamma\right)\upsilon+\varphi\mathcal{P}-\left(\Gamma+\mu\right)\mathcal{S} ,\\ \Gamma\mathcal{S}-\left(\omega+\varrho+\mu\right)\mathcal{I} , \\ \varrho\mathcal{I}-\mu\mathcal{T} \end{array} \right),$

and $v_0(t)=\left(\mathcal{P}(0), \mathcal{S}(0), \mathcal{I}(0), \mathcal{T}(0)\right)$ are initial conditions. In addition to this, the Lipschitz condition is fulfilled by $\mathcal{U}$ as

$\begin{eqnarray} \|\mathcal{U}(t, v_1(t))-\mathcal{U}(t, v_2(t))\|&\leq&\mathcal{N}\|v_1(t)-v_2(t)\|. \end{eqnarray}$

(5.2)

In the upcoming step, the existence and uniqueness of the fractional order dynamical model (3.3) will be proved.

Theorem 5.1. The solution of fractional order (FO) model of equations (3.3) will be unique if the below mentioned condition satisfies:

$\begin{equation} \frac{(1-\jmath)}{ABC(\jmath)}\mathcal{N}+\frac{\jmath}{ABC(\jmath)\Gamma(\jmath)}T^{\jmath}_{max}\mathcal{N} < 1. \end{equation}$

(5.3)

Proof. In order to demonstrate the desired result, we utilize the Atangana-Beleanu (AB) fractional integral, as defined in (2.3), on system (5.1). This yields a non-linear Volterra integral equation as follows:

$\begin{eqnarray} v(t)=v_0+\frac{1-\jmath}{ABC(\jmath)}\mathcal{U}(t,r(t))+\frac{\jmath}{ABC(\jmath)\Gamma(\jmath)}\int_0^t(t-\kappa)^{\jmath-1}\mathcal{U}(\kappa,v(\kappa))d\kappa. \end{eqnarray}$

(5.4)

Take $I=(0, T)$ and the operator $\Phi:\mathcal{G}(I,R^4)\rightarrow \mathcal{G}(I,R^4)$ defined by

$\begin{eqnarray} \Phi[v(t)]=v_0+\frac{1-\jmath}{ABC(\jmath)}\mathcal{U}(t,v(t))+\frac{\jmath}{ABC(\jmath)\Gamma(\jmath)}\int_0^t(t-\kappa)^{\jmath-1}\mathcal{U}(\kappa,v(\kappa))d\kappa. \end{eqnarray}$

(5.5)

Equation (5.4) takes the following form

$\begin{eqnarray} v(t)=\Phi[v(t)], \end{eqnarray}$

(5.6)

and the supremum norm on I is indicated by $\|.\|_I$ and is

$\begin{eqnarray} \|v(t)\|_I=\sup_{t\in I}\|v(t)\|,\;\;v(t)\in\mathcal{G}. \end{eqnarray}$

(5.7)

Certainly, $\mathcal{G}(I,R^4)$ with norm $\|.\|_I$ becomes a Banach space and

$\begin{eqnarray} \bigg\|\int_0^t \mathcal{F} (t,\kappa)v(\kappa)d\kappa\bigg\|\leq T\|\mathcal{F} (t,\kappa)\|_I \|v(t)\|_I, \end{eqnarray}$

(5.8)

with $v(t)\in \mathcal{G}(I,R^4)$ , $\mathcal{F}(t,\kappa)\in \mathcal{G}(I^2,R)$ in a manner that

$\begin{eqnarray} \|\mathcal{F}(t,\kappa)\|_I=\sup_{t,\kappa \in I}|\mathcal{F}(t,\kappa)|. \end{eqnarray}$

(5.9)

Here, utilizing the definition of Φ shown in (5.6), we have

$\begin{eqnarray} \|\Phi [v_1(t)]-\Phi [v_2(t)]\|_I&\leq& \bigg \| \frac{(1-\jmath)}{ABC(\jmath)} (\mathcal{U} (t,v_1(t))-\mathcal{U}(t,v_2(t))+\frac{\jmath}{ABC(\jmath)\Gamma(\jmath)}\times\cr\cr&&\int^t_0(t-\kappa)^{\jmath-1}(\mathcal{U}(\kappa,v_1(\kappa))-\mathcal{U}(\kappa,v_2(\kappa)))d\kappa\bigg\|_I \end{eqnarray}$

(5.10)

Additionally, by applying the Lipschitz condition (5.2) and the triangular inequality, in combination with the result in (5.8), we arrive at the below expression after simplification:

$\begin{eqnarray} \|\Phi [v_1(t)]-\Phi [v_1(t)]\|_I&\leq& \bigg(\frac{(1-\jmath)\mathcal{N}}{ABC(\jmath)}+\frac{\jmath}{ABC(\jmath)\Gamma(\jmath)}\mathcal{N}T^{\jmath}_{max}\bigg)\| v_1(t)-v_2(t)\|_I. \end{eqnarray}$

(5.11)

Consequently, we have the following:

$\begin{eqnarray} \|\Phi [v_1(t)]-\Phi [v_1(t)]\|_I&\leq& B\| v_1(t)-v_2(t)\|_I, \end{eqnarray}$

(5.12)

in which

$B=\frac{(1-\jmath)\mathcal{N}}{ABC(\jmath)}+\frac{\jmath}{ABC(\jmath)\Gamma(\jmath)}\mathcal{N}T^{\jmath}_{max} .$

It is evident that when condition (5.3) is met, Φ becomes a contraction, implying that the fractional order dynamical system (5.1) possesses a unique solution.

6. Iterative scheme and results

In this section, we will present an iterative scheme for the numerical solution of the recommended model (3.3) of the infection. First, an iterative scheme will be developed and then the the solution pathways of the system will be represented through the scheme to understand the dynamics. We use the newly established numerical approach (iterative) proposed for the approximation of the AB integral operator [39]. We briefly describe and apply the aforesaid approach to our dynamical system (3.3) in order to show the impact of different parameters on the infection.

Here, rewriting system (5.4) of the infection into the fractional integral equation form using the fundamental theorem of fractional calculus:

$\begin{eqnarray} s(t)-s(0)&=& \frac{(1-\jmath)}{ABC(\jmath)} \mathcal{U}(t,s(t)) +\frac{\jmath}{ABC(\jmath)\times\Gamma(\jmath)}\int^t_0\mathcal{U}(\kappa,x(\xi))(t-\kappa)^{\jmath-1}d\xi. \end{eqnarray}$

(6.1)

At $t = t_{ξ + 1}$ , $ξ = 0,1,2,...,$ we have

$\begin{eqnarray}\label{N1} s(t_{\xi+1})-s(0)&=&\frac{1-\jmath}{ABC(\jmath)}\mathcal{U}(t_\xi,s(t_\xi))+\cr\cr&&\frac{\jmath}{ABC(\jmath)\times\Gamma(\jmath)}{\int_0^{{t}_{\xi+1}}} \mathcal{U}(\kappa,s(\kappa))(t_{\xi+1}-\kappa)^{\jmath-1}d\jmath,\cr\cr &=&\frac{1-\jmath}{ABC(\jmath)}\mathcal{U}( t_\xi,s(t_\xi))+\cr\cr&&\frac{\jmath}{ABC(\jmath)\times\Gamma(\jmath)}\sum^{\xi}_{j=0}{\int_{t_j}^{{t}_{j+1}}} \mathcal{U}(\kappa,s(\kappa))(t_{\xi+1}-\kappa)^{\jmath-1}d\jmath . \end{eqnarray}$

(6.2)

The function $\mathcal{U} (\kappa,s(\kappa))$ can be estimated over the interval $[t_{j}, t_{j + 1}]$ , and we apply the interpolation polynomial

$\begin{eqnarray} \mathcal{U} (\kappa,s(\kappa)) &\cong &\frac{\mathcal{U}(t_j, s(t_j) )}{\mathrm{h}} (t-t_{j-1} )-\frac{\mathcal{U}(t_{j-1}, s(t_{j-1})) }{\mathrm{h}} (t-t_j ), \end{eqnarray}$

(6.3)

then substituting in (6.2), we get

$\begin{eqnarray} s(t_{\xi+1})&=&s(0)+\frac{1-\jmath}{ABC(\jmath)}\mathcal{U}(t_\xi,s(t_\xi))+\frac{\jmath}{ABC(\jmath)\times\Gamma(\jmath)}\cr\cr&&\sum^{\xi}_{j=0}\Big(\frac{\mathcal{U}(t_j,s(t_j))}{\mathrm{h}}{\int^{t_{j+1}}_{{t_j}}} (t-t_{j-1})(t_{\xi+1}-t)^{\jmath-1}dt\cr\cr&&-\frac{\mathcal{U}( t_{j-1},s(t_{j-1}))}{\mathrm{h}}{\int^{t_{j+1}}_{{t_j}}}(t-t_j)(t_{\xi+1}-t)^{\jmath-1}d t\Big). \end{eqnarray}$

(6.4)

After computing these integrals, the approximate answer is as follows:

$\begin{eqnarray} s(t_{\xi+1})&=&s(t_0)+\frac{1-\jmath}{ABC(\jmath)}\mathcal{U}( t_\xi,s(t_\xi))+\frac{\jmath}{ABC(\jmath)}\sum^{\xi}_{j=0}\cr\cr&&\Big(\frac{\mathrm{h}^{\jmath}\mathcal{U}( t_j,s(t_j))}{\Gamma(\jmath+2)}((\xi+1-j)^\jmath(\xi-j+2+\jmath)-(\xi-j)^{\jmath}(\xi-j+2+2\jmath))\cr\cr&&-\frac{\mathrm{h}^\jmath \mathcal{U}(t_{j-1},s(t_{j-1}))}{\Gamma(\jmath+2)}((\xi+1-j)^{\jmath+1}-(\xi-j)^{\jmath}(\xi-j+1+\jmath))\Big). \end{eqnarray}$

(6.5)

Finally, for the suggested model, we achieved the following recursive formulae:

$\begin{eqnarray} \mathcal{P}(t_{\xi+1}) &=& \mathcal{P}(t_0)+\frac{1-\jmath}{ABC(\jmath)}\mathcal{U}_1( t_\xi,s(t_\xi))+\frac{\jmath}{ABC(\jmath)}\sum^{\xi}_{j=0}\cr\cr&&\Big(\frac{\mathrm{h}^{\jmath}\mathcal{U}_1( t_j,s(t_j))}{\Gamma(\jmath+2)}((\xi+1-j)^\jmath(\xi-j+2+\jmath)-(\xi-j)^{\jmath}(\xi-j+2+2\jmath))\cr\cr&&-\frac{\mathrm{h}^\jmath \mathcal{U}_1(t_{j-1},s(t_{j-1}))}{\Gamma(\jmath+2)}((\xi+1-j)^{\jmath+1}-(\xi-j)^{\jmath}(\xi-j+1+\jmath))\Big) \cr\cr \mathcal{S}(t_{\xi+1})&=&\mathcal{S}(t_0)+\frac{1-\jmath}{ABC(\jmath)}\mathcal{U}_2( t_\xi,s(t_\xi))+ \frac{\jmath}{ABC(\jmath)}\sum^{k}_{j=0}\cr\cr&&\Big(\frac{f^{\jmath}\mathcal{U}_2( t_j,s(t_j))}{\Gamma(\jmath+2)}((\xi+1-j)^\jmath(\xi-j+2+\jmath)-(\xi-j)^{\jmath}(\xi-j+2+2\jmath))\cr\cr&&-\frac{\mathrm{h}^\jmath \mathcal{U}_2(t_{j-1},s(t_{j-1}))}{\Gamma(\jmath+2)}((\xi+1-j)^{\jmath+1}-(\xi-j)^{\jmath}(\xi-j+1+\jmath))\Big) \cr\cr \mathcal{I}(t_{\xi+1})&=&\mathcal{I}(t_0)+\frac{1-\jmath}{ABC(\jmath)}\mathcal{U}_3( t_\xi,s(t_\xi))+ \frac{\jmath}{ABC(\jmath)}\sum^{\xi}_{j=0}\cr\cr&&\Big(\frac{\mathrm{h}^{\jmath}\mathcal{U}_3( t_j,s(t_j))}{\Gamma(\jmath+2)}((\xi+1-j)^\jmath(\xi-j+2+\jmath)-(\xi-j)^{\jmath}(\xi-j+2+2\jmath))\cr\cr&&-\frac{\mathrm{h}^\jmath \mathcal{U}_3(t_{j-1},s(t_{j-1}))}{\Gamma(\jmath+2)}((\xi+1-j)^{\jmath+1}-(\xi-j)^{\jmath}(\xi-j+1+\jmath))\Big) \cr\cr \mathcal{T}(t_{\xi+1})&=&\mathcal{T}(t_0)+\frac{1-\jmath}{ABC(\jmath)}\mathcal{U}_4( t_{\xi},s(t_\xi))+\frac{\jmath}{ABC(\jmath)}\sum^{\xi}_{j=0}\cr\cr&&\Big(\frac{f^{\jmath}\mathcal{U}_4( t_j,s(t_j))}{\Gamma(\jmath+2)}((\xi+1-j)^\jmath(\xi-j+2+\jmath)-(\xi-j)^{\jmath}(\xi-j+2+2\jmath))\cr\cr&&-\frac{f^\jmath \mathcal{U}_4(t_{j-1},s(t_{j-1}))}{\Gamma(\jmath+2)}((n+1-j)^{\jmath+1}-(\xi-j)^{\jmath}(\xi-j+1+\jmath))\Big). \end{eqnarray}$

(6.6)

Figure 1. Visualization of the tracking path behavior of the suggested model under varying fractional parameter values j, i.e., j = 0.5, 0.6, 0.7, 0.8.

DownLoad: Full-Size Img PowerPoint

Figure 2. Illustration of the solution pathways of the suggested model with varying values of the losing rate of protection φ, i.e., φ = 0.40, 0.45, 0.50, 0.55.

DownLoad: Full-Size Img PowerPoint

Figure 3. Illustration of the solution pathways of the suggested model with varying values of the treatment rate ρ, i.e., ρ = 0.02, 0.04, 0.06, 0.08.

DownLoad: Full-Size Img PowerPoint

Figure 4. Representing the time series of the class of the proposed system of the infection with varying values of the transmission probability σ, i.e., σ = 0.40, 0.50, 0.60, 0.70.

DownLoad: Full-Size Img PowerPoint

Figure 5. Visualization of the tracking path behavior of the suggested model under varying fractional parameter values j, i.e., j = 0.85, 0.90, 0.95, 1.00.

DownLoad: Full-Size Img PowerPoint

We will utilize the above alterative approach to investigate the dynamical behaviour of the recommended fractional model of tungiasis. It is well-known that the graphical view analysis of epidemic models provide a comprehensive exploration of the model's behavior, leading to a better understanding of the epidemic's spread and potential control strategies. Here, we will perform different simulations to visualize the impact of different input factors on the system.

The values of the input parameter and state variables of the system are assumed for numerical purposes. For our investigation, we conducted a series of simulations to analyze the behavior of the system under different conditions. In these simulations, we varied specific parameters and observed their impact on the solution pathways of the system. First, in Figure 1 and Figure 5, we examined the effect of the fractional parameter j on the dynamics of the system. By changing the fractional order, we aimed to understand its role in shaping the infection spread in the society. Notably, we observed that decreasing the fractional parameter led to a reduction in the infection rate within the community. This finding indicates that the fractional parameter can serve as a crucial control parameter for managing the epidemic. In Figure 2, we focused on the losing rate of protection and its influence on the dynamics of the tungiasis infection. The results clearly indicate that this parameter plays a pivotal role in determining the risk of infection within the population. A higher losing rate of protection was found to increase the susceptibility to infection, posing a potential threat to public health.

The treatment rate was the primary parameter of interest in Figure 3. By varying the treatment rate ρ, we explored its impact on the transmission dynamics of the system. The simulation revealed that the treatment rate significantly affects the control of the infection. Higher values of ρ demonstrated a more effective reduction in infection rates, suggesting its importance in implementing effective intervention strategies. Lastly, Figure 4 focused on understanding the influence of the transmission probability on the tungiasis dynamics. Through this simulation, we visualized how changes in the transmission probability affect the overall dynamics of the infection. The results clearly illustrated that the transmission probability is a critical factor in shaping the infection spread within the community.

Our numerical analysis provided valuable insights into the control and dynamics of tungiasis. The fractional parameter, losing rate of protection, treatment rate, and transmission probability all emerged as significant factors that influence the infection's behavior. In the realm of tungiasis control and intervention, our discovered insights provide a foundational framework for the formulation of effective strategies aimed at managing and mitigating the spread of the infection within affected populations. It is, however, crucial to acknowledge that the values employed in this analysis were assumptions crafted for numerical expediency and may not faithfully replicate real-world scenarios. Consequently, for heightened precision in predictions and practical implications, an imperative next step involves the meticulous validation and refinement of our model using authentic real-world data. This process is indispensable for aligning our theoretical framework with the intricacies of tangible scenarios and enhancing the applicability of our findings in the context of tungiasis control and management.

7. Conclusion

In this study, we formulated an epidemic model for tungiasis disease in the framework of fractional Atangana-Baleanu derivative to conceptualize the transmission route of the infection. The biologically meaningful steady-states of the system are investigated through analytic methods. We determined the basic reproduction number of our fractional model, symbolized by $\mathcal{R}_0$ . The existence and uniqueness of the model's solution has been demonstrated. We have shown that the infection-free steady-states of the recommended model are locally asymptotically stable if $\mathcal{R}_0<1$ and unstable in other cases, while the endemic steady-state is locally asymptotically stable if $\mathcal{R}_0>1$ and unstable in other circumstances. A numerical scheme is presented to illustrate the solution pathways of the recommended system of the infection. The dynamical behaviour of the model is presented with the variation of different input parameters. We recommended the most critical factors of the system for the control and subsequent of tungiasis. It is acknowledged that delays play a fundamental role in capturing the temporal dynamics of systems [40]. Incorporating delays into mathematical models enhances their predictive power and allows for a more accurate representation of real-world phenomena [41]–[43]. In the future work, we will incorporate time delay in the transmission dynamics of the disease to comprehend the dynamics of the disease and provide more accurate information for prevention.

Use of AI tools declaration

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

References

[1]	Ackermann J (2019) Europe suffers from the sorry state of its banks. The Financial Times. Available from: https://www.ft.com/content/497d4d9e-0f3a-11e9-b2f2-f4c566a4fc5f
[2]	Altunbas Y, Carbo S, Gardener EPM, et al. (2007) Examining the Relationships between Capital, Risk and Efficiency in European Banking. Eur Financ Manage 13: 49-70. doi:10.1111/j.1468-036x.2006.00285.x doi: 10.1111/j.1468-036X.2006.00285.x
[3]	Anginer D, Demirguc-Kunt A (2014) Bank Capital and Systemic Stability. Policy Research Working Papers. doi:10.1596/1813-9450-6948
[4]	Ashraf B, Arshad S, Hu Y (2016) Capital Regulation and Bank Risk-Taking Behavior: Evidence from Pakistan. Int J Financ Stud 4: 16. doi:10.3390/ijfs4030016 doi: 10.3390/ijfs4030016
[5]	Banker RD, Chang H, Lee SY (2010). Differential impact of Korean banking system reforms on bank productivity. J Bank Financ 34: 1450-1460. doi:10.1016/j.jbankfin.2010.02.023 doi: 10.1016/j.jbankfin.2010.02.023
[6]	Barra C, Zotti R (2018) Bank performance, financial stability, and market concentration: Evidence from cooperative and non-cooperative banks. Ann Public Coop Econ 90: 103-139. doi:10.1111/apce.12217 doi: 10.1111/apce.12217
[7]	Bashir A, Hassan A (2017) Interrelationship Among Basel Capital Regulation, Risk, and Efficiency in Pakistani Commercial Banks. Bus Econ Rev 9: 165-186. doi:10.22547/ber/9.2.7 doi: 10.22547/BER/9.2.7
[8]	Ben Salah Mahdi I, Boujelbene Abbes M (2018) Relationship between capital, risk and liquidity: a comparative study between Islamic and conventional banks in MENA region. Res Int Bus Financ 45: 588-596. doi:10.1016/j.ribaf.2017.07.113 doi: 10.1016/j.ribaf.2017.07.113
[9]	Binham C, Noonan L (2015) Bad loans at Europe banks double that of the U.S. The Financial Times. Available from: https://www.ft.com/content/3ff8b5a0-92cf-11e5-94e6-c5413829caa5
[10]	Bitar M, Pukthuanthong K, Walker T (2018) The effect of capital ratios on the risk, efficiency and profitability of banks: Evidence from OECD countries. J Int Financ Markets Inst Money 53: 227-262. doi:10.1016/j.intfin.2017.12.002 doi: 10.1016/j.intfin.2017.12.002
[11]	Casu B, Girardone C (2009) Does Competition Lead to Efficiency? The Case of EU Commercial Banks. Ssrn Electron J. doi:10.2139/ssrn.1200362
[12]	Charnes A, Cooper WW, Rhodes E (1978) Measuring the efficiency of decision making units. Eur J Oper Res 2: 429-444. doi:10.1016/0377-2217(78)90138-8 doi: 10.1016/0377-2217(78)90138-8
[13]	Christopoulos AG, Dokas IG, Katsimardou S, et al. (2020) Assessing banking sectors' efficiency of financially troubled Eurozone countries. Res Int Bus Financ 52: 101121. doi:10.1016/j.ribaf.2019.101121 doi: 10.1016/j.ribaf.2019.101121
[14]	Dahl D, Shrieves RE (1990) The impact of regulation on bank equity infusions. J Bank Financ 14: 1209-1228. doi:10.1016/0378-4266(90)90010-y doi: 10.1016/0378-4266(90)90010-Y
[15]	Deelchand T, Padgett C (2009) The Relationship between Risk, Capital and Efficiency: Evidence from Japanese Cooperative Banks. Ssrn Electron J. doi:10.2139/ssrn.1525423
[16]	Ding D, Sickles RC (2019) Capital Regulation, Efficiency, and Risk Taking: A Spatial Panel Analysis of US Banks. Panel Data Econometrics, 405-466. doi:10.1016/b978-0-12-815859-3.00013-5 doi: 10.1016/B978-0-12-815859-3.00013-5
[17]	European Central Bank (2007) Financial Stability Review.
[18]	Federal Reserve Bank (2019) Return on average assets for all U.S. banks. Available from: https://fred.stlouisfed.org/series/USROA
[19]	Feng G, Wang C (2018) Why European banks are less profitable than U.S. banks: A decomposition approach. J Bank Financ 90: 1-16. doi:10.1016/j.jbankfin.2018.02.017
[20]	Fiordelisi F, Girardone C, Radic N (2011a) Price Competition, Efficiency and Riskiness in Investment Banking. Ssrn Electron J. doi:10.2139/ssrn.1786753
[21]	Fiordelisi F, Marques-Ibanez D, Molyneux P (2011b) Efficiency and Risk Taking in European Banking. Ssrn Electron J. doi:10.2139/ssrn.1512619
[22]	Ghosh S (2014) Risk, capital and financial crisis: Evidence for GCC banks. Borsa Istanbul Rev 14: 145-157. doi:10.1016/j.bir.2014.06.003 doi: 10.1016/j.bir.2014.06.003
[23]	Golin J, Delhaise P (2013) The Bank Credit Analysis Handbook: A Guide For Analysts, Bankers And Investors, Hoboken, NJ: John Wiley & Sons.
[24]	Goodhart C, Schoenmaker D (2016) The Global Investment Banks are now all Becoming American: Does that Matter for Europeans? J Financ Regul 2:163-181. doi: 10.1093/jfr/fjw012
[25]	Haq M, Heaney R (2012) Factors determining European bank risk. J Int Financ Markets Inst Money 22: 696-718. doi:10.1016/j.intfin.2012.04.003 doi: 10.1016/j.intfin.2012.04.003
[26]	Hou X, Wang Q, Zhang Q (2014) Market structure, risk taking, and the efficiency of Chinese commercial banks. Emerg Markets Rev 20: 75-88. doi:10.1016/j.ememar.2014.06.001 doi: 10.1016/j.ememar.2014.06.001
[27]	Jacques K, Nigro P (1997) Risk-based capital, portfolio risk, and bank capital: A simultaneous equations approach. J Econ Bus 49: 533-547. doi:10.1016/s0148-6195(97)00038-6 doi: 10.1016/S0148-6195(97)00038-6
[28]	Jenkins P (2015) What has delayed Europe's bank recovery? The Financial Times. Available from: https://www.ft.com/content/9e6f450c-86d8-11e5-90de-f44762bf9896
[29]	Kabir MN, Worthington AC (2017) The "competition-stability/fragility" nexus: A comparative analysis of Islamic and conventional banks. Int Rev Financ Anal 50: 111-128. doi:10.1016/j.irfa.2017.02.006 doi: 10.1016/j.irfa.2017.02.006
[30]	Kim D, Santomero AM (1988) Risk in Banking and Capital Regulation. J Financ 43: 1219. doi:10.2307/2328216 doi: 10.1111/j.1540-6261.1988.tb03966.x
[31]	Kok C, Móré C, Petrescu M (2016) Recent Trends In Euro Area Banks' Business Models And Implications For Banking Stability Sector. Financial Stability Review 1.
[32]	Kwan S, Eisenbeis RA (1997) Bank Risk, Capitalization and Inefficiency. SSRN Electron J. doi:10.2139/ssrn.1188
[33]	Lakhani K, Heid J, Templeman L (2019) How to fix European banking… and why it matters. Deutsche Bank.
[34]	Le T (2018) Bank Risk, Capitalisation and Technical Efficiency in the Vietnamese Banking System. Aust Account Bus Financ J 12: 41-61. doi:10.14453/aabfj.v12i3.4
[35]	Lin J, Chen S, Tsai J (2019) How does soft information about small business lending affect bank efficiency under capital regulation? Quant Financ Econ 3: 53-74. doi:10.3934/qfe.2019.1.53 doi: 10.3934/QFE.2019.1.53
[36]	Maity S, Ganguly D (2019) Is demonetization really impact efficiency of banking sector-An empirical study of banks in India. Asian J Multidimensional Res (AJMR) 8: 315. doi:10.5958/2278-4853.2019.00108.3 doi: 10.5958/2278-4853.2019.00108.3
[37]	Marjanović I, Stanković JJ, Popović Ž (2018) Efficiency Estimation of Commercial Banks Based on Financial Performance: Input Oriented DEA CRS/VRS Models. Econ Themes 56: 239-252. doi:10.2478/ethemes-2018-0014 doi: 10.2478/ethemes-2018-0014
[38]	McLannahan B, Arnold N (2017) US bank's profits leave European rivals in the shade. Financ Times. Available from: https://www.ft.com/content/d5b38c74-d9e0-11e6-944b-e7eb37a6aa8e
[39]	Miah MD, Sharmeen K (2015) Relationship between capital, risk and efficiency: A comparative study between Islamic and conventional banks in Bangladesh. Int J Islamic Middle East Financ Manage 8: 203-221. doi:10.1108/imefm-03-2014-0027 doi: 10.1108/IMEFM-03-2014-0027
[40]	Molyneux P, Reghezza A, Xie R (2019) Bank margins and profits in a world of negative rates. J Bank Financ 107: 105613. doi:10.1016/j.jbankfin.2019.105613 doi: 10.1016/j.jbankfin.2019.105613
[41]	Mosko A, Bozdo A (2016) Modeling the Relationship between Bank Efficiency, Capital and Risk in Albanian Banking System. Procedia Econ Financ 39: 319-327. doi:10.1016/s2212-5671(16)30330-6 doi: 10.1016/S2212-5671(16)30330-6
[42]	Moyo B (2018) An analysis of competition, efficiency and soundness in the South African banking sector. South Afr J Econ Manage Sci 21. doi:10.4102/sajems.v21i1.2291
[43]	Nguyen TP, Nghiem SH (2015) The interrelationships among default risk, capital ratio and efficiency. Managerial Financ 41: 507-525. doi:10.1108/mf-12-2013-0354 doi: 10.1108/MF-12-2013-0354
[44]	Peura S, Keppo J (2006) Optimal Bank Capital with Costly Recapitalization. J Bus 79: 2163-2201. doi:10.1086/503660 doi: 10.1086/503660
[45]	Reboul P, Florentin J, Robien E, et al. (2018) European corporate and investment banking. Outlook 2018. Eurogroup Conculting.
[46]	Repullo R (2002) Capital requirements, market power, and risk-taking in banking. J Financ Intermediation 13: 156-182. doi:10.1016/j.jfi.2003.08.005 doi: 10.1016/j.jfi.2003.08.005
[47]	Rime B (2001) Capital requirements and bank behaviour: Empirical evidence for Switzerland. J Bank Financ 25: 789-805. doi:10.1016/s0378-4266(00)00105-9 doi: 10.1016/S0378-4266(00)00105-9
[48]	Saeed M, Izzeldin M (2016) Examining the relationship between default risk and efficiency in Islamic and conventional banks. J Econ Behav Organ 132: 127-154. doi:10.1016/j.jebo.2014.02.014 doi: 10.1016/j.jebo.2014.02.014
[49]	Salas V, Saurina J (2003) Deregulation, market power and risk behaviour in Spanish banks. Eur Econ Rev 47: 1061-1075. doi:10.1016/s0014-2921(02)00230-1 doi: 10.1016/S0014-2921(02)00230-1
[50]	Schildbach J (2017) Where do European banks stand? Deutsche Bank Res.
[51]	Sealey CW, Lindley JT (1977) Inputs, Outputs, and a Theory of Production and Cost at Depository Financial Institutions. J Financ 32: 1251. doi:10.2307/2326527 doi: 10.1111/j.1540-6261.1977.tb03324.x
[52]	Shah AA, Wu D, Korotkov V (2019) Are Sustainable Banks Efficient and Productive? A Data Envelopment Analysis and the Malmquist Productivity Index Analysis. Sustainability 11: 2398. doi:10.3390/su11082398
[53]	Shim J (2013) Bank capital buffer and portfolio risk: The influence of business cycle and revenue diversification. J Bank Financ 37: 761-772. doi:10.1016/j.jbankfin.2012.10.002 doi: 10.1016/j.jbankfin.2012.10.002
[54]	Simm J, Besstremyannaya G (2016) CRAN-Package rDEA. Available from: https://CRAN.R-project.org/package=rDEA
[55]	Sufian F (2016) Determinants of Efficiency in the Malaysian Banking Sector: Evidence from Semi-parametric Data Envelopment Analysis Method. Stud Microecon 4: 151-172. doi:10.1177/2321022216636443 doi: 10.1177/2321022216636443
[56]	Tan Y, Floros C (2013) Risk, capital and efficiency in Chinese banking. J Int Financ Markets Inst Money 26: 378-393. doi:10.1016/j.intfin.2013.07.009 doi: 10.1016/j.intfin.2013.07.009
[57]	Teixeira JC, Silva FJ, Fernandes AV, et al. (2014) Banks' capital, regulation and the financial crisis. North Am J Econ Financ 28: 33-58. doi:10.1016/j.najef.2014.01.002 doi: 10.1016/j.najef.2014.01.002
[58]	Teixeira JC, Matos TF, Da Costa GL, et al. (2019) Investor protection, regulation and bank risk-taking behavior. North Am J Econ Financ, 101051. doi:10.1016/j.najef.2019.101051
[59]	Teixeira JC, Silva FJ, Costa FA, et al. (2019) Banks' profitability, institutions, and regulation in the context of the financial crisis. Int J Financ Econ doi:10.1002/ijfe.1753
[60]	Weigand RA (2016) The performance and risk of banks in the U.S., Europe and Japan post-financial crisis. Investment Manage Financ Innovations 13: 75-94. doi:10.21511/imfi.13(4).2016.07 doi: 10.21511/imfi.13(4).2016.07
[61]	Williams J (2004) Determining management behaviour in European banking. J Bank Financ 28: 2427-2460. doi:10.1016/j.jbankfin.2003.09.010 doi: 10.1016/j.jbankfin.2003.09.010
[62]	Zheng C, Moudud-Ul-Huq S, Rahman MM, et al. (2017) Does the ownership structure matter for banks' capital regulation and risk-taking behavior? Empirical evidence from a developing country. Res Int Bus Financ 42: 404-421. doi:10.1016/j.ribaf.2017.07.035 doi: 10.1016/j.ribaf.2017.07.035
[63]	Zhang J, Jiang C, Qu B, et al. (2013) Market concentration, risk-taking, and bank performance: Evidence from emerging economies. Int Rev Financ Anal 30: 149-157. doi:10.1016/j.irfa.2013.07.016 doi: 10.1016/j.irfa.2013.07.016

This article has been cited by:

Dongfei Gao, A deep learning approach for energy management systems in smart buildings towards a low-carbon economy, 2025, 20, 1748-1325, 1136, 10.1093/ijlct/ctaf063

Reader Comments

Your name:*

Email:*
© 2020 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)