1. Masalah yang diberikan adalah persamaan rekursif untuk $u^{n+1}$ yang melibatkan beberapa variabel dan operator matriks. 2. Persamaan yang diberikan adalah: $$ u^{n+1} = u^n + \frac{k}{6} \left(R_u^1 + 2R_u^2 + 2R_u^3 + R_u^4\right) $$ Dimana: $$ R_u^1 = \left(\frac{\sigma^2}{2}\right) B^{-1}(Au + f_u) + \xi_\tau w - ru $$ $$ R_u^2 = \left(\frac{\sigma^2}{2}\right) B^{-1}(Aw + f_w) + \xi_\tau B^{-1}(Au + f_u) - rw $$ $$ R_u^3 = \left(\frac{\sigma^2}{2}\right) B^{-1}(Ay + f_y) + \xi_\tau B^{-1}(Aw + f_w) - ry $$ $$ R_u^4 = \left(\frac{\sigma^2}{2}\right) B^{-1}(Az + f_z) + \xi_\tau B^{-1}(Ay + f_y) - rz $$ 3. Untuk menyelesaikan persamaan ini, kita harus menghitung setiap $R_u^i$ terlebih dahulu dengan memasukkan nilai-nilai dari $u, w, y, z$ dan parameter lainnya. 4. Setelah mendapatkan nilai $R_u^1, R_u^2, R_u^3, R_u^4$, substitusikan ke dalam persamaan utama untuk mendapatkan $u^{n+1}$. 5. Proses ini biasanya dilakukan secara numerik karena melibatkan invers matriks $B^{-1}$ dan operator matriks $A$. 6. Jika nilai-nilai numerik dari $u, w, y, z, A, B, f_u, f_w, f_y, f_z, \sigma, \xi_\tau, r, k$ diberikan, maka langkah-langkah perhitungan dapat dilakukan secara langsung. 7. Secara ringkas, solusi dari persamaan ini adalah iterasi nilai $u$ menggunakan rumus rekursif di atas dengan menghitung setiap $R_u^i$ secara berurutan.