In this article, a numerical method integrated with statistical data simulation technique is introduced to solve a nonlinear system of ordinary differential equations with multiple random variable coefficients. The utilization of Monte Carlo simulation with central divided difference formula of finite difference (FD) method is repeated n times to simulate values of the variable coefficients as random sampling instead being limited as real values with respect to time. The mean of the n final solutions via this integrated technique, named in short as mean Monte Carlo finite difference (MMCFD) method, represents the final solution of the system. This method is proposed for the first time to calculate the numerical solution obtained fo

The aims of this paper is investigating the spread of AIDS both within-host, through the contact between healthy cells with free virus inside the body, and between-host, through sexual contact among individuals and external sources of infectious. The outbreak of AIDS is described by a mathematical model consisting of two stages. The first stage describes the within-host spread of AIDS and is represented by the first three equations. While the second stage describes the between-host spread of AIDS and represented by the last four equations. The existence, uniqueness and boundedness of the solution of the model are discussed and all possible equilibrium points are determined. The local asymptotic stability (LAS) of the model is studied, wh

In this paper, a harvested prey-predator model involving infectious disease in prey is considered. The existence, uniqueness and boundedness of the solution are discussed. The stability analysis of all possible equilibrium points are carried out. The persistence conditions of the system are established. The behavior of the system is simulated and bifurcation diagrams are obtained for different parameters. The results show that the existence of disease and harvesting can give rise to multiple attractors, including chaos, with variations in critical parameters.

In this paper, the dynamical behavior of a three-dimensional fractional-order prey-predator model is investigated with Holling type III functional response and constant rate harvesting. It is assumed that the middle predator species consumes only the prey species, and the top predator species consumes only the middle predator species. We also prove the boundedness, the non-negativity, the uniqueness, and the existence of the solutions of the proposed model. Then, all possible equilibria are determined, and the dynamical behaviors of the proposed model around the equilibrium points are investigated. Finally, numerical simulations results are presented to confirm the theoretical results and to give a better understanding of the dynami

Coronavirus disease 2019 (COVID-19) is a flu-like infection caused by a novel virus known as Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2). After the widespread around the world, it was announced by the World Health Organization (WHO) as a global pandemic. The symptoms of COVID-19 may arise within 2 weeks and the severity ranged from mild with signs of respiratory infection to severe cases of organ failure and even death. Management of COVID-19 patients includes supportive treatment and pharmacological medications expected to be effective with no definitive cure of the disease. The aims of this study are highlighting the management protocol and supportive therapy especially vitamin D and manifesting the clinical symptoms b

This article examines and proposes a dietary chain model with a prey shelter and alternative food sources. It is anticipated that mid-predators' availability is positively correlated with the number of refuges. The solution's existence and exclusivity are examined. It is established that the solution is bounded. It is explored whether all potential equilibrium points exist and are locally stable. The Lyapunov approach is used to investigate the equilibrium points' worldwide stability. Utilizing a Sotomayor theorem application, local bifurcation is studied. Numerical simulation is used to better comprehend the dynamics of the model and define the control set of parameters.

This paper deals with the blow-up properties of positive solutions to a parabolic system of two heat equations, defined on a ball in  associated with coupled Neumann boundary conditions of exponential type. The upper bounds of blow-up rate estimates are derived. Moreover, it is proved that the blow-up in this problem can only occur on the boundary.