Мой способ использует физические соображения как наводящие, но потом всё будет обосновано.
1 часть (эвристическая). Возьмём три точки:
источник с координатами
![$(0,L), \; L>0$ $(0,L), \; L>0$](https://dxdy-03.korotkov.co.uk/f/6/0/c/60c06fa7460b02c2d70005da6e50171382.png)
;
подвижную точку
![$(x, y(x))$ $(x, y(x))$](https://dxdy-02.korotkov.co.uk/f/1/d/d/1dd74b74121ae65818b46a57d17e599d82.png)
на отражателе конечных размеров;
фокус с координатами
![$(0,F)$ $(0,F)$](https://dxdy-04.korotkov.co.uk/f/b/2/4/b24a139a42fdb20337edf4f7b5e1d64982.png)
.
Пусть
![$L$ $L$](https://dxdy-02.korotkov.co.uk/f/d/d/c/ddcb483302ed36a59286424aa5e0be1782.png)
настолько велико, что отрезок «источник — точка на отражателе» почти параллелен оси ординат. Тогда его длина равна
![$L-y(x)$ $L-y(x)$](https://dxdy-01.korotkov.co.uk/f/4/9/1/491e0ec4b685f407fb4eee7f595ac61f82.png)
.
Длина отрезка «точка на отражателе — фокус»
![$r(x)=\sqrt{x^2+(y(x)-F)^2}$ $r(x)=\sqrt{x^2+(y(x)-F)^2}$](https://dxdy-01.korotkov.co.uk/f/8/4/a/84a794e3d2877f2da0134b22e4356c2582.png)
.
Исходя из физики, их сумма не меняется при движении точки по отражателю:
![$r(x)+L-y(x)=\operatorname{const}$ $r(x)+L-y(x)=\operatorname{const}$](https://dxdy-03.korotkov.co.uk/f/a/9/3/a93b4e3195a3eb95d0a8931ca6c87dfd82.png)
, или
![$r(x)-y(x)=\operatorname{const}.$ $r(x)-y(x)=\operatorname{const}.$](https://dxdy-01.korotkov.co.uk/f/4/5/1/451b86f8cbd56c728eda6bbbddb1cf2a82.png)
2 часть. Покажем, что последнее условие приводит к Вашему уравнению. Дифференцируем его по
![$x$ $x$](https://dxdy-04.korotkov.co.uk/f/3/3/2/332cc365a4987aacce0ead01b8bdcc0b82.png)
:
![$(r-y)'=\dfrac{x+(y-F)y'}{r}-y'=0$ $(r-y)'=\dfrac{x+(y-F)y'}{r}-y'=0$](https://dxdy-01.korotkov.co.uk/f/8/a/b/8ab74d5b26b0c8b831bda426ccd0919b82.png)
![$x+(y-F)y'=ry'$ $x+(y-F)y'=ry'$](https://dxdy-03.korotkov.co.uk/f/2/8/3/2838f1a694f1830342bba2036636d57c82.png)
Возводим в квадрат:
![$x^2+2x(y-F)y'+\begin{xy}*{(y-F)^2y'^2};p+UL;+DR**h@{-}\end{xy}=x^2y'^2+\begin{xy}*{(y-F)^2y'^2};p+UL;+DR**h@{-}\end{xy}$ $x^2+2x(y-F)y'+\begin{xy}*{(y-F)^2y'^2};p+UL;+DR**h@{-}\end{xy}=x^2y'^2+\begin{xy}*{(y-F)^2y'^2};p+UL;+DR**h@{-}\end{xy}$](https://dxdy-04.korotkov.co.uk/f/7/d/6/7d696a181f85a7256615e294500c2ce782.png)
![$2(y-F)y'=x(y'^2-1)$ $2(y-F)y'=x(y'^2-1)$](https://dxdy-03.korotkov.co.uk/f/6/1/5/6155bddabf4c4c51668c75f1e7efe5c082.png)
и всё.
Вам остаётся все действия 2 части проделать в обратном порядке и получить
![$r\pm y=\operatorname{const}$ $r\pm y=\operatorname{const}$](https://dxdy-04.korotkov.co.uk/f/b/5/f/b5f7c14498e8ea839af39cf2c8a6768582.png)
.