WOLFRAM NOTEBOOK

WOLFRAM|DEMONSTRATIONS PROJECT

Bifurcation Diagram for the Three-Variable Autocatalator

bifurcation parameter μ

0.15

This Demonstration shows how one can compute the bifurcation diagram for a nonlinear chemical system such as the three-variable autocatalator (see Details).

Indeed, in order to obtain such a diagram, one has to locate the maxima of the time series for all values of the bifurcation parameter,

, which can be readily done using Mathematica's built-in function FindArgMax.

The Demonstration illustrates the dynamics of the concentrations

, and

for various values of the bifurcation parameter

. The time series option gives a plot of

logγ

versus time and shows the loci of all maxima.

Try the following values of

: 0.1, 0.14, 0.15, 0.151, and 0.155 to observe period 1, 2, 4, 8, and 5 behaviors, respectively. For

μ=0.153

, chaotic behavior occurs. When

is large enough, you can observe a reversed sequence leading back to period 1 behavior. These results are confirmed by the bifurcation diagram (a remerging Feigenbaum tree), first given in [1] and reproduced in the present Demonstration.

You are using a browser not supported by the Wolfram Cloud

Supported browsers include recent versions of Chrome, Edge, Firefox and Safari.

I understand and wish to continue anyway »