Как было показано ранее - см. выражение (2.6), - в системе главных диэлектрических осей объемная плотность электрической энергии поля определяется двумя аналогичными выражениями:
.
Введя новые переменные
;
;
и
;
;
, где
, получим уравнения двух эллипсоидов
, (2.21)
, (2.22)
из которых первый называется лучевым (рис. 11), а второй - волновым (рис. 12) эллипсоидами Френеля. Координаты точек на поверхностях эллипсоидов представляют собой с точностью до постоянной

значения компонент векторов

и

соответственно. Отметим, что полуоси лучевого эллипсоида пропорциональны соответствующим главным скоростям распространения

,

,

а полуоси волнового эллипсоида обратно пропорциональны им.
Рассмотрим решение волнового уравнения Френеля для волны, распространяющейся вдоль одной из координатных осей, например вдоль оси
. В этом случае
,
. Из (2.17) получим
.
Решения этого уравнения очевидны:
,
. С другой стороны, из рассмотрения частных случаев мы ужe знаем, что с главной скоростью
распространяется волна, поляризованная вдоль оси
, со скоростью
- волна, поляризованная вдоль оси
. Из рис. 12 следует, что указанные направления колебаний векторов
и
совпадают с полуосями эллипса центрального сечения волнового эллипсоида Френеля плоскостью, ортогональной
, т.е. плоскостью
. Аналогично для волны, распространяющейся вдоль оси
из (2.20) для групповой скорости получим 
,
т.е.
. Направления колебаний векторов
и
совпадают с осями
и
,т.е. с полуосями эллипса центрального сечения лучевого эллипсоида (см. рис. 11). Рассмотренный алгоритм определения направлений колебания векторов
,
и
,
можно распространить и на общий случай положения векторов
и
в пространстве. Этот алгоритм формулируется следующим образом: для того чтобы определить направление колебаний векторов
при заданном значении вектора
распространения луча (волнового фронта) в кристалле, необходимо построить центральное сечение лучевого (волнового) эллипсоида Френеля плоскостью, ортогональной вектору
. Тогда направления колебаний векторов
будут совпадать с полуосями полученного эллипса сечения. Лучевые (фазовые) скорости соответствующих волн будут прямо (обратно) пропорциональны длинам полуосей.
Указанный алгоритм нагляден (см. рис. 11 и 12) и очень удобен для анализа особенностей распространения световых волн в кристаллах. Подробное доказательство его приведено в работе [2]. Как известно из аналитической геометрии, в любом эллипсоиде общего вида могут быть найдены не более двух центральных сечений, представляющих собой окружность; это означает, что в любом кристалле существуют не более двух направлений
, ортогональных соответствующим сечениям, в которых кристалл проявляет себя как изотропная среда. Такие направления называют лучевыми (волновыми) оптическими осями, а соответствующий кристалл - двухосным. При равенстве двух из трех значений главных диэлектрических проницаемостей эллипсоиды Френеля являются эллипсоидами вращения. В этом случае оптические оси вырождаются в одну оптическую ось - ось вращения, а кристалл называется одноосным.