Предложенные схемы переносяться на случай приближенного решения систем линейных дифференциальных уравнений с постоянными коэффициентами. На их основе программно анализируется устойчивость решений в смысле Ляпунова. Характер устойчивости выявляется непосредственно по матрице коэффициентов без учета информации о корнях ее характеристического многочлена.