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