В работе рассматриваются вопросы применения метода объёмных интегральных уравнений для расчёта магнитных систем. В пакете программ GFUN, использующем интегральную постановку задачи магнитостатики, для дискретизации уравнений используются метод коллокаций и кусочно-постоянная аппроксимация неизвестных в пределах элемента разбиения области. Ограниченность применимости данного подхода связана с сингулярностью ядра интегральных уравнений.