Численная интеграция с c ++ на заданной сетке с фиксированной дискретизацией констант

У меня есть следующая проблема:

Мой код на C ++ может вычислить две функции

f1 (i1, i2, i3, i4)

f2 (j1, j2)

для каждого набора {i1, i2, i3, i4 } Я получаю некоторое значение f1, и для каждого набора {j1, j2} я получаю некоторое значение f2.

наборы {i1, i2, i3, i4} и {j1, j2} задаются на ФИКСИРОВАННОЙ сетке с некоторым постоянным шагом дискретизации "h".

Мне нужно вычислить на математическом языке интеграл F3 (x1, x3) = Integral [f1 (x1, x2, x3, x4) * f2 (x3, x4) dx3 dx4]

Простое суммирование недостаточно, так как у f2 много прыжков.

Есть ли какая-нибудь библиотека C ++, которая может выполнять такую ​​интеграцию? Или какой-нибудь алгоритм, который легко реализовать (я не очень хорошо разбираюсь в c ++)

большое спасибо