бифуркации в дифуре и получение упорядоченных узоров

[dikiy]

Простите за несколько сумбурное название темы, но я не знаю, как это правильно на русском будет :)

В общем имее следующую систему уравнений:

\begin{equation} \begin{split}\frac{\partial n}{\partial t}=\frac{\gamma w}{1+\sigma w}n-n^2-\mu n+\Delta n\\ \frac{\partial w}{\partial t}=p-(1-\ro n)w-w^2n+D\Delta(w-\beta n) \end{split}\end{equation}
$\beta,\ro,D,\sigma,\gamma,\mu$ - фиксированные параметры. p - управляющий параметр. Анализ на несабильность показывает, что при p>0.15 должна возникать нестабильность. Я моделирую это с помощью обычного FTCS метода (в лом по Эйлеру). начальные параметры для $n$ и $w$ задаю случайным генератором. Через 200 итераций картинка стабилизируется, но никаких узоров я не наблюдаю. Что я не так делаю? Или чего я не понимаю?

вот программа для octave. Но в matlabе тоже должна работать.
http://dpaste.com/1763636/

[Lia]

Тема перемещена из форума «Математика (общие вопросы)» в форум «Помогите решить / разобраться (М)»
Причина переноса: не указана.