Masalah perpindahan dan perombakan pencemar dalam air tanah dapat dimodelkan dalam bentuk sistem persamaan diferensial parsial parabolik nonlinier. Perombakan yang dimaksud adalah menghilangkan atau mengurangi pencemar sehingga tidak membahayakan bagi lingkungannya. Variable-variabel dalam sistem persamaan diferensial parsial tersebut teridiri dari pencemar, oksigen, dan mikroorganisme. Dalam sistem ini menunjukkan adanya interaksi antara pencemar, oksigen dan mikroorganisme. Interaksinya adalah dengan adanya oksigen, mikroorganisme akan berkembang biak dan mendegradasi pencemar sehingga konsentrasi pencemar akan berkurang dari keadaan sebelumnya. Untuk menyelesaikan sistem persamaan diferensial parsial tersebut dilakukan diskritisasi. Diskritisasi ini menggunakan metoda beda hingga (finite difference). Diskritisasi dilakukan pada variabel spasial persamaan diferensial parsial, yang menghasilkan sistem persamaan diferensial biasa. Sistem persamaan diferensial biasa yang dihasilkan adalah sistem persamaan diferensial biasa nonlinier yang berukuran besar. Penyelesaian sistem persamaan diferensial biasa di atas dilakukan dengan pendekatan numerik. Integrator yang digunakan adalah metoda multistep, prediktor-korektor. Beberapa nilai solusi awal diperlukan disediakan oleh metoda Rung-Kutta Implisit Diagonal dari Viz. Penyelesaian ini menggunakan besar langkah yang adaptif. Besar langkah adaptif diperlukan untuk mengatasi besar kesalahan lokal (local trancation error) yang berubah-ubah pada setiap langkahnya. Implementasi dilakukan dengan menggunakan bahasa pemrograman FORTRAN 77. Eksperimen dilakukan pada komputer pribadi dengan clock rate 120 Mhz, dengan nilai awal pencemar berkonsentrasi tinggi di bagian tengahnya yang ditunjukkan oleh ekspresi 4* v [1-x*x-y*y] , oksigen dan mikroorganisme berkonsentrasi merata (homogen) pada seluruh penampang (aquifer). Konsentrasi oksigen dan mikroorganisme masing-masing adalah 4,5 dan 0,1. Waktu pengamatan selama 130 hari dan area terbagi atas grid 10x10 dengan toleransi 10^-4. Hasil komputasi penyelesaian dengan metoda prediktor-korektor menunjukkan hasil yang efektif, dan waktu komputasi untuk menyelesaiakan tiga titik pertama (metoda Runge-Kutta Implisit Diagonal) mendekati dua kali dari waktu yang digunakan oelh metoda prodiktor-korektor untuk mnyelesaiakn titik berikutnya sampai titik terakhir (tf).