Gillespie Algorithm for Turing patterns on networks

3 次查看（过去 30 天）
Alexander V.
Alexander V. 2020-1-25
Hello everyone,
I tried to implement the Gillespie algorithm for the study of Turing patterns from a three-species model on a network structure. However, the simulation results do not match my expectations. First, a pattern also occurs if diffusion is switched off (epsilon =0). Second, the resulting pattern / node differentiation seems to be purely random and does not show any similarity to the deterministic (temporally oscillating) pattern for large system sizes O. I don't see were my implementation is wrong....
The (deterministic) model was taken from here. Maybe I made a mistake with the derivation of the reaction rates? Or maybe something is wrong with the for-loops? The code is given in the attached file.
I would be very grateful for any help!
Alex

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Modeling 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by