Системы ДУ произвольного порядка

[FermaYails]

Доброго времени суток!
Подскажите, как решать (или где прочитать) систему ДУ вида:
$x^{(n)} = A_{n-1}x^{(n-1)} +...+ A_{0}x$.
Тут вместо степени дифференцирование. Переменные есть вектор-функции. Буквы при иксах - квадратные матрицы.
Плюс ещё решение с краевыми условиями.
Заранее благодарю!

[DeBill]

FermaYails
Ну, в принципе, задача сводится (стандартными заменами $y_1=x, y_2=x',...$) к здоровенной системе $Y'=AY$, где $A$- блочная матрица (над диагональю стоят единичные матрицы, а в нижней - блочной - строке - все те матрицы $A_j$). Фундаментальная матрица решений тогда есть $e^{At}$...Проблема - в вычислении этой экспоненты... Без каких-нить хороших ограничений на матрицы (типа - они все попарно коммутируют; или все одновременно диагонализируются...) надежд на получение компактного ответа мало. Впрочем, если размеры их невелики, и $n$ небольшо, то справится все же можно - тупо приводя $A$ к жордановой нормальной форме (или даже - найти ейные с.значения, и искать решение в виде комбинации соответствующих экспонент)

[Singular]

Если у Вас есть конкретная система уравнений, то лучше обсудить её.