Численное решение систем линейных обыкновенных дифференциальных уравнений с постоянными коэффициентами

Рассматривается метод решения систем линейных обыкновенных дифференциальных уравнений (о.д.у.) с постоянными коэффициентами, основанный на вычислении матричной экспоненты и интеграла от нее с помощью рекуррентных соотношений. Приводятся формулы для вычисления глобальной погрешности численного решения.