
- точечный источник сферической волны. Вокруг него нарисуем сферу

радиуса

.
Мысленно поместим еще источник в точку

. И тоже окружим сферой

, нового радиуса

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

выполнено

,

, где

- расстояния между центрами сфер.
Натянем еще одну сферу, содержащую внутри обе эти сферы. Назовем ее

и радиус ее соответственно

.
Теперь проинтегрируем дивергенцию имеющегося выражения по внутренности поверхности, ограниченной сферой

и не содержащей сферы

и

, т.е.

Дальше надо разбить на три интеграла и два из них оценить,
радиус

устремить к бесконечности
радиус поверхности от воображаемого источника к нулю
Отсюда и получится Ваш принцип Гюйгенса-Френеля.

интеграл берется по поверхности

, хотя в доказательстве мы говорили о сфере, это вовсе не обязательно. Сферу вокруг воображаемого источника мы стянули в точку, а большую сферу наоборот растягивали на бесконечность, поэтому считать их сферами или нет тоже не имеет большого значения.