A discrete fractional-order mathematical model for in vitro fertilization dynamics

Iqbal M. Batiha; Nidal Anakira; Irianto Irianto; Abed Al-Rahman M. Malkawi; Tala Sasa; Shaher Momani; Iqbal M. Batiha; Nidal Anakira; Irianto Irianto; Abed Al-Rahman M. Malkawi; Tala Sasa; Shaher Momani

doi:10.3934/math.2026622

AIMS Mathematics

2026, Volume 11, Issue 5: 15120-15142. doi: 10.3934/math.2026622

Previous Article Next Article

Research article

A discrete fractional-order mathematical model for in vitro fertilization dynamics

1.
Department of Mathematics, Al Zaytoonah University of Jordan, Amman 11733, Jordan
2.
Nonlinear Dynamics Research Center (NDRC), Ajman University, Ajman 346, United Arab Emirates
3.
Mathematics Education Program, Faculty of Education and Arts, Sohar University, Sohar 311, Oman
4.
Department General Education, Faculty of Resilience, Rabdan Academy, Abu Dhabi, United Arab Emirates
5.
Department of Mathematics, Faculty of Arts and Science, Amman Arab University, Amman 11953, Jordan
6.
Department of Mathematics, Faculty of Science, Private Applied Science University, Amman, Jordan
7.
Department of Mathematics, University of Jordan, Amman 19392, Jordan

Received: 27 March 2026 Revised: 13 May 2026 Accepted: 18 May 2026 Published: 29 May 2026
MSC : 39A30, 26A33, 92D25

Infertility affects millions of couples worldwide, and in vitro fertilization (IVF) is the foremost assisted reproductive technology. However, IVF success rates vary considerably due to the memory-dependent nature of hormonal regulation, embryo development, and repeated treatment cycles. Clinical data are typically collected at discrete monthly intervals, yet few mathematical models simultaneously capture both the discrete time structure and the memory effects inherent in IVF. This work develops and rigorously analyses a discrete fractional-order (FO) IVF model using the dual Caputo nabla fractional difference (CNFD) operator, which naturally incorporates long-range memory while aligning with the cycle-based format of medical records. The model stratifies the IVF process into six compartments: infertile couples, patients under treatment, high-, medium-, and low-quality embryos, and positive pregnancy outcomes. Existence, uniqueness, and non-negativity of solutions are proved. Equilibrium analysis yields a unique positive endemic equilibrium, and sufficient conditions for global asymptotic stability (GAS) and Mittag-Leffler stability (MLS) are established via a novel Volterra-type discrete Lyapunov function (LF). Numerical simulations, performed with representative parameter values, confirm the theoretical results and show that lower FOs introduce stronger memory effects and slower convergence, thereby reproducing realistic, protracted IVF dynamics. The theoretical findings are not tied to a specific experimental dataset and thus provide a general framework that could assist clinicians in optimizing treatment decisions and personalizing patient management.
- discrete fractional calculus,
- Caputo nabla fractional difference,
- in vitro fertilization model,
- global asymptotic stability,
- Mittag-Leffler stability,
- memory effects
Citation: Iqbal M. Batiha, Nidal Anakira, Irianto Irianto, Abed Al-Rahman M. Malkawi, Tala Sasa, Shaher Momani. A discrete fractional-order mathematical model for in vitro fertilization dynamics[J]. AIMS Mathematics, 2026, 11(5): 15120-15142. doi: 10.3934/math.2026622

Related Papers:

Abstract

Infertility affects millions of couples worldwide, and in vitro fertilization (IVF) is the foremost assisted reproductive technology. However, IVF success rates vary considerably due to the memory-dependent nature of hormonal regulation, embryo development, and repeated treatment cycles. Clinical data are typically collected at discrete monthly intervals, yet few mathematical models simultaneously capture both the discrete time structure and the memory effects inherent in IVF. This work develops and rigorously analyses a discrete fractional-order (FO) IVF model using the dual Caputo nabla fractional difference (CNFD) operator, which naturally incorporates long-range memory while aligning with the cycle-based format of medical records. The model stratifies the IVF process into six compartments: infertile couples, patients under treatment, high-, medium-, and low-quality embryos, and positive pregnancy outcomes. Existence, uniqueness, and non-negativity of solutions are proved. Equilibrium analysis yields a unique positive endemic equilibrium, and sufficient conditions for global asymptotic stability (GAS) and Mittag-Leffler stability (MLS) are established via a novel Volterra-type discrete Lyapunov function (LF). Numerical simulations, performed with representative parameter values, confirm the theoretical results and show that lower FOs introduce stronger memory effects and slower convergence, thereby reproducing realistic, protracted IVF dynamics. The theoretical findings are not tied to a specific experimental dataset and thus provide a general framework that could assist clinicians in optimizing treatment decisions and personalizing patient management.

References

[1]	B. O. Ibitoye, O. Akadiri, F. O. Ibitoye, O. A. Fasasi, WHO fact sheet on infertility (definition of infertility), published in Global Reproductive Health (2021), Global Reprod. Health, 8 (2023), e60. https://doi.org/10.1097/GRH.0000000000000060 doi: 10.1097/GRH.0000000000000060
[2]	H. Fertilisation, E. Authority, The current status of IVF: are we putting the needs of the individual first, eClinicalMedicine, 65 (2023), 102343. https://doi.org/10.1016/j.eclinm.2023.102343 doi: 10.1016/j.eclinm.2023.102343
[3]	W. M. Hirst, A. Vail, D. R. Brison, S. A. Roberts, Prognostic factors influencing fresh and frozen IVF outcomes: an analysis of the UK national database, Reprod. BioMed. Online, 22 (2011), 437–448. https://doi.org/10.1016/j.rbmo.2010.12.008 doi: 10.1016/j.rbmo.2010.12.008
[4]	Y. Bentov, J. Schenker, IVF and pregnancy outcomes: the triumphs, challenges, and unanswered questions, J. Ovarian Res., 18 (2025), 228. https://doi.org/10.1186/s13048-025-01692-5 doi: 10.1186/s13048-025-01692-5
[5]	A. Al-Nana, I. M. Batiha, Modeling adaptive memory in neural dynamics: a dual nabla fractional-order FitzHugh–Nagumo framework, Ricerche Mat., 2026. https://doi.org/10.1007/s11587-026-01067-y doi: 10.1007/s11587-026-01067-y
[6]	A. K. Alsharidi, M. Qousini, Dynamics and stability of soliton structures for the generalized nonlinear fractional $(3+1)$-dimensional wave model in computational physics, Fractal Fract., 9 (2025), 806. https://doi.org/10.3390/fractalfract9120806 doi: 10.3390/fractalfract9120806
[7]	N. Hafizi, P. Tulay, Mathematical modeling in reproduction and infertility, In: I. Ozsahin, D. U. Ozsahin, Applied machine learning and multi-criteria decision-making in healthcare, Bentham Science Publishers, (2021), 214–234. https://doi.org/10.2174/9781681088716121010014
[8]	B. Yeolekar, N. Shukla, J. Shukla, M. Yeolekar, S. Patil, P. K. Santra, Optimal control model for IVF treatment in women, Int. J. Biomath., 2024, 2450101. https://doi.org/10.1142/S1793524524501018 doi: 10.1142/S1793524524501018
[9]	M. Ghannam, B. Kaymakamzade, M. Farman, K. S. Nisar, F. Aljuaydi, A. Sambas, Investigating the success of IVF with fractional-order model: analysis and case study, Fractals, 34 (2026), 2540255. https://doi.org/10.1142/S0218348X25402558 doi: 10.1142/S0218348X25402558
[10]	A. N. Anber, Z. Dahmani, The Laplace decomposition method for solving nonlinear conformable fractional evolution equations, Int. J. Open Probl. Comput. Math., 17 (2024), 67.
[11]	I. H. Jebril, Linear and nonlinear control for complete synchronization of fractional-order discrete reaction-diffusion systems, Int. J. Anal. Appl., 24 (2026), 26. https://doi.org/10.28924/2291-8639-24-2026-26 doi: 10.28924/2291-8639-24-2026-26
[12]	P. R. Ashok, S. A. Sayed, A comprehensive review of fractional differential equations and their role in modeling complex dynamical systems, Int. J. Eng. Sci. Humanit., 14 (2024), 71–80.
[13]	F. M. Atıcı, M. Atıcı, N. Nguyen, T. Zhoroev, G. Koch, A study on discrete and discrete fractional pharmacokinetics-pharmacodynamics models for tumor growth and anti-cancer effects, Comput. Math. Biophys., 7 (2019), 10–24. https://doi.org/10.1515/cmb-2019-0002 doi: 10.1515/cmb-2019-0002
[14]	M. Liu, F. Stella, A. Hommersom, P. J. Lucas, L. Boer, E. Bischoff, A comparison between discrete and continuous time Bayesian networks in learning from clinical time series data with irregularity, Artif. Intell. Med., 95 (2019), 104–117. https://doi.org/10.1016/j.artmed.2018.10.002 doi: 10.1016/j.artmed.2018.10.002
[15]	I. H. Jebril, I. Batiha, A stable and convergent implicit finite difference scheme for variable-order time-fractional convection–diffusion equations, Int. J. Open Probl. Comput. Math., 19 (2026), 73–94.
[16]	M. Mohammadi, M. Kavousi, T. Madani, P. Amini, A. Ghaheri, Joint modeling of in vitro fertilization outcomes among a population of Iranian infertile couples: a historical cohort study, Int. J. Fertil. Steril., 17 (2023), 306–311. https://doi.org/10.22074/IJFS.2023.562653.1374 doi: 10.22074/IJFS.2023.562653.1374
[17]	F. H. Alawad, A. A. Alharbi, K. I. Mayoof, H. M. R. Hawsawi, N. N. Alharthi, R. A. Ahmed, Comparative retrospective analysis of clinical and hormonal profiles in PCOS patients with and without infertility, J. Adv. Trends Med. Res., 1 (2024), 513–518. https://doi.org/10.4103/ATMR.ATMR_21_24 doi: 10.4103/ATMR.ATMR_21_24
[18]	L. Sterling, J. Liu, N. Okun, A. Sakhuja, S. Sierra, E. Greenblatt, Pregnancy outcomes in women with polycystic ovary syndrome undergoing in vitro fertilization, Fertil. Steril., 105 (2016), 791–797. https://doi.org/10.1016/j.fertnstert.2015.11.019 doi: 10.1016/j.fertnstert.2015.11.019
[19]	A. Templeton, J. K. Morris, W. Parslow, Factors that affect outcome of in-vitro fertilisation treatment, Lancet, 348 (1996), 1402–1406. https://doi.org/10.1016/S0140-6736(96)05291-9 doi: 10.1016/S0140-6736(96)05291-9
[20]	L. G. Cooney, M. D. Sammel, I. Lee, M. A. Clapp, M. Goldsammler, E. Scott, et al., The details matter: personalized prediction of live birth after in vitro fertilization in women with polycystic ovary syndrome, Fertil. Steril., 121 (2024), 1010–1019. https://doi.org/10.1016/j.fertnstert.2024.01.033 doi: 10.1016/j.fertnstert.2024.01.033
[21]	L. Sterling, J. Liu, N. Okun, A. Sakhuja, S. Sierra, E. Greenblatt, Pregnancy outcomes in women with polycystic ovary syndrome undergoing in vitro fertilization, Fertil. Steril., 105 (2016), 791–797. https://doi.org/10.1016/j.fertnstert.2015.11.019 doi: 10.1016/j.fertnstert.2015.11.019
[22]	M. Batool, M. Farman, A. Ahmad, K. S. Nisar, Mathematical study of polycystic ovarian syndrome disease including medication treatment mechanism for infertility in women, AIMS Public Health, 11 (2023), 19–25. https://doi.org/10.3934/publichealth.2024002 doi: 10.3934/publichealth.2024002
[23]	Y. Wei, Y. Wei, Y. Chen, Y. Wang, Mittag–Leffler stability of nabla discrete fractional-order dynamic systems, Nonlinear Dyn., 101 (2020), 407–417. https://doi.org/10.1007/s11071-020-05776-3 doi: 10.1007/s11071-020-05776-3
[24]	T. Abdeljawad, On delta and nabla Caputo fractional differences and dual identities, Discrete Dyn. Nat. Soc., 2013 (2013), 406910. https://doi.org/10.1155/2013/406910 doi: 10.1155/2013/406910
[25]	X. Wu, X. Yang, Q. Song, X. Chen, Stability analysis on nabla discrete distributed-order dynamical system, Fractal Fract., 6 (2022), 429. https://doi.org/10.3390/fractalfract6080429 doi: 10.3390/fractalfract6080429
[26]	Q. Huang, Z. Tu, Discrete-time Mittag–Leffler state estimation for fractional-order quaternion memristive neural networks, Sci. Rep., 15 (2025), 38082. https://doi.org/10.1038/s41598-025-21865-y doi: 10.1038/s41598-025-21865-y
[27]	Y. Wei, Y. Wei, Y. Chen, Y. Wang, Mittag–Leffler stability of nabla discrete fractional-order dynamic systems, Nonlinear Dyn., 101 (2020), 407–417. https://doi.org/10.1007/s11071-020-05776-3 doi: 10.1007/s11071-020-05776-3
[28]	A. Al-Nana, I. M. Batiha, Modeling adaptive memory in neural dynamics: a dual nabla fractional-order FitzHugh–Nagumo framework, Ricerche Mat., 2026. https://doi.org/10.1007/s11587-026-01067-y doi: 10.1007/s11587-026-01067-y
[29]	A. Ghaheri, A. Rasekhi, R. O. Samani, E. Hajizadeh, Modeling in vitro fertilization data considering multiple outcomes observed among Iranian infertile women, Int. J. Fertil. Steril., 12 (2018), 27–30. https://doi.org/10.22074/ijfs.2018.5187 doi: 10.22074/ijfs.2018.5187
[30]	I. S. Jamadar, K. Kumar, S. A. Khan, A. Khan, M. N. Akhtar, E. A. Bakar, Quantum pressure and memory effects in cancer modeling: a fractional calculus neural network approach, Results Eng., 27 (2025), 106080. https://doi.org/10.1016/j.rineng.2025.106080 doi: 10.1016/j.rineng.2025.106080
[31]	J. Xu, D. O'Regan, Existence and uniqueness of solutions for a first-order discrete fractional boundary value problem, RACSAM, 112 (2018), 1005–1016. https://doi.org/s13398-017-0406-7
[32]	Z. Lv, C. Wu, D. O'Regan, J. Xu, Solvability of a boundary value problem involving fractional difference equations, Axioms, 12 (2023), 650. https://doi.org/10.3390/axioms12070650 doi: 10.3390/axioms12070650
[33]	J. Xu, C. S. Goodrich, Y. Cui, Positive solutions for a system of first-order discrete fractional boundary value problems with semipositone nonlinearities, RACSAM, 113 (2019), 1343–1358. https://doi.org/10.1007/s13398-018-0551-7 doi: 10.1007/s13398-018-0551-7
[34]	W. Cheng, J. Xu, Y. Cui, Q. Ge, Positive solutions for a class of fractional difference systems with coupled boundary conditions, Adv. Differ. Equ., 2019 (2019), 249. https://doi.org/10.1186/s13662-019-2184-3 doi: 10.1186/s13662-019-2184-3
[35]	N. Gopal, J. M. Jonnalagadda, Multiple positive solutions of discrete fractional boundary value problems, Kragujevac J. Math., 50 (2026), 25–38. https://doi.org/10.46793/KgJMat2601.025G doi: 10.46793/KgJMat2601.025G
[36]	I. Bendib, I. Batiha, A. Ouannas, N. Anakira, O. Ogilat, T. Sasa, On general fractional-order discrete-time reaction-diffusion systems: finite time stability and simulations, Int. J. Rob. Control Syst., 5 (2025), 2995–3015. https://doi.org/10.31763/ijrcs.v5i6.1979 doi: 10.31763/ijrcs.v5i6.1979
[37]	Y. H. Liang, K. J. Wang, The modified variational principles of the fractal Rosenau–Burgers equation, Fractals, 2026. https://doi.org/10.1142/S0218348X26500532 doi: 10.1142/S0218348X26500532
[38]	Y. Zhang, C. Zhang, K. J. Wang, Qualitative analysis of the fractal complex Hirota-dynamical model for optical fiber communication: variational principle, Hamiltonian and diverse soliton solutions, Fractals, 34 (2026), 2650022. https://doi.org/10.1142/S0218348X26500222 doi: 10.1142/S0218348X26500222
[39]	K. J. Wang, An effective computational approach to the local fractional low-pass electrical transmission lines model, Alexandria Eng. J., 110 (2025), 629–635. https://doi.org/10.1016/j.aej.2024.07.021 doi: 10.1016/j.aej.2024.07.021

Reader Comments

Your name:*

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