В работе развит алгоритм, основанный на векторном методе конечных элементов для численного моделирования квазистационарных уравнений Максвелла. Этот алгоритм позволяет естественным образом учесть непрерывность тангенциальных компонент электрического поля на границах сред с различными физическими свойствами. Представлены результаты тестовых расчетов.