Розглянуто оптимальні алгоритми для реалізації кубатурної формули Сімпсона з застосуванням принципу подвійного перерахунку під час обчислення багатократних інтегралів. Порівняно запропонований алгоритм з вбудованими функціями пакета розширень системи комп'ютерної математики (СКМ) на тестовому прикладі обчислення інтегральних тригонометричних функцій. Розширено функціональні можливості використання СКМ Mathematica та Maple.