https://tnsmb.org/journal/index.php/ijmam/issue/feed International Journal of Mathematical Analysis and Modelling 2024-11-24T07:18:56+00:00 Daniel Okuonghae (Ph.D.), Editor-in-Chief ijmam.nsmb@gmail.com Open Journal Systems <p>The International Journal of Mathematical Analysis and Modelling (IJMAM), formerly known as the Journal of the Nigerian Society for Mathematical Biology, is the official journal of the Nigerian Society for Mathematical Biology, designed to disseminates original research findings in all areas of pure and applied mathematics.</p> <p>The aim of the journal is to cause to be published, original research that shows the connection between Pure and Applied Mathematicians, Statisticians, Operations Researchers, Physicists, Engineers, Computer Scientists and specialists in any field of human endeavour such as in Agriculture, Medicine, Pharmacy, Life Sciences, Physical Sciences, Engineering, Social Sciences, Humanities, Arts, etc, at the crossroad of theoretical, experimental and field studies with the Mathematical Sciences, in order to provide solutions to problems that cut across these fields through exchange of ideas and collaborations.</p> <p>Submitted papers should provide rigorous mathematical analysis of any real-world problem being discussed; at the same time, the results obtained should be useful and easily understandable by non-mathematicians. The journal will accept different kinds of articles, ranging from review articles to special invitations. All submitted manuscript will be assessed by the Editor-in-Chief before being sent out for peer-review.</p> <p>Areas covered in the IJMAM include, Mathematical Modelling, Numerical Analysis, Fluid Mechanics, Real and Complex Analysis, Functional Analysis, Operator Theory, Statistics, Algebra, Computational Mathematics and Differential Equations. Articles that show the connections between different areas of Mathematics when solving real-world problems are especially welcome in the IJMAM. </p> <p><a href="https://tnsmb.org/journal/index.php/ijmam/about/editorialTeam">Click to see our Editorial Board members</a></p> https://tnsmb.org/journal/index.php/ijmam/article/view/169 An effective new one-step block integrator for the numerical solution of system of first order differential equation arising from epidemiology models 2024-11-24T06:28:42+00:00 M.I. Modebei gmarc.ify@gmail.com O.O. Olaiya oolaiya.o@gmail.com O.E. Faniyi ofaniyi@aul.edu.ng O.F. Anorue onyinyechianorue@gmail.com <p>This paper proposes an effective new one-step block integrator for the solution of system of first order ordinary differential equation that arises from models in epidemiology in life science. Showing the performance of this method, SIR, SEIR and more compartmental models are solved numerically and thus comparing with the classical Runge-Kutta to reveal its accuracy. The numerical results obtained clearly shows that the proposed method is more efficient and robust for obtaining the numerical solution of mathematical models in epidemiology.</p> 2024-11-24T00:00:00+00:00 Copyright (c) 2024 M.I. Modebei, O.O. Olaiya, O.E. Faniyi, O.F. Anorue https://tnsmb.org/journal/index.php/ijmam/article/view/153 Numerical solution approximation of nonlinear functional integral by Homotopy perturbation method 2024-09-24T10:01:45+00:00 C.E. Chika∗ charles.chika@unn.edu.ng <p>A numerical method for solving nonlinear functional integral is presented. Numerical examples are presented to show the validity of the numerical algorithm developed using Homotopy Perturbation Method (HPM). This method is simple to implement and gives good approximation at first iteration.</p> 2024-09-25T00:00:00+00:00 Copyright (c) 2024 C.E. Chika∗ https://tnsmb.org/journal/index.php/ijmam/article/view/167 Fractional order model of dynamical behavior and qualitative analysis of Anthrax with infected vector and saturation 2024-11-13T10:31:37+00:00 A.C. Loyinmi loyinmiac@tasued.edu.ng A.L. Ijaola‡ loyinmiac@tasued.edu.ng <p>We present a fractional order model to capture the transmission dynamics of anthrax infection with a nonlinear force of infection and its long-term impacts. An analysis with the standard Caputo–Fabrizio approach established that the model is well posed, and for the investigation on the stability of the model, an anthrax-free equilibrium as connected to the basic reproduction number R0 &lt; 1 and other important conditions that consolidate epidemiological feasibility of the model were established. Using a normalized forward sensitivity index, we conducted the sensitivity analysis of important variables, the contact rates, and the recruitment rate of the vectors to examine the effects of their variation on the dynamics of anthrax. These allow us to identify the most sensitive parameters that healthcare professionals need to pay attention to. The results from simulations also demonstrated that the presence of saturation instantly causes the system to approach an anthrax-free equilibrium (DFE), and our findings from the qualitative analysis compel us to recommend maximum hygiene practices for humans and domestic animals, early intervention, and the treatment of anthrax with a vaccine (as well as all other medical interventions) on the infected and periodic evaluation of these practices in households and the community at large to achieve a disease-free community.</p> 2024-11-13T00:00:00+00:00 Copyright (c) 2024 A.C. Loyinmi, A.L. Ijaola‡ https://tnsmb.org/journal/index.php/ijmam/article/view/151 The influence of aquatic mosquitoes on Lymphatic Filariasis transmission dynamics: mathematical model approach 2024-09-23T18:01:52+00:00 A. Abokwara∗† agathamicheal@gmail.com C.E. Madubueze‡ agathamicheal@gmail.com <p>The disease burden of lymphatic filariasis (LF) is a concern worldwide as it is one of the leading<br>causes of permanent disability. In this study, a mathematical model of lymphatic filariasis is formulated with interest in the aquatic mosquito population (egg, larva and pupa) on the metamorphosis<br>of the mosquito life cycle. The existence of the disease-free and the endemic equilibrium states and<br>the basic reproduction number are established. In addition, the sensitivity indices indicate that the<br>biting rate of mosquitoes, the recruitment rate of the aquatic mosquito population, the probability<br>of transmitting infection during contact between humans and mosquitoes and the mortality rate<br>of mosquitoes are the sensitive parameters in the disease dynamics. Through numerical simulations, it is observed that the sensitive parameters have an effect in reducing the burden of lymphatic<br>filariasis in the population. Thus, public health education and the usage of larvicides should be<br>encouraged as these would cause a drastic reduction in the maturation rate of mosquitoes, especially<br>in communities where stagnant waters are endemic.</p> 2024-09-25T00:00:00+00:00 Copyright (c) 2024 A. Abokwara∗†, C.E. Madubueze‡ https://tnsmb.org/journal/index.php/ijmam/article/view/159 Detection and correction of heteroscedasticity and its effect on modelling of Nigerian economic data 2024-09-25T11:34:42+00:00 E.O. Idowu*† godsbabe2002@gmail.com E.M. Ikegwu‡ ademola.fadiji@yabatech.edu.ng A.A. Fadiji†§ ademola.fadiji@yabatech.edu.ng M.U. Evro†** merit.evro@gmail.com <p>The study centred on detection and correction of heteroscedasticity and its effect on modelling of economic data. Data were collected from CBN Statistical bulletin on five economic variables namely; Gross Domestic Product (GDP), Inflation Rate, Exchange Rate (Exch_rate), Balance of Payment and Government Debt from 1987 - 2017. The data was analysed using the Ordinary Least Square method and variance stability techniques was used to correct heteroscedasticity using log transformation of the dependent variable while Breusch-Pagan, White and HarveyGodfrey tests were used to detect the presence of heteroscedasticity in the variables at 5% level of significance with the aid of E-views 9. The regression model fitted to the data set was GDP = 37242.63 - 256.7745 Inflr – 40.196 Exchr + 252.364 BoP + 3.147 GovD and it showed that only inflation rate and Government debt were discovered to be significant predictors (p &lt; 0.05) while and heteroscedasticity was detected. Result showed that after log transforming the variables, heteroscedasticity was eliminated as shown by Breusch-Pagan, White and HarveyGodfrey test respectively. It is concluded therefore that transforming economic data helps correct and eliminate detected heteroscedasticity in modelled economic variables.</p> 2024-09-25T00:00:00+00:00 Copyright (c) 2024 E.O. Idowu*†, E.M. Ikegwu‡, A.A. Fadiji†§, M.U. Evro†** https://tnsmb.org/journal/index.php/ijmam/article/view/172 Effect of thermal radiation and chemical reaction on unsteady MHD plane – Poiseuille flow of fourth-grade fluid in horizontal parallel plates 2024-11-24T07:16:28+00:00 K.M. Joseph kpop.moses@kasu.edu.ng E.A. Andi eaandi@nda.edu.ng E.J. Adoyi jeadoyi@nda.edu.ng C.O. Akusu cakusu@nda.edu.ng C. Onwubuoua onwubuoyacletus@novenauniversity.edu.ng <p>A fourth-grade fluid is an important subclass of differential type that is capable of describing shear thinning and shear thickening effects (examples are ketchup, blood, paint, cream, nail polish, etc.). This sort of model is used to explain the flow behaviour of non-Newtonian fluids which are considered vital and applicable in many industrial production processes such as in the drilling of oil and gas wells, polymer extrusion from dye, glass fibre, paper production and draining of plastics films, etc. There is a dearth of knowledge in this area, therefore, this study investigated the thermal radiation and chemical reaction effects on unsteady magnetohydrodynamics plane – poiseuille flow of fourth-grade fluid in horizontal parallel plates channel. The equations that governed the flow were the momentum equation; which is a coupled nonlinear differential equation, the energy equation and concentration equation. The partial differential equations were solved using He – Laplace method which is a combination of Homotopy perturbation method and Laplace transformation method. The findings of the study revealed that: (i) Velocity and temperature fields rise due to the increment of thermal radiation parameter. (ii) For upsurging data of chemical reaction, velocity and concentration fields diminish. (iii) Velocity profile goes up when third and fourth-grade parameters get to raise; Velocity and skin friction fields decline due to the increment of magnetic parameter. (iv) Increasing Prandtl number tend to diminish the velocity and temperature profiles. (v) Strong values of Schmidt number decrease the boundary layer of the Sherwood number field. The results of this work are applicable to industrial processes such polymer extrusion of dye, draining of plastic films etc.</p> 2024-11-24T00:00:00+00:00 Copyright (c) 2024 K.M. Joseph, E.A. Andi, E.J. Adoyi, C.O. Akusu, C. Onwubuoua https://tnsmb.org/journal/index.php/ijmam/article/view/157 Mathematical model for transmission dynamics of Malaria disease: impact of invasive alien plants 2024-09-25T11:01:29+00:00 C.J. Alhassan∗ charity.alhassan@edouniversity.edu.ng D. Okuonghae daniel.okuonghae@uniben.edu K.O. Achema‡ achema.kenneth@uam.edu.ng <p>A mathematical model to study the impact of invasive alien plants on the dynamics of malaria transmission and its analyses is studied in this work. The resulted model equations are divided into homogeneous and non-homogeneous equations. The homogeneous equations are solved to determine its disease free equilibrium (DFE) and their stabilities. A basic reproduction number was determined from the DFE. It was found that when R0 &lt; 1, the disease will die out, when R0 = 1, the model undergoes a backward bifurcation, R0 = 0, the model undergoes forward bifurcation and whenever R0 &gt; 1, the disease will persist in the population. At R0 &gt; 1, the global analysis of the model was carried out and it was found to be globally and asymptotically<br>stable (GAS). A sensitivity analysis of the model parameters was also investigated to determine the parameters that are sensitive for malaria transmission. A travelling wave equations and solutions were also provided for possible understanding of the behaviour of mosquito’s mobility in human environment.</p> 2024-09-25T00:00:00+00:00 Copyright (c) 2024 C.J. Alhassan∗, D. Okuonghae, K.O. Achema‡ https://tnsmb.org/journal/index.php/ijmam/article/view/170 Mathematical model showing preventive measure as a key to reducing Coronavirus disease (COVID-19) spread 2024-11-24T06:49:07+00:00 C. Okoye chukwudi.okoye@unn.edu.ng M.B. Okofu mary.okofu@unn.edu.ng N.K. Nubila nkoyo.nubila@unn.edu.ng C.L. Ejikeme chioma.ejikeme@unn.edu.ng D.F. Agbebaku dennis.agbebaku@unn.edu.ng <p>A set of non-linear Mathematical equations of COVID-19 model, presented in a flow diagram was proposed to study and ascertain the impact of preventive measures (self isolation and social distancing) to curtailing the spread of the disease. To this effect, we computed the threshold quantity Reff (Effective Reproduction Number), carried out the stability analysis on the equilibrium points obtained using Jacobian technique (Eigenvalue approach), Castilo-Chavez method and a suitable Lyapunov function. We discovered that the Disease Free equilibrium point failed the long term test whereas the Endemic Equilibrium was globally asymptotically stable, thus implying that the COVID-19 spread could be curtailed or curbed in the presence of preventive measures (self isolation and social distancing) as advocated.</p> 2024-11-24T00:00:00+00:00 Copyright (c) 2024 C. Okoye, M.B. Okofu, N.K. Nubila, C.L. Ejikeme , D.F. Agbebaku https://tnsmb.org/journal/index.php/ijmam/article/view/154 A theoretical investigation of the impact of a 2-dose vaccination strategy on Monkeypox transmission dynamics at population level 2024-09-24T10:07:35+00:00 I.I. Ako∗ ignatius.ako@uniben.edu <p>In order to qualitatively evaluate the influence of a two-dose monkeypox (MPV) vaccination strategy when decreased human-to-human infection is accounted for, a new deterministic model for<br>MPV is created and implemented. The backward bifurcation phenomena is demonstrated to occur in the model as a result of the parameter that describes how quickly susceptible individuals receive their first dose of the MPV vaccination. Additionally, a distinct cutoff point for the rate at which vulnerable individuals receive the first dose of the MPV vaccination was discovered An example<br>of the model demonstrated that, in the absence of the rate at which susceptible individuals receive the first dose of the MPV vaccine, the disease-free equilibrium was globally asymptotically stable. The endemic equilibrium point of the distinguished case was also shown to globally asymptotically stable.</p> 2024-09-25T00:00:00+00:00 Copyright (c) 2024 I.I. Ako∗ https://tnsmb.org/journal/index.php/ijmam/article/view/168 An infection-age-structured modelling for Monkey-Pox disease dynamics incorporating control measures 2024-11-13T11:17:51+00:00 N. O. Lasisi *† nurudeenlasisi2009@yahoo.com O.J. Ejiwole ‡ enitioluwafe2012@gmail.com E. Azuaba†§ emmanuel.azuaba@binghamuni.edu.ng A. A. Ibrahim†** ibrahimabdukareema@fedponam.edu.ng.com <p>Monkey-pox is known as pathogens affecting livestock animals and humans and belongs to the orthopox virus. The pathogen causes lymph nodes to swell and increasing transmission risk associated with factors involving introduction of virus to the oral mucosa. In this paper, we developed an Age-Structured model for Monkey-pox disease in a population with vital dynamics, incorporating standard incidence rate and vaccination. We showed the existence and uniqueness of the solution of the model. We obtained the Disease-Free Equilibrium state and shown the effective reproduction number of the model. We proved the conditions for Local and Global Stability of the Disease-Free Equilibrium (DFE) State and we found that the disease free equilibrium state is locally asymptotically stable if 1 )0( G (RE &lt; 1) and Globally Asymptotically Stable (GAS) in Ω if RE ≤ 1 while unstable if RE ≥ 1.</p> 2024-11-13T00:00:00+00:00 Copyright (c) 2024 N. O. Lasisi *†, O.J. Ejiwole ‡, E. Azuaba†§, A. A. Ibrahim†** https://tnsmb.org/journal/index.php/ijmam/article/view/152 Approximate common solution of variational inequality and fixed point problems for family of nonexpansive semigroup 2024-09-24T09:52:59+00:00 A.H. Ugwu∗ albertugwu5@gmail.com I.G. Ezugorie† ikegodezugorie@esut.edu.ng E.C. Agbo‡ agk.ejike@gmail.com E.U. Ofoedu§ euofoedu@yahoo.com <p>In this paper, we prove a path convergence theorem and introduce a new iterative scheme for approximation of common element of fixed point sets of family of nonexpansive semigroup and the set of solutions of variational inequality problem for an α-inverse strongly accretive mapping in strictly convex q-uniformly smooth real Banach space. As an application, solution is given to the problem of finding a common element of fixed point sets of family of nonexpansive semigroup and the set of zeros of an α-inverse strongly accretive operator. The result obtained augment and generalize several existing results.&nbsp;</p> 2024-09-25T00:00:00+00:00 Copyright (c) 2024 A.H. Ugwu∗, I.G. Ezugorie†, E.C. Agbo‡, E.U. Ofoedu§ https://tnsmb.org/journal/index.php/ijmam/article/view/160 Options pricing using combination of Brownian motion with a random time component 2024-09-25T11:45:32+00:00 C.A. Onyegbuchulem*† ochialuka@yahoo.com B. O. Onyegbuchulem‡ ochialuka@yahoo.com U. Alwell† ochialuka@yahoo.com <p>This paper provides an efficient option pricing technique for European call options. The technique is a process that combines Geometric Brownian motion with random time component. The additional parameter to the Geometric Brownian motion is the volatility of the time change. This additional parameter provides control over the skewness and kurtosis of the return distribution. Partial integro-differential equation (PIDE) for pricing options under the proposed process was derived. Its values were priced numerically and the pricing performance of the new model is compared on option data sourced from Bloomberg on 28th April, 2023 with the Geometric Brownian motion. The result revealed that the additional parameter was able to address the pricing biases of the Geometric Brownian motion. Recommendation was made that more work should be done in the advancement of quick and effective procedures for calculating PIDE’s having in view to calibrating stochastic processes to a surface of option prices.</p> 2024-09-25T00:00:00+00:00 Copyright (c) 2024 C.A. Onyegbuchulem*†, B. O. Onyegbuchulem‡, U. Alwell† https://tnsmb.org/journal/index.php/ijmam/article/view/150 Utility maximization and portfolio management with administrative fee and dividend 2024-09-23T17:47:58+00:00 E.E. Akpanibah∗† edemae@fuotuoke.edu.ng L.E. George† edemae@fuotuoke.edu.ng W.V. German† edemae@fuotuoke.edu.ng <p>The study of utility maximization cannot be over emphasized in the study of portfolio management<br>as it plays crucial role in the determination of investor’s strategy over any investment. In this<br>research work, we investigate some utility functions and their applications to portfolio optimization<br>and risk management. Also, a portfolio consisting of one risk free asset and a risky asset which<br>follows the geometric Brownian motion was considered. We took into consideration transaction<br>cost, dividend and tax on invested funds. By applying Ito’s lemma and maximum principle, we<br>obtained an optimization problem which comprise of a nonlinear PDE (Hamilton Jacobi Bellman<br>equation) and a control problem (optimal investment strategy). These problems are functions of the value function which is dependent on the utility of the investor at the expiration of such investment. Furthermore, the exponential and logarithm utility were used in solving for the optimal investment strategies and some numerical results also presented. It was observed for both cases that the optimal investment strategies (OIS) were directly proportional to appreciation rate, dividend and inversely proportional to the administrative charges, instantaneous volatility, tax and risk averse coefficient. Also, the OIS under logarithmic utility does not depend on the risk averse coefficient and tax while the optimal investment strategy under exponential utility does not depend on wealth.&nbsp;</p> 2024-09-25T00:00:00+00:00 Copyright (c) 2024 E.E. Akpanibah∗† , L.E. George†, W.V. German† https://tnsmb.org/journal/index.php/ijmam/article/view/158 On the nilpotents of the semigroup of partial contractions of a finite chain 2024-09-25T11:17:02+00:00 M.M. Zubairu∗† mmzubairu.mth@buk.edu.ng B Ali bali@nda.edu.ng <p>Let [n] = {1, 2, . . . , n} denote a finite chain. Consider the semigroups Pn and OIn, consisting of all partial, and all injective and order-preserving transformations on [n], respectively. Furthermore, let CPn = {α ∈ Pn : |xα − yα| ≤ |x − y| ∀ x, y ∈ Dom α} and OCIn = {α ∈ OIn : |xα − yα| ≤ |x − y| ∀ x, y ∈ Dom α}. It follows that CPn and OCIn are subsemigroups of Pn and OIn, respectively. In this paper, we characterize the nilpotent elements in CPn, compute the number of nilpotents of height (n − 1) in CPn, and determine the number of nilpotents in OCIn of a specific type</p> 2024-09-25T00:00:00+00:00 Copyright (c) 2024 M.M. Zubairu∗†, B Ali https://tnsmb.org/journal/index.php/ijmam/article/view/171 On the vector Lyapunov functions and eventual stability of nonlinear impulsive differential equations 2024-11-24T07:02:30+00:00 J.E. Ante jackson.ante@topfaith.edu.ng A.B. Inyang ab.inyang@topfaith.edu.ng E.J. Oduobuk ej.oduobuk@topfaith.edu.ng U.P. Akai up.akai@topfaith.edu.ng <p>In this paper, the eventual stability of nonlinear impulsive differential equations with fixed moments of impulse is examined using the vector Lyapunov functions which is generalized by a class of piecewise continuous functions. The novelty in the use of the vector form of the Lypunov function lies in the fact that the "restrictions" encountered via the use of scalar Lyapunov function is safely handled especially for large scale dynamical systems, since it involves splitting the Lyapunov functions into components so that each of the components can easily describe the behaviour of the solution state Together with comparison results, sufficient conditions for eventual stability are presented. Results obtained are extensions and improvements on existing results.</p> 2024-11-24T00:00:00+00:00 Copyright (c) 2024 J.E. Ante, A.B. Inyang, E.J. Oduobuk, U.P. Akai https://tnsmb.org/journal/index.php/ijmam/article/view/155 Complex Ebola Virus Disease fractional order model in the Caputo sense with Laplace Adomian Decomposition Method model solution 2024-09-24T10:27:59+00:00 Q.O. Ahman ∗ williamsatokolo@gmail.com R.O. Aja†† williamsatokolo@gmail.com W. Atokolo ‡§ williamsatokolo@gmail.com D. Omale¶ williamsatokolo@gmail.com N.E. Didigwu williamsatokolo@gmail.com C.S. Ugwu∗∗ williamsatokolo@gmail.com A.L. Ozioko†† williamsatokolo@gmail.com S.O. Joseph‡‡ williamsatokolo@gmail.com V.I. Ezaegu∗ williamsatokolo@gmail.com R.V. Paul† williamsatokolo@gmail.com G.C.E. Mbah‡ williamsatokolo@gmail.com <p>Ebola Virus Disease is a life-threatening disease that is transmitted between humans and animals. The aim of this study is to develop a fractional order Ebola Virus Disease model for the dynamics of Ebola considering the control measures;quarantine,isolation,treatment, vaccine and protection (condom use) with a view to study the effects of these control measures on the spread of Ebola in a given population of humans and animals in the Caputo Sense. The fractional order Ebola Virus Disease model in the Caputo sense to control the spread of the disease with eighteen compartments is considered in this paper.This model considered vaccination, condom use, quarantine, use of treatment drugs and isolation as control measures combined together. We analyzed the model where the validity of the model was proven by establishing a region of invariance and the positivity of solutions. We also established the model disease-free equilibrium (DFE)and endemic equilibrium (EE). The reproduction number (a threshold value) (B) was obtained using the next generation matrix. The result of the study produced an expression for B, where if the reproduction number B, is less than 1 then the model DFE point is stable, so the Ebola virus disease dies out. If the reproduction number B is greater than 1, then the model DFE point is unstable, so the Ebola Virus disease persists in the populations.We also analysed the local stability of the DFE point and found out it will be stable when the reproduction number B &lt; 1. The novel application of the Laplace–Adomian decomposition method to the model obtained a result of the complex fractional model in an infinite series form which converges further to its exact value. Comparing the solutions of the fractional Ebola virus disease model with that of the classical case using simulation plots we found out that the case of fractional order has more degree of freedom in such a way that can be varied as the fraction (α) could be varied to get different results. We proved convergence of the Laplace Adomian Decomposition method using the Cauchy-Kovalevskaya theorem for differential<br>equations with analytic vector fields and then obtained a new result on the convergence rate of the Laplace Adomian Decomposition Method.&nbsp;</p> 2024-09-25T00:00:00+00:00 Copyright (c) 2024 Q.O. Ahman ∗, R.O. Aja††, W. Atokolo ‡§, D. Omale¶, N.E. Didigwu, C.S. Ugwu∗∗, A.L. Ozioko††, S.O. Joseph‡‡, V.I. Ezaegu∗, R.V. Paul†, G.C.E. Mbah‡