Применительно к сопряжённой задаче теплопроводности и химической кинетики рассматриваются метод расщепления по физическим процессам и два варианта глобально-неявного метода: с простыми итерациями и с итерациями по дискретному методу Ньютона. Установлены особенности из применения. Показана возможность использования глобально-неявного метода с итерациями по дискретному методу Ньютона в пакетах прикладных программ общего назначения. Исходные тексты разработанных программ опубликованы в свободном доступе.