Ну хорошо, попытаемся решить. Чтобы
единичный вектор

был направляющим для этой прямой, нужно, чтобы на нём достигался максимум квадратичной формы

, элементы матрицы которой суть

, а

-- это вектор, соединяющий центр масс с

-той точкой. Это означает, что

-- это собственный вектор матрицы

, отвечающий её наибольшему собственному числу. В нашем случае:

,

и

. Векторы (после сокращения на общий множитель

):




Матрица:

Её собственные числа (после сокращения на четвёрку):

и

. Собственный вектор, отвечающий наибольшему собственному числу

:

(он не нормирован, но здесь это и не обязательно -- нормировка нужна была только при формулировании алгоритма). Уравнение прямой:

Конечно, это решение не для школьников, но зато именно "геометрично" и универсально.
--------------------------------------------------------------
Да, а почему прямая проходит через центр масс -- это просто. Пусть на данной прямой достигается минимум. Развернём координатные оси так, чтобы прямая оказалась параллельной одной из осей -- например, иксов. Пусть она проходит через

(в новых координатах, естественно). Тогда сумма квадратов расстояний -- это

. Очевидно, что она минимальна именно при

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