Numerical Modeling of Oil Palm Plantation Area Growth in Indonesia Using Verhulst Model
DOI:
https://doi.org/10.32665/james.v9i2.6337Keywords:
Backward Differentiation Formula, Milne Method, Oil Palm Plantation, Runge-Kutta Method, Verhulst ModelAbstract
The expansion of oil palm plantation area in Indonesia exhibits nonlinear growth behavior influenced by environmental, spatial, and policy constraints. This study aims to model the growth of oil palm plantation area using the Verhulst (logistic) growth model and to obtain numerical solutions of the resulting nonlinear differential equation. A quantitative approach is employed, in which the fourth-order Runge–Kutta (RK4) method is used to generate the required initial values, followed by two predictor–corrector schemes, namely the Adams–Bashforth–Backward Differentiation Formula (AB4–BDF4) and the Milne–Backward Differentiation Formula (Milne–BDF4), which are compared to approximate the long-term growth dynamics. The analysis utilizes annual oil palm plantation area data from 2017 to 2024 obtained from Statistics Indonesia (BPS). Numerical results indicate that all methods produce stable solutions, with relative errors remaining bounded and non-divergent throughout the simulation, on the order of 10-6. On average, the AB4–BDF4 scheme yields a relative error of 3.2477330 × 10-6, compared to 5.9206896 × 10-6 for the Milne–BDF4 scheme, corresponding to an accuracy improvement of approximately 1.8 times. The estimated intrinsic growth rate is r≈0.15, and the estimated carrying capacity is 20 million hectares, representing the maximum sustainable oil palm plantation area under current growth trends.
References
Abanum, G. C. (2020). Comparison Of Numerical Methods In Predicting The Growth Of Normal Agricultural Assets. Journal of Mathematical Sciences & Computational Mathematics, 1(4), 502–523. https://doi.org/10.15864/jmscm.1408
Adebisi, A. F., Aremu, S., Ogunniran, M. O., Tijani, K., & Ajileye, A. M. (2025). Enhanced Milne-Simpson’s Methods For Autonomous And Singular Differential Equations. Mathematics and Computational Sciences, 6(2), 125–140. https://doi.org/10.30511/mcs.2025.2035400.1205
Anggreini, D., Afifah, D. S. N., Hayuhantika, D., Hastari, R. C., & Puspasari, R. (2019). Population Projection Using The Implementation Of Differential Equation Of Logistic Models. ICBLP 2019: Proceedings of the 1st International Conference on Business, Law And Pedagogy, ICBLP 2019, 13-15 February 2019, Sidoarjo, Indonesia, 326. https://doi.org/10.4108/eai.13-2-2019.2286525
Anie, N. I. (2016). On The Accuracy Of The Backward Difference Formula Of Order Four And Runge-Kutta Of Order Four In Solving A First Order Non-Linear Ordinary Differential Equation. Kampala International University.
Audu, K. J., Garba, J., Tiamiyu, A. T., & Thomas, B. A. (2022). Application of Backward Differentiation Formula on Fourth-Order Differential Equations. Journal of Science and Technology, 14(2), 52–65. https://doi.org/10.30880/jst.2022.14.02.006
BPS Indonesia. (2025). Statistik Tanaman Perkebunan Tahunan Indonesia 2024 (Kelapa Sawit, Kopi, Kakao, Karet, Teh, dan Komoditas Perkebunan Unggulan). Bps.Go.Id. https://www.bps.go.id/id/publication/2025/08/29/8d2a6ab3510f9828daf73191/statistik-tanaman-perkebunan-tahunan-indonesia-2024--kelapa-sawit--kopi--kakao--karet--teh--dan-komoditas-perkebunan-unggulan-.html
Hidayat, T., & Fajri, S. (2026). Analysis of Public Policies in Indonesia’s Palm Oil Industry to Support Sustainable Development. ARISTO, 14(2), 998–1020. https://doi.org/10.24269/ars.v14i2.12130
Ibrahim, Z. B., Ijam, H. M., Aksah, S. J., & Abd Rasid, N. (2024). Enhancing accuracy and efficiency in stiff ODE integration using variable step diagonal BBDF approaches. Malaysian Journal of Fundamental and Applied Sciences, 20(5), 1083–1100. https://doi.org/10.11113/mjfas.v20n5.3530
Isra, M., Sitepu, A. C., Panjaitan, A. I. A., Juniarta, M. R., Safitri, N., Hasibuan, U. B., & Harahap, V. W. L. (2025). Masa Depan Industri Kelapa Sawit Tantangan dan Peluang. Jurnal Riset Rumpun Ilmu Tanaman, 4(1), 275–282. https://doi.org/10.55606/jurrit.v4i1.5673
Mahtta, R., Fragkias, M., Güneralp, B., Mahendra, A., Reba, M., Wentz, E. A., & Seto, K. C. (2022). Urban land expansion: the role of population and economic growth for 300+ cities. Npj Urban Sustainability, 2(1), 5. https://doi.org/10.1038/s42949-022-00048-y
Maulidya, A., Damayanti, A., Indra, T. L., & Dimyati, M. (2021). Prediction of Land Change for Oil Palm Plantations in Penajam Subdistrict, Penajam Paser Utara Regency, East Kalimantan Province. Journal of Physics: Conference Series, 1811(1), 12072. https://doi.org/10.1088/1742-6596/1811/1/012072
Okarda, B., Purnomo, H., Juniyanti, L., Kusumadewi, S. D., & Nadhira, S. (2024). Indonesian palm oil towards sustainability: a system dynamic approach. IOP Conference Series: Earth and Environmental Science, 1379(1), 12037. https://doi.org/10.1088/1755-1315/1379/1/012037
Olatunji, O., & Akeju, A. (2025). Comparative Analysis of Euler and Order Four Runge-Kutta Methods in Adams-Bashforth-Moulton Predictor-Corrector Method. International Journal of Mathematics, Statistics, and Computer Science, 3, 276–293. https://doi.org/10.59543/ijmscs.v3i.10471
Raming, I., Wirawan, A. S., Putri, A. A., Sahputra, D. R., Dala, M. A. D., & Avrilia, M. P. (2024). Pemetaan pertumbuhan penduduk di Kota Samarinda melalui pemodelan logistik dengan metode Adams-Bashforth-Moulton. Journal of Mathematics, Computations and Statistics, 7(1), 133–143. https://doi.org/10.35580/jmathcos.v7i1.1945
Robbaniyyah, N. A., Hanani, H., Ihwani, I. L., Awalushaumi, L., & Romdhini, M. U. (2025). Numerical Estimation of Rice Harvest Growth in West Lombok regency Using Adam-Bashforth-Moulton Method. IOP Conference Series: Earth and Environmental Science, 1493(1), 12010. https://doi.org/10.1088/1755-1315/1493/1/012010
Romadhon, R. (2025). Prediksi Produksi Daging Ayam Kampung di Provinsi Lampung Menggunakan Metode Adams-Bashforth-Moulton dan Metode Milne-Simpson. Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Lampung.
Shahputra, M. A., & Zen, Z. (2018). Positive and negative impacts of oil palm expansion in Indonesia and the prospect to achieve sustainable palm oil. IOP Conference Series: Earth and Environmental Science, 122, 12008. https://doi.org/10.1088/1755-1315/122/1/012008
Suryani, I., Suprianto, A., & Wartono, W. (2023). Solusi Numerik Model Verhulst Pada Estimasi Pertumbuhan Produksi Padi Menggunakan Metode Milne-Simpson dan Metode Adams-Bashforth-Moulton. Jurnal Sains Matematika Dan Statistika, 9(1), 27–36. https://doi.org/10.24014/jsms.v9i1.19694
Susanti, A., & Maryudi, A. (2016). Development narratives, notions of forest crisis, and boom of oil palm plantations in Indonesia. Forest Policy and Economics, 73, 130–139. https://doi.org/10.1016/j.forpol.2016.09.009
Tukan, H., Maure, O. P., & Ina, K. T. D. (2024). Analisis Keakuratan Metode Numerik Dalam Menyelesaikan Turunan Persamaan Nonlinier. Leibniz: Jurnal Matematika, 4(2), 10–22. https://doi.org/10.59632/leibniz.v4i02.446
Downloads
Published
Issue
Section
Categories
License
Copyright (c) 2026 Journal of Mathematics Education and Science

This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.
Authors who publish with this journal agree to the following terms:
- Authors retain copyright and grant the journal right of first publication with the work simultaneously licensed under a Creative Commons Attribution License that allows others to share the work with an acknowledgment of the work's authorship and initial publication in this journal.
- Authors are able to enter into separate, additional contractual arrangements for the non-exclusive distribution of the journal's published version of the work (e.g., post it to an institutional repository or publish it in a book), with an acknowledgment of its initial publication in this journal.
- Authors are permitted and encouraged to post their work online (e.g., in institutional repositories or on their website) before and during the submission process, as it can lead to productive exchanges, as well as earlier and greater citation of published work
PDF Download: 17












