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