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