From Theory to Code: Transforming Classical Root-Finding Methods into Efficient Python Implementations
Article Information
Abstract
This study conducts a comparative evaluation of seven numerical methods for finding roots of nonlinear equations: Bisection, Regula-Falsi, Fixed-Point Iteration, Newton-Raphson, Secant, Aitken's \(\Delta^2\), and Steffensen. The aim is to analyze the effectiveness of each method based on convergence speed, numerical accuracy, stability, and computational time efficiency. Algorithm implementation was carried out in the Python programming language using NumPy, SymPy, Pandas, and Matplotlib libraries. Test functions included polynomial, trigonometric, exponential, and mixed functions to represent diverse functional characteristics. The results indicate that Steffensen and Newton-Raphson achieved the fastest convergence in terms of iteration count, while Secant excelled in execution time efficiency. Bracketing methods such as Bisection and Regula-Falsi guaranteed convergence stability despite being slower. Fixed-Point Iteration was highly sensitive to the choice of iteration function, whereas Aitken's \(\Delta^2\) served effectively as an accelerator. The study concludes that method selection should be tailored to function characteristics and practical needs, such as derivative availability, computational complexity, and error tolerance. The implication is that this research provides an empirical guide for researchers and practitioners in selecting optimal numerical algorithms for scientific and engineering applications.
Graphical Abstract
Keywords
Data Availability Statement
Funding
Conflicts of Interest
Ethical Approval and Consent to Participate
References
- Chaudhary, R., Rawat, S., & Chauhan, P. (2024). Unraveling the roots: A comprehensive review of numerical methods for root finding. Journal of Emerging Technologies and Innovative Research, 11(1).
[Google Scholar] - Qureshi, S., Argyros, I. K., Soomro, A., Gdawiec, K., Shaikh, A. A., & Hincal, E. (2024). A new optimal root-finding iterative algorithm: Local and semilocal analysis with polynomiography. Numerical Algorithms, 95(4), 1715–1745.
[CrossRef] [Google Scholar] - Kelley, C. T. (2018). Numerical methods for nonlinear equations. Acta Numerica, 27, 207–287.
[CrossRef] [Google Scholar] - Gezerlis, A. (2023). Numerical methods in physics with Python (Vol. 1). Cambridge, UK: Cambridge University Press.
[CrossRef] [Google Scholar] - Zaguskin, V. L. (2014). Handbook of numerical methods for the solution of algebraic and transcendental equations. Elsevier.
[Google Scholar] - Cordero, A., Hueso, J. L., Martínez, E., & Torregrosa, J. R. (2012). Steffensen type methods for solving nonlinear equations. Journal of Computational and Applied Mathematics, 236(12), 3058–3064.
[CrossRef] [Google Scholar] - Saad, Y. (2025). Acceleration methods for fixed-point iterations. Acta Numerica, 34, 805-890.
[CrossRef] [Google Scholar] - Van Der Walt, S., Colbert, S. C., & Varoquaux, G. (2011). The NumPy array: A structure for efficient numerical computation. Computing in Science & Engineering, 13(2), 22–30.
[CrossRef] [Google Scholar] - Azure, I., Aloliga, G., & Doabil, L. (2019). Comparative study of numerical methods for solving non-linear equations using manual computation. Mathematics Letters, 5(4), 41-46.
[CrossRef] [Google Scholar] - Thakur, G., & Saini, J. K. (2021). Comparative study of iterative methods for solving non-linear equations. Journal of University of Shanghai for Science and Technology, 23(7), 858–866.
[CrossRef] [Google Scholar] - Ahmad, I., Ullah, M. A., & Uddin, J. (2024). Analysis and Enhancement of Newton Raphson Method. The Sciencetech, 5(2), 13-23.
[Google Scholar] - Bhardwaj, R. (2025). Numerical Simulation of Nonlinear Equations by Modified Bisection and Regula Falsi Method. Proceedings of the Pakistan Academy of Sciences: A. Physical and Computational Sciences, 62(1), 11-19.
[CrossRef] [Google Scholar] - Olanrewaju, A. F., Fadugba, S. E., & Adeyemi, O. P. (2024). Comparative study of numerical solutions to non-linear equations using regula-falsi and newton-raphson method. In 2024 International Conference on Science, Engineering and Business for Driving Sustainable Development Goals (SEB4SDG) (pp. 1–6). IEEE.
[CrossRef] [Google Scholar] - Hanief, M. (2021). Combined regula–falsi and newton–raphson method. In Recent Advances in Mathematical Methods and Economic Theory (pp. 123–132). Springer.
[CrossRef] [Google Scholar] - Fernández-Díaz, J. M., & Menéndez-Pérez, C. O. (2023). A common framework for modified regula falsi methods and new methods of this kind. Mathematics and Computers in Simulation, 205, 678–696.
[CrossRef] [Google Scholar] - Buldaev, A., & Kazmin, I. (2024). Fixed point methods for solving boundary value problem of the maximum principle. Journal of Mathematical Sciences, 279(6), 763–775.
[CrossRef] [Google Scholar] - Khuri, S. A., & Louhichi, I. (2021). A new fixed point iteration method for nonlinear third-order BVPs. International Journal of Computer Mathematics, 98(11), 2220–2232.
[CrossRef] [Google Scholar] - Thota, S., & Srivastav, V. K. (2018). Quadratically convergent algorithm for computing real root of non-linear transcendental equations. BMC Research Notes, 11(1), 909.
[CrossRef] [Google Scholar] - Argyros, I. K., Argyros, C., Ceballos, J., & González, D. (2022). Extended comparative study between Newton’s and Steffensen-like methods with applications. Mathematics, 10(16), 2851.
[CrossRef] [Google Scholar] - Dixit, N. D., & Mathur, P. K. (2021). Comparision of Numerical Accuracy of Bisection, Newton Raphson, Falsi-Position and Secant Methods. Advances in Mathematics: Scientific Journal, 10.
[Google Scholar] - Chen, C., & Liao, X. (2025). Solving non-linear equations by fixed point iteration method and its accelerating approach. Thermal Science, 29(3 Part A), 2031–2039.
[CrossRef] [Google Scholar] - Gulshan, G., Budak, H., Hussain, R., & Sadiq, A. (2023). Generalization of the bisection method and its applications in nonlinear equations. Advances in Continuous and Discrete Models, 2023(1), 18.
[CrossRef] [Google Scholar] - Mueen, H. A., & Shiker, M. A. K. (2024). Comparison Newton method with bisection method for solving nonlinear equations. In 2024 8th International Symposium on Multidisciplinary Studies and Innovative Technologies (ISMSIT) (pp. 1–4). IEEE.
[CrossRef] [Google Scholar] - Wang, S. (2025). A Comparative Study of Five Root-Finding Methods for Applying the Black-Scholes Model in Real Markets. Available at SSRN 5381527.
[CrossRef] [Google Scholar] - Cywiak, M., & Cywiak, D. (2021). Multi-Platform Graphics Programming with Kivy (pp. 173-190). Apress.
[CrossRef] [Google Scholar] - Steele, J. S., & Grimm, K. J. (2024). Using SymPy (Symbolic Python) for understanding structural equation modeling. Structural Equation Modeling: A Multidisciplinary Journal, 31(6), 1104–1115.
[CrossRef] [Google Scholar] - Ali, A., & Khusro, S. (2024). SA-MEAS: Sympy-based automated mathematical equations analysis and solver. SoftwareX, 25, 101596.
[CrossRef] [Google Scholar] - Harris, C. R., Millman, K. J., Van Der Walt, S. J., Gommers, R., Virtanen, P., Cournapeau, D., ... & Oliphant, T. E. (2020). Array programming with NumPy. nature, 585(7825), 357-362.
[CrossRef] [Google Scholar] - Gupta, P., & Bagchi, A. (2024). Introduction to NumPy. In Essentials of Python for Artificial Intelligence and Machine Learning (pp. 127-159). Cham: Springer Nature Switzerland.
[CrossRef] [Google Scholar] - Unpingco, J. (2021). Numpy. In Python programming for data analysis (pp. 103-126). Cham: Springer International Publishing.
[CrossRef] [Google Scholar] - Fuhrer, C., Solem, J. E., & Verdier, O. (2021). Scientific Computing with Python: High-performance scientific computing with NumPy, SciPy, and pandas. Packt Publishing Ltd.
[Google Scholar] - Gupta, P., & Bagchi, A. (2024). Introduction to pandas. In Essentials of python for artificial intelligence and machine learning (pp. 161-196). Cham: Springer Nature Switzerland.
[CrossRef] [Google Scholar] - Roy, S. S., Chouhan, A., & Khera, N. (2025). Python Fast Track: A Complete Guide to Rapidly Mastering and Applying Python Programming. Elsevier.
[CrossRef] [Google Scholar] - Mattingly, W. J. B. (2023). Introduction to Pandas. In Introduction to Python for Humanists (pp. 77–84). Chapman and Hall/CRC.
[CrossRef] [Google Scholar] - Lavanya, A., Gaurav, L., Sindhuja, S., Seam, H., Joydeep, M., Uppalapati, V., ... & SD, V. S. (2023). Assessing the performance of Python data visualization libraries: a review. Int. J. Comput. Eng. Res. Trends, 10(1), 28-39.
[CrossRef] [Google Scholar] - Haslwanter, T. (2016). An introduction to statistics with python. With applications in the life sciences. Switzerland: Springer International Publishing.
[Google Scholar] - Stančin, I., & Jović, A. (2019, May). An overview and comparison of free Python libraries for data mining and big data analysis. In 2019 42nd International convention on information and communication technology, electronics and microelectronics (MIPRO) (pp. 977-982). IEEE.
[CrossRef] [Google Scholar] - Gan, G. (2025). The Matplotlib Library. In Data Clustering with Python (pp. 66–74). Chapman and Hall/CRC.
[CrossRef] [Google Scholar] - Amat, S., Ezquerro, J. A., & Hernández-Verón, M. A. (2016). On a Steffensen-like method for solving nonlinear equations. Calcolo, 53(2), 171–188.
[CrossRef] [Google Scholar] - Zhou, Y., Ding, F., Alsaedi, A., & Hayat, T. (2021). Aitken-based acceleration estimation algorithms for a nonlinear model with exponential terms by using the decomposition. International Journal of Control, Automation and Systems, 19(11), 3720–3730.
[CrossRef] [Google Scholar] - Fika, P., & Mitrouli, M. (2017). Aitken’s method for estimating bilinear forms arising in applications. Calcolo, 54(1), 455–470.
[CrossRef] [Google Scholar] - Zhao, Y., Fu, S., Zhang, L., & Huang, H. (2025). Aitken optimizer: An efficient optimization algorithm based on the Aitken acceleration method. The Journal of Supercomputing, 81(1), 264.
[CrossRef] [Google Scholar]
Cited By (1)
-
Richa Sharma, Virendra Singh Chouhan. Computational performance of classical and hybrid root-finding methods for real-world nonlinear models.
Discover Computing, 2026 , 29 (1).
[CrossRef]
Cite This Article
TY - JOUR AU - Romadani, Moh. Nafis Husen PY - 2025 DA - 2025/12/14 TI - From Theory to Code: Transforming Classical Root-Finding Methods into Efficient Python Implementations JO - ICCK Journal of Applied Mathematics T2 - ICCK Journal of Applied Mathematics JF - ICCK Journal of Applied Mathematics VL - 1 IS - 3 SP - 154 EP - 189 DO - 10.62762/JAM.2025.840767 UR - https://www.icck.org/article/abs/JAM.2025.840767 KW - numerical methods KW - nonlinear equations KW - root-finding algorithms KW - python implementation KW - computational mathematics AB - This study conducts a comparative evaluation of seven numerical methods for finding roots of nonlinear equations: Bisection, Regula-Falsi, Fixed-Point Iteration, Newton-Raphson, Secant, Aitken's \(\Delta^2\), and Steffensen. The aim is to analyze the effectiveness of each method based on convergence speed, numerical accuracy, stability, and computational time efficiency. Algorithm implementation was carried out in the Python programming language using NumPy, SymPy, Pandas, and Matplotlib libraries. Test functions included polynomial, trigonometric, exponential, and mixed functions to represent diverse functional characteristics. The results indicate that Steffensen and Newton-Raphson achieved the fastest convergence in terms of iteration count, while Secant excelled in execution time efficiency. Bracketing methods such as Bisection and Regula-Falsi guaranteed convergence stability despite being slower. Fixed-Point Iteration was highly sensitive to the choice of iteration function, whereas Aitken's \(\Delta^2\) served effectively as an accelerator. The study concludes that method selection should be tailored to function characteristics and practical needs, such as derivative availability, computational complexity, and error tolerance. The implication is that this research provides an empirical guide for researchers and practitioners in selecting optimal numerical algorithms for scientific and engineering applications. SN - 3068-5656 PB - Institute of Central Computation and Knowledge LA - English ER -
@article{Romadani2025From,
author = {Moh. Nafis Husen Romadani},
title = {From Theory to Code: Transforming Classical Root-Finding Methods into Efficient Python Implementations},
journal = {ICCK Journal of Applied Mathematics},
year = {2025},
volume = {1},
number = {3},
pages = {154-189},
doi = {10.62762/JAM.2025.840767},
url = {https://www.icck.org/article/abs/JAM.2025.840767},
abstract = {This study conducts a comparative evaluation of seven numerical methods for finding roots of nonlinear equations: Bisection, Regula-Falsi, Fixed-Point Iteration, Newton-Raphson, Secant, Aitken's \(\Delta^2\), and Steffensen. The aim is to analyze the effectiveness of each method based on convergence speed, numerical accuracy, stability, and computational time efficiency. Algorithm implementation was carried out in the Python programming language using NumPy, SymPy, Pandas, and Matplotlib libraries. Test functions included polynomial, trigonometric, exponential, and mixed functions to represent diverse functional characteristics. The results indicate that Steffensen and Newton-Raphson achieved the fastest convergence in terms of iteration count, while Secant excelled in execution time efficiency. Bracketing methods such as Bisection and Regula-Falsi guaranteed convergence stability despite being slower. Fixed-Point Iteration was highly sensitive to the choice of iteration function, whereas Aitken's \(\Delta^2\) served effectively as an accelerator. The study concludes that method selection should be tailored to function characteristics and practical needs, such as derivative availability, computational complexity, and error tolerance. The implication is that this research provides an empirical guide for researchers and practitioners in selecting optimal numerical algorithms for scientific and engineering applications.},
keywords = {numerical methods, nonlinear equations, root-finding algorithms, python implementation, computational mathematics},
issn = {3068-5656},
publisher = {Institute of Central Computation and Knowledge}
}
Publisher's Note
ICCK stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and Permissions
Copyright © 2025 by the Author(s). Published by Institute of Central Computation and Knowledge. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.