А вот то второе что вы назвали я что то не уверен что интеграл даже в главном значении сойдется.
Похоже, действительно, строгого решения не существует.
Morkonwen, ну Вы же знаете, что уравнение

не из пальца высосано. Поэтому решение со смыслом обязано быть.
Поставим такой вопрос: как изменится правая часть уравнения

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

по отношению к быстро осциллирующим функциям. Он их банально дифференцирует. Это, например, можно понять из его матричного представления. Поэтому искажения дельта-функций правой части не вылезут за пределы своих бесконечно малых окрестностей. Теперь начнем уменьшать вес дельта-функций частного решения. По тем же причинам можно утверждать,
что найдется такое размазывание этих дельта-функций, при котором в правой части произойдут изменения только в указанных бесконечно малых окрестностях. В целом, будет уменьшаться и вес правой части. В пределе можно подойти к нулевому весу правой части, но, принципиально, останутся сингулярности и решения, и правой части на концах отрезка интегрирования.
Таким образом, решение уравнения с описанными сингулярностями в правой части имеет вид

,
где

при

и

.
Функцию

можно задать бесконечным числом способов (как и дельта-функцию). Вот, например, один вид решения
![$\tau(x)= \frac {\tau_0} {\sqrt[n] {1-x^2}}$ $\tau(x)= \frac {\tau_0} {\sqrt[n] {1-x^2}}$](https://dxdy-03.korotkov.co.uk/f/6/c/e/6ce1b39446192c68eeef843f210afb5982.png)
, где

.
Можете придумать свой вид. И нет особого смысла, искать конкретный вид правой части. Так как она заведомо слабее правой части исходного уравнения

. Главное, что теперь можно записать его общее решение.