Purpose: This research intends to compare two numerical method models, which are Runge-Kutta and by using HYDRUS software. Theoretical reference: One of the consequences of human population that grows rapidly is the increase of human activity, along with the increase of pollutant that flows into water bodies. Water pollution is a significant environmental issue that threatens the health of ecosystems and human populations worldwide. Method: To determine a set of parameters of the mechanism in a constructed wetland, a numerical model may be of aid to acquire parameters which represent the mechanism of pollutant removal in it. Result and Conclusion: This study yields different efficiency of pollutant removal when simulated with HYDRUS and with Runge-Kutta. The differences in the effectivity results in HYDRUS and Runge-Kutta might be caused by the detail of simulated mechanism, where HYDRUS considers the mechanism happened in soil column, while Runge-Kutta only focuses on the water column. Implication of research: In order to observe the effectivity of a constructed wetland, a numerical model may be used. The parameters needed to represent the environment and properties of both the plants and soil can be obtained by building a demonstration plot. Originality/ value: to compare two numerical method models, which are Runge-Kutta and by using HYDRUS software.