Я немного изменю обозначения. Первая точка

имеет декартовы координаты

, вторая точка

— координаты

. То есть у меня индекс — это не номер точки, а номер координаты. И ещё я не различаю случайную величину и её значение.
Пусть

— некоторая интегрируемая функция одной переменной. Запишем интеграл по четырёхмерной области

:

С помощью замены

перейдём от набора независимых переменных

к набору

. Якобиан перехода равен

.
Теперь с помощью замены

перейдём к набору

. Здесь якобиан перехода равен

. Очевидно,

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

От порядка зависят пределы интегрирования. В общем случае они в каждом интеграле определяются значениями предыдущих переменных интегрирования. У нас

. Область, которую пробегает точка

при фиксированных

, обозначим

. Она определяется условиями

. Это та самая область, о которой я спрашивал
здесь.
С учётом всего этого

Внутренний двойной интеграл равен площади области

. Вы нашли, что эта площадь равна
![$$R^2\left[2\arccos(s/2R) - \frac{s}{R}\sqrt{1-\frac{s^2}{(2R)^2}} \right]=\frac{\pi R^4}{2s}f_l(s)$$ $$R^2\left[2\arccos(s/2R) - \frac{s}{R}\sqrt{1-\frac{s^2}{(2R)^2}} \right]=\frac{\pi R^4}{2s}f_l(s)$$](https://dxdy-03.korotkov.co.uk/f/6/b/f/6bf1befdbdf88200c0808776c2351bb382.png)
Подставляем это вместо двойного интеграла, много чего сокращается, получаем

Остался один шаг, где надо использовать смысл интегралов (2), (3) и произвольность функции

.