g______d, если не ошибаюсь, общая формула многочлена получается такая:

, где

- биномиальные коэффициенты,

. Промоделировал, полиномы получаются такие же, как и при решении системы линейных уравнений "в лоб", только теперь можно рассчитать до степени порядка

, дальше коэффициенты выходят за диапазон. Насчет сходимости - посмотрел изменение коэффициента при первой степени аргумента - монотонно возрастает сильнее чем логарифм.
UPD насчет сходимости. У меня сильное подозрение, что эти полиномы сходятся к идеальной ступеньке. Рассмотрев подынтегральную функцию при увеличении

, и разбив интеграл на сумму двух по областям
![$[0, \varepsilon] $[0, \varepsilon]](https://dxdy-02.korotkov.co.uk/f/1/c/9/1c9e1c5b57326b291cd8b45105c2221d82.png)
и
![[\varepsilon, x]$ [\varepsilon, x]$](https://dxdy-03.korotkov.co.uk/f/a/3/3/a338561b4320f9f8b3a591bda16f310582.png)
(рассматривая без уменьшения общности положительные аргументы), можно показать, что для любого

существует такое

, что отношение интеграла по второй области к интегралу по первой меньше любого наперед заданного числа. Отсюда следует, что (при нормировке интеграла) для любого

существует такое

, что значение интеграла будет отличаться от единицы на любую заранее заданную малую величину уже на интервале
![$[0, \varepsilon]$ $[0, \varepsilon]$](https://dxdy-04.korotkov.co.uk/f/3/2/8/3285dfee1ba8536f367e305406fc2b9b82.png)
, из чего следует, что с увеличением

мы можем как угодно близко приблизиться к единице при любом сколь угодно малом значении аргумента - чистая ступенька в пределе.