Анализ дискретных динамических систем

[PAK]

Доброго времени суток!
Интересуюсь следующей проблемой, которую, формально, можно поместить и в раздел Computer Science (на усмотрение модераторов).
Даны две вещественные квадратные матрицы $N \times N$: $S_x$ и $S_y$. К слову, эти матрицы суть $x$-овая и $y$-овая компоненты некоторого двухмерного векторного поля (на физическом языке это звучит как усредненное поперечное поле Пойнтинга), вычисленного в фиксированных узлах квадратной сетки соответствующей размерности.
Собственно, задача состоит в том, чтобы, имея лишь матрицы $S_x$ и $S_y$ восстановить структуру (фазовый портрет) следующей двумерной динамической системы:

\[ \dfrac{\mathrm{d}}{\mathrm{d}\,\tau}\,\boldsymbol{r}=\boldsymbol{S}(\boldsymbol{r}), \]

где $\boldsymbol{r}=(x,y)^{T}$ и $\boldsymbol{S}=\big(S_x(\boldsymbol{r}),S_y(\boldsymbol{r})\big)^{T}$. Необходимо также уметь выделять особые траектории и точки ДС такие, как сепаратрисы, предельные циклы, фокусы и седла.
Может быть кто-нибудь подскажет, есть ли какие-либо библиотеки (лучше на Фортране), где можно отыскать процедуры для подобного анализа «дискретных» динамических систем. Если ничего похожего нет, посоветуйте с чего начинать.

Заранее благодарен