Обобщенная квадратурная формула для кратных интегралов

Излагается метод приближенного вычисления для кратных неберущихся интегралов. Рассуждения, ради простоты, проводятся для двойных интегралов
$$ \int_a^b\int_c^df(x,y)g(x,y)\,dx\,dy, $$
причем $f$, $g$ удовлетворяют определенным условиям. Дается точная оценка погрешности метода и устанавливается сходимость процесса.