Pengintegralan numerik adalah cabang ilmu komputasi yang bertujuan menghitung nilai integral suatu fungsi ketika solusi analitik tidak tersedia atau terlalu rumit untuk dihitung secara manual. Metodemetode numerik menggantikan proses integrasi kontinu dengan operasi aljabar yang dapat dijalankan oleh komputer.
Berbagai bidang ilmu memerlukan integral, antara lain fisika (hukum gravitasi, dinamika fluida), teknik (analisis rangkaian listrik, mekanika struktural), ekonomi (model pertumbuhan, opsi finansial), dan biologi (model populasi). Pada banyak kasus fungsi yang terlibat bersifat nonlinier, memiliki singularitas, atau didefinisikan secara empiris (misalnya data eksperimen). Pada situasi inilah pengintegralan numerik menjadi solusi yang praktis.
Integral tentu didefinisikan sebagai limit dari jumlah tak terhingga area kecil di bawah kurva. Pengintegralan numerik mendekati limit tersebut dengan membagi interval integrasi [a, b] menjadi n subinterval berukuran tetap atau tidak tetap, kemudian menjumlahkan nilai fungsi pada titiktitik khusus dalam setiap subinterval.
Jika I = _a^b f(x)dx, maka pendekatan numerik umumnya memiliki bentuk:
I w_if(x_i)
di mana x_i adalah titik evaluasi dan w_i adalah bobot yang ditentukan oleh metode tertentu.
Terbagi menjadi tiga jenis utama: leftendpoint, rightendpoint, dan midpoint. Hanya cocok untuk fungsi yang halus dan interval kecil.
Mengganti area di bawah kurva dengan trapesium. Rumus dasar untuk satu subinterval:
I (ba)[f(a)+f(b)]/2
Dengan n subinterval, hasilnya menjadi penjumlahan semua trapesium.
Memanfaatkan polinomial kuadratik pada setiap pasangan subinterval. Formula umum (untuk n genap):
I (h/3)[f(x0) + 4 f(x_odd) + 2 f(x_even) + f(xn)]
di mana h = (ba)/n. Simpson biasanya memberikan akurasi yang jauh lebih tinggi dibanding trapesium untuk fungsi yang halus.
Memilih titiktitik x_i dan bobot w_i secara optimal sehingga integrasi dari polinomial derajat hingga 2n1 dapat dihitung secara tepat. Contoh paling terkenal adalah GaussLegendre dengan titiktitik akar polinomial Legendre.
Berbasis statistik: mengestimasi integral dengan mengambil sampel acak pada domain integrasi. Sangat berguna untuk integrasi berdimensi tinggi, di mana metode deterministik menjadi tidak praktis.
Setiap metode memiliki estimasi kesalahan yang bergantung pada turunan fungsi yang tidak diketahui. Misalnya, kesalahan Trapesium dapat dituliskan sebagai:
|E_T| (ba)^3 / (12 n^2) max|f''()|
Dimana berada di dalam [a,b]. Dengan menambah n, kesalahan berkurang secara kuadratik. Simpson memiliki orde kesalahan O(h^4), sehingga konvergen lebih cepat bila fungsi cukup halus.
Berikut contoh sederhana implementasi Metode Trapesium dalam JavaScript:
function trapesium(f, a, b, n) { const h = (b - a) / n; let sum = f(a) + f(b); for (let i = 1; i < n; i++) { sum += 2 * f(a + i * h); } return (h / 2) * sum;}// contoh penggunaanconst f = x => Math.sin(x);const hasil = trapesium(f, 0, Math.PI, 1000);console.log('Integral sin(x) dari 0 sampai ', hasil); Konsep yang sama dapat diterapkan dalam Python, MATLAB, atau C++ dengan menyesuaikan sintaks.
Pengintegralan numerik adalah alat penting bagi ilmuwan dan insinyur modern. Dengan memahami prinsip dasar, pilihan metode, serta cara memperkirakan kesalahan, pengguna dapat menyelesaikan masalah integrasi yang sebelumnya tak terjangkau. Selalu lakukan verifikasi dengan contoh yang sudah diketahui nilai tepatnya, dan pilih metode yang seimbang antara kebutuhan akurasi dan sumber daya komputasi.
