Numerical Error Comparison of Euler, Heun, and Fourth-Order Runge–Kutta Methods for a SICR Tuberculosis Transmission Model in South Batipuh, West Sumatera
DOI:
https://doi.org/10.64570/jamm.v2i1.63Keywords:
Euler method, Heun method, mathematical modeling, numerical error, Runge–Kutta method, SICR model, tuberculosisAbstract
Tuberculosis (TB) remains a major infectious-disease problem and requires quantitative approaches for understanding transmission and treatment dynamics. This study formulates a SICR (Susceptible–Infected–Care–Recovered) compartmental model for TB transmission in South Batipuh District, Tanah Datar Regency, West Sumatra, Indonesia, and compares the numerical accuracy of the Euler, Heun, and fourth-order Runge–Kutta (RK4) methods. Secondary data for 2020–2024 were obtained from population records and TB case records from the local health service. The model consists of susceptible, infected, care, and recovered compartments and is represented by a nonlinear ordinary differential equation system. Model parameters were estimated by curve fitting, and all three numerical schemes were implemented using the same initial conditions and time-step setting. The observed data show that the infected population increased from 69 cases in 2020 to 166 cases in 2023 before decreasing to 132 cases in 2024, while the care and recovered compartments also fluctuated. Numerical simulations reproduce the qualitative dynamics of declining susceptible individuals, an initial increase in infected individuals, a subsequent increase in the care compartment, and growth of the recovered compartment. The average relative errors were 0.10843 for Euler, 0.08342 for Heun, and 0.08317 for RK4. Therefore, RK4 produced the smallest average error and was the most accurate method among the three methods for the studied SICR system. The findings support the use of higher-order numerical integration when accurate approximation of nonlinear TB transmission dynamics is required.
References
Aldila, D., Ramadhan, D. A., Chukwu, C. W., Handari, B. D., Shahzad, M., & Kamalia, P. Z. (2024). On the role of early case detection and treatment failure in controlling tuberculosis transmission: A mathematical modeling study. Communications in Biomathematical Sciences, 7(1), 61–86. https://doi.org/10.5614/cbms.2024.7.1.4
Burden, R. L., & Faires, J. D. (2011). Numerical analysis (9th ed.). Brooks/Cole, Cengage Learning.
Chapra, S. C., & Canale, R. P. (2015). Numerical methods for engineers (7th ed.). McGraw-Hill Education.
Citra, D. Y., Glorya, E. B., Helmina, K., Ambarwati, L., & Wiraningsih, E. D. (2024). Perbandingan nilai galat pada model penyakit tuberkulosis menggunakan metode Euler dan Heun di DKI Jakarta. JMT (Jurnal Matematika dan Terapan), 6(2). https://doi.org/10.21009/jmt.6.2.2
Fatmawati, Purwati, U. D., Utoyo, M. I., Alfiniyah, C., & Prihartini, Y. (2020). The dynamics of tuberculosis transmission with optimal control analysis in Indonesia. Communications in Mathematical Biology and Neuroscience, 2020:25. https://doi.org/10.28919/cmbn/4605
Jannah, M. (2026). Perbandingan nilai galat pada model SICR penyebaran penyakit TBC dengan metode Euler, Heun dan Runge-Kutta orde empat: Studi kasus data TBC Kecamatan Batipuh Selatan [Undergraduate thesis, Universitas Dharma Andalas].
Kanwal, S. (2022). Analysis of the epidemic biological model of tuberculosis (TB) via numerical schemes. Complexity, 2022, 5147951. https://doi.org/10.1155/2022/5147951
Kasbawati, Tungga, R. A., Jaya, A. K., & Kalondeng, A. (2021). A mathematical study of tuberculosis infections using a deterministic model in comparison with continuous Markov chain model. Communications in Mathematical Biology and Neuroscience, 2021, 25. https://doi.org/10.28919/cmbn/5180
Muniroh, M. A. (2025). Simulation modeling incomplete treatment impact on tuberculosis transmission. Jurnal Matematika, Statistika dan Komputasi, 21(2), 456–466. https://doi.org/10.20956/j.v21i2.36825
Permatasari, A. H., & Utomo, R. H. S. (2023). Mathematical analysis of a tuberculosis model with two different stages of infection. Journal of Fundamental Mathematics and Applications, 6(1), 80–89. https://doi.org/10.14710/jfma.v5i2.15723
Simorangkir, G., Aldila, D., Rizka, A., Tasman, H., & Nugraha, E. S. (2021). Mathematical model of tuberculosis considering observed treatment and vaccination interventions. Journal of Interdisciplinary Mathematics, 24(6), 1717–1737. https://doi.org/10.1080/09720502.2021.1958515
Suddin, S., & Bano, E. N. (2021). Mathematical modeling and simulation to control the spread of multidrug-resistant tuberculosis. International Journal of Computer Science and Applied Mathematics, 7(1), 1–7.
Ullah, I., Ahmad, S., Al-Mdallal, Q., Khan, Z. A., Khan, H., & Khan, A. (2020). Stability analysis of a dynamical model of tuberculosis with incomplete treatment. Advances in Difference Equations, 2020, 499. https://doi.org/10.1186/s13662-020-02950-0
World Health Organization. (2024). Consolidated guidance on tuberculosis data generation and use: Module 1: Tuberculosis surveillance. https://www.who.int/publications/i/item/9789240075290
World Health Organization. (2024). Global tuberculosis report 2024. https://www.who.int/publications/i/item/9789240101531
