Вычислить двойной интеграл от функции f=f(x, y) по множеству M={(x, y)

Вычислить двойной интеграл от функции f=f(x,y) по множеству M={(x,y): a<x<b, p(x)<y<q(x)} используя метод повторного интегрирования и квадратурные формулы для вычисления одномерных интегралов.
Кратный интеграл сводиться к повторному интегрированию одномерных интегралов [1, Старченко А.В., Берцун В.Н. стр. 151-152], для вычисления которых используются одномерные квадратурные формулы средних прямоугольников [1, Старченко А.В., Берцун В.Н. стр. 142-143] и трапеции [1, Старченко А.В., Берцун В.Н. стр. 144-145].
1. Написать программу, реализующую скалярный алгоритм вычисления интеграла при заданном числе N разбиения отрезка интегрирования (одинаковое в двух направлениях).
2. Преобразовать скалярную программу средствами OpenMP [2, Антонов А.С.] в программу, реализующую параллельный алгоритм для произвольного числа R скалярных процессоров (работников).
3. Провести численные расчёты для возрастающего набора значений чисел разбиения отрезка интегрирования с замером времени T. Составить таблицы с зависимостью времени от числа разбиений и числа работников. Для контроля получаемых результатов использовать две квадратурные формулы.