Добрый день.
Передо мной встала задача вычисления значений функции определяемой рекурсивно через интегрирования самой себя порядком ниже.
В упрощеном виде всё выглядит примерно так:

Функции

и

- обычные, не рекурсивные. Индекс

у

призван показать что для разных

она вычисляется по-разному.
Интегралы не берутся, значения нахожу численно, квадратурными формулами.
У меня есть с чем сравнивать полученные результаты (некоторые таблицы значений данной фунции представлены в литературе, также я посчитал её методом Монте-Карло) и результаты выглядят корректно.
Однако у меня вы получается оценить погрешность вычисления.
Предположим, для определенности, что мы интегрируем методом средних прямоугольников. Оценить погрешность при вычислении лишь самого верхнего в рекурсии интеграла будет неверным - ведь значение функции в каждом из узлов на которые разбиваем интервал интегрирования так же получено с погрешностью.
Оценим погрешность уровня

как

где
![$M_n = \max_{[a, b]}{f''_n}$ $M_n = \max_{[a, b]}{f''_n}$](https://dxdy-01.korotkov.co.uk/f/8/3/0/830214f206fbd274bb6eca6f9672a58082.png)
а

, где

- количество узлов.
Но значение функции в каждом из

узлов получено с погрешностью

, что дает суммарную погрешность от вычисления функции в узлах

И т.д.
Каким же образом верно оценить погрешность?
Перемножать все эти погрешности, видимо, будет неверным, сделав это получим огромное число в

Насколько обосновано в таком случае пользоваться методом двойного счета?
Подсчитав интегралы с

и

узлами на
каждом уровне имеем ли мы право считать относительную погрешность как

?
Посоветуйте что можно предпринять в данной задаче.
Спасибо.