Метод Симпсона для многомерных интегралов

[aaastalm]

Здравствуйте возник небольшой вопрос. Метод Симпсона очень прост для одномерных интегралов. Так вот возник вопрос, как его использовать для многомерных интегралов. Для одномерного случая интегральная формула при разбиении на 2N частей имеет вид:

$$I=b-a/6 N (f_0+4(f_1+f_3+..+f_{2N-1})+2(f_2+f_4..+f_{2N-2})+f_{2N})$$

Понятно, что для многомерного интеграла надо брать разбиение на

$$(2N)^k$$

частей. Но как при этом будет выглядеть формула? А точнее интересует какие значения функции будут считаться четными, а какие нет? А точнее какие из них умножать на четыре, а какие на два? Если есть какая литература по этому методу посоветуйте пожалуйста.