Forward Bifurcation with Hysteresis Phenomena from Atherosclerosis Mathematical Model


  • Dipo Aldila Department of Mathematics, Universitas Indonesia, Depok 16424, Indonesia
  • Arthana Islamilova Department of Mathematics, Universitas Indonesia, Depok 16424, Indonesia
  • Sarbaz H.A. Khosnaw Department of Mathematics, University of Raparin, Ranya, Kurdistan Region, Iraq
  • Bevina D. Handari Department of Mathematics, Universitas Indonesia, Depok 16424, Indonesia
  • Hengki Tasman Department of Mathematics, Universitas Indonesia, Depok 16424, Indonesia



atherosclerosis, forward bifurcation, backward bifurcation, hysteresis, optimal control


Atherosclerosis is a non-communicable disease (NCDs) which appears when the blood vessels in the human body become thick and stiff. The symptoms range from chest pain, sudden numbness in the arms or legs, temporary loss of vision in one eye, or even kidney failure, which may lead to death. Treatment in cases with severe symptoms requires surgery, in which the number of doctors or hospitals is limited in some countries, especially countries with low health levels. This article aims to propose a mathematical model to understand the impact of limited hospital resources on the success of the control program of atherosclerosis spreads. The model was constructed based on a deterministic model, where the hospitalization rate is defined as a time-dependent saturated function concerning the number of infected individuals. The existence and stability of all possible equilibrium points were shown analytically and numerically, along with the basic reproduction number. Our analysis indicates that our model may exhibit various types of bifurcation phenomena, such as forward bifurcation, backward bifurcation, or a forward bifurcation with hysteresis depending on the value of hospitalization saturation parameter and the infection rate for treated infected individuals. These phenomenon triggers a complex and tricky control program of atherosclerosis. A forward bifurcation with hysteresis auses a possible condition of having more than one stable endemic equilibrium when the basic reproduction number is larger than one, but close to one. The more significant value of hospitalization saturation rate or the infection rate for treated infected individuals increases the possibility of the stable endemic equilibrium point even though the disease-free equilibrium is stable. Furthermore, the Pontryagin Maximum Principle was used to characterize the optimal control problem for our model. Based on the results of our analysis, we conclude that atherosclerosis control interventions should prioritize prevention efforts over endemic reduction scenarios to avoid high intervention costs. In addition, the government also needs to pay great attention to the availability of hospital services for this disease to avoid the dynamic complexity of the spread of atherosclerosis in the field.

Author Biography

Dipo Aldila, Department of Mathematics, Universitas Indonesia, Depok 16424, Indonesia


Ross, R., The pathogenesis of atherosclerosis- an update, New England journal of medicine, 314(8), pp. 488-500, 1986.

Libby, P., Ridker, P. M., and Maseri, A., Inflammation and atherosclerosis, Circulation, 105(9), pp. 1135-1143, 2002.

Falk, E., Pathogenesis of atherosclerosis, Journal of the American College of cardiology, 47(8S), pp. C7-C12, 2006.

Haas, D. C., Davidson, K. W., Schwartz, D. J., Rieckmann, N., Roman, M. J., Pickering, T. G., ... and Schwartz, J. E., Depressive symptoms are independently predictive of carotid atherosclerosis, The American journal of cardiology, 95(4), pp. 547-550, 2005.

Vinkers, D. J., Stek, M. L., Van der Mast, R. C., De Craen, A. J. M., Le Cessie, S., Jolles, J., ... and Gussekloo, J., Generalized atherosclerosis, cognitive decline, and depressive symptoms in old age, Neurology, 65(1), pp. 107-112, 2005.

Chong, P. H., and Bachenheimer, B. S., Current, new and future treatments in dyslipidaemia and atherosclerosis, Drugs, 60(1), pp. 55-93, 2000.

Khatib, N.E., Kafi, O., Sequeira, A., Simakov, S., Vassilevski Yu, Volpert, V., Mathematical modelling of atherosclerosis, Math. Model. Nat. Phenom., 14 (2019), p. 603, 2019.

Calvez, V., Ebde, A., Meunier, N., Raoult, A., Mathematical and numerical modeling of the atherosclerotic plaque formation, ESAIM Proceedings, 28, pp.1?12, 2009.

Silva, T. Sequeira, A., Santos, R.F., Tiago, J., Mathematical Modeling of Atherosclerotic Plaque Formation Coupled with a Non-Newtonian Model of Blood Flow, Conference Papers in Science, vol. 2013, Article ID 405914.

Islamilova, A., Aldila, D., Gayatri, W., and Tasman, H., Modelling the spread of atherosclerosis considering relapse and linear treatment, Journal of Physics: Conference Series, 1722(1), p. 012039, 2021.

Allen, L.N. and Feigl, A.B., Reframing non-communicable diseases as socially transmitted conditions, The Lancet Global Health, 5(7), pp. e644-e646, 2017.

Kompas, S., Aldila, D., and Handari, B.D., Modelling the Spread of Diabetes Transmission through Social Contact, AIP Conf. Proceeding, 1722(1), p. 012039, 2020.

Delavani, Aldila, D., Handari, B.D., Effect of healthy life campaigns on controlling obesity transmission: A Mathematical Study, Journal of Physics: Conference Series, 1747(1), p. 012003, 2021.

Al-Tuwairqi, S.M., Matbouli, R.T., Modeling dynamics of fast food and obesity for evaluating the peer pressure effect and workout impact, Advances in Difference Equations, 2021(1), p. 59, 2021.

Jodar,L.,Santonja,F.J.,Gonzalez-Parra,G., Modeling dynamics of infant obesity intheregion of Valencia, Spain, Comput. Math. Appl., 56(3), pp. 679?689, 2008.

Diekmann, O., Heesterbeek, J. A. P., and Roberts, M. G., The construction of next-generation matrices for compartmental epidemic models, Journal of Royal Society Interface, 7(47), pp. 873-885, 2010.

Aldila, D., Ndii, M.Z., Samiadji, B.M., Optimal control on COVID-19 eradication program in Indonesia under the effect of community awareness, MBE, 17(6), pp. 6355-6389, 2020.

Simorangkir, G., Aldila, D., Rizka, A., Tasman, H., Nugraha, E.S., Mathematical model of tuberculosis considering observed treatment and vaccination interventions, Journal of Interdisciplinary Mathematics, 24(6), pp. 1717?1737, 2021.

Balya, M.A., et al, Investigating the Impact of Social Awareness and Rapid Test on A COVID-19 Transmission Model, Commun. Biomath. Sci., 4(1), pp. 46-64, 2021.

Maimunah, Aldila, D., Mathematical model for HIV spreads control program with ART treatment, Journal of Physics: Conference Series, 974(1), p. 012035, 2018.

Driessche, P.V., Watmough, J., Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math Biosci., 180(1?2), pp. 29?48, 2002.

Castillo-Chavez, C., Song, B., Dynamical models of tuberculosis and their applications, Math. Biosci. Eng., 1(2), pp. 361?404, 2004.

Gulbudak, H., Martcheva, M., Forward Hysteresis and Backward Bifurcation Caused by Culling in an Avian Influenza Model, Mathematical Biosciences, 246(1), pp. 202-212, 2013.

Pontryagin, L.S., Boltyanski, V.G., Gamkrelidze, R.V., Mischenke, E.V., The Mathematical Theory of Optimal Process, Wiley, New Jersey, 1962.

Aldila, D., Optimal control for dengue eradication program under the media awareness effect, International Journal of Nonlinear Sciences and Numerical Simulation, 2021.

Aldila, D., Angelina, M., Optimal control problem and backward bifurcation on malaria transmission with vector bias, Heliyon, 7(4), p. e06824, 2021.