solution

Tính tích phân kép $I=\int_{a_1}^{b_1}\int_{a_2}^{b_2} (x^2+y^2)\,dy\,dx$ trên miền chữ nhật $[a_1,b_1]\times[a_2,b_2]$ bằng Simpson 1/3 kép: với mỗi mốc $x_i$ trên $[a_1,b_1]$ ($n_x$ đoạn con, $n_x$ chẵn), áp dụng Simpson 1/3 theo $y$ với $n_y$ đoạn con ($n_y$ chẵn) để tính $g(x_i)=\int_{a_2}^{b_2}(x_i^2+y^2)\,dy$ xấp xỉ, sau đó áp dụng Simpson 1/3 theo $x$ trên các giá trị $g(x_i)$ vừa tính.