Całki podwójne po przystających trójkątach

[ferrm]

Może to nie jest czysta analiza, ale mam pytanie. Jak wygląda algorytm do wyznaczania \(\displaystyle{ \iint\limits_D f(x,y)dxdy}\), gdzie D jest trójkątem z zastosowaniem złożonej kwadratury (trójkąt dzielimy na \(\displaystyle{ n^2}\) trójkątów przystających) rzędu drugiego (środki ciężkości)