Статья посвящена разработке методов численного решения и тестирования жесткой интегро-дифференциальной задачи моделирования изоэлектрического фокусирования (ИЭФ) в так называемых «аномальных» режимах. Были решены следующие задачи: аналитическое преобразование интегро-дифференциальной задачи математической модели к виду стандартной краевой задачи, пригодной для численного решения методом Рунге—Кутта; преодоление неконтролируемого накопления вычислительной погрешности, вызванной «жесткостью задачи», путем представления решения в экспоненциальной форме и составления оптимизационных алгоритмов. Для тестирования модели разработаны два численных метода: асимптотический и метод касательных. Оба метода показали высокую степень соответствия расчетных и асимптотических решений задачи