Показано применение геометрических методов электродинамики к расчету оптических устройств, таких как линзы Максвелла и Люнеберга. Рассмотрено уравнение эйконала, которое с помощью метода характеристик было преобразовано к системе ОДУ. Полученная система записана для случая линзы Максвелла и Люнеберга и решена стандартными численными методами. Описаны детали реализации и приведены изображения траекторий лучей и фронтов волн.