Komputasi untuk Sains dan Teknik - Universitas Indonesia

7.2. METODE RUNGE KUTTA 111terlibat dalam perhitungan ini. Dimulai menghitung k 1 kembalik 1 = hf(q 1 )= h(m 1 − q 1 m 2 )= 0, 1((1, 5 × 10 −5 ) − (0, 14814 × 10 −5 )(0, 25))= 0, 14630 × 10 −5lalu menghitung k 2k 2 = hf(q 1 + k 12 )= h[(m 1 − (q 1 + k 12 )m 2)]= 0, 1[(1, 5 × 10 −5 − ((0, 14814 × 10 −5 ) += 0, 14447 × 10 −50, 14630 × 10−5)(0, 25)]2dilanjutkan dengan k 3k 3 = hf(q 1 + k 22 )= h[(m 1 − (q 1 + k 22 )m 2)]= 0, 1[(1, 5 × 10 −5 − ((0, 14814 × 10 −5 ) += 0, 14449 × 10 −50, 14447 × 10−5)(0, 25)]2kemudian k 4k 4 = hf(q 1 + k 3 )= h[(m 1 − (q 1 + k 3 )m 2 )]= 0, 1[(1, 5 × 10 −5 − ((0, 14814 × 10 −5 ) + 0, 14449 × 10 −5 )(0, 25)]= 0, 14268 × 10 −5akhirnya diperoleh q 2q 2 = q 1 + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 )= 0, 14814 × 10 −5 + 1 (0, 14630 + 2(0, 14447) + 2(0, 14449) + 0, 14268) × 10−56= 0, 29262 × 10 −5Dengan cara yang sama, q 3 , q 4 , q 5 dan seterusnya dapat dihitung. Tabel di atas menunjukkanhasil perhitungannya. Kolom q exact diperoleh dari persamaan (7.15).Luar biasa!! Tak ada error sama sekali. Mungkin, kalau kita buat 7 angka dibelakang koma,errornya akan terlihat. Tapi kalau anda cukup puas dengan 5 angka dibelakang koma, hasil ini