Пусть у нас есть замкнутая поверхность с заданной на ней функцией

, где

-поверхность.
И надо найти распределение

внутри и снаружи это поверхности (
внешняя и внутренняя задача Дирихле), которое удовлетворяет уравнению Лапласа

(или Кирхгофа, если рассматривать оптику).
Решение имеет вид

, где

- функция Грина оператора Лапласа (потенциал единичного точечного заряда).
Несмотря на подчёркнутое мной в цитате (и на название темы), я не думаю, что здесь идёт речь о решении задачи Дирихле. Я предлагаю другую интерпретацию. Написанный интеграл — потенциал двойного слоя, а функция

, где

, — его плотность. Далее,

(или

), а не функция Грина для данной области.
Автор пишет, что функция

задана, следовательно, можно найти

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

с плотностью

:

Однако ни

, ни

(предельные значения потенциала снаружи и изнутри) не равны

. Сам потенциал и плотность его источников на границе — совершенно разные величины.