Аннотация:Рассмотрен метод решения жёстких систем обыкновенных дифференциальных уравнений на основе вычисления матричной экспоненты. Дополнительно предложен метод контроля точности в случае линейных задач. Представлены результаты численных расчётов для линейных и нелинейных задач в зависимости от значений параметров как самих задач, так и численного метода.