Рассматривается кубатурная формула приближенного вычисления интегралов от быстроосциллирующей функции трех переменных с использованиемлагранжевой полиномиальной интерлинации функций с оптимальным выбором узловых плоскостей для приближения неосциллирующего множителя. На классе дифференцируемых действительных функций получена оценка погрешности кубатурной формулы.