...производные в левой части пренебрежимо малы по сравнению с дельтами...
Не обязательно. Особенности искомой функции вполне могут сократиться с дельтами. В одномерном случае функция для этого должна иметь излом в начале координат.
В многомерном случае это бессмыслица.
Я пытался действовать по аналогии с одномерным случаем: проинтегрировать уравнение по малой области

, охватывающей начало координат:
![$$\[\iint\limits_{S}{\delta \left( x \right)\delta \left( y \right)f\left( x,y \right)dxdy}=f\left( 0,0 \right)\]$$ $$\[\iint\limits_{S}{\delta \left( x \right)\delta \left( y \right)f\left( x,y \right)dxdy}=f\left( 0,0 \right)\]$$](https://dxdy-02.korotkov.co.uk/f/9/c/e/9cead59f1c1243a2bfdc31c4d8fac7e682.png)
![$$\[\iint\limits_{S}{\left( \frac{{{\partial }^{2}}}{\partial {{x}^{2}}}+\frac{{{\partial }^{2}}}{\partial {{y}^{2}}} \right)f\left( x,y \right)dxdy}=\int\limits_{\partial S}{\left( -\frac{\partial f}{\partial y}dx+\frac{\partial f}{\partial x}dy \right)}\]$$ $$\[\iint\limits_{S}{\left( \frac{{{\partial }^{2}}}{\partial {{x}^{2}}}+\frac{{{\partial }^{2}}}{\partial {{y}^{2}}} \right)f\left( x,y \right)dxdy}=\int\limits_{\partial S}{\left( -\frac{\partial f}{\partial y}dx+\frac{\partial f}{\partial x}dy \right)}\]$$](https://dxdy-03.korotkov.co.uk/f/a/5/3/a53ddf1338828f3f58e6edc4ab4a28ba82.png)
В последнем интеграле по контуру

интегрируется повёрнутый на 90° градиент функции. Отличие этого интеграла от нуля означает, что функция имеет особенность типа

или

или что-то типа того, то есть, что градиент направлен от центра (или к нему) вне зависимости от того, в какое направление от начала координат отступить (острая вершина конуса или расходимость). Проблемы начинаются тогда, когда искомая функция в нуле расходится.
Неужели математики ещё не придумали, как решать такие задачи?