The Brusselator is a simple model for oscillatory chemical systems such as the Belousov-Zhabotinski reaction. The time evolution of the concentration of two chemical species, x and y, is described by the ODE's: dx/dt = A+x2y - (B+1)x AND dy/dt = Bx - x2y,
Where A>0, B>0 are constants.
[A]. Find the single steady state of this model (do this portion on paper and attach to solution).
[B]. Write a program (MATLAB) to compute x(t) and y(t) for a given initial condition x(0)≥0 and y(0) ≥0. Plot a few trajectories in the xy plane, and mark the location of the steady state. Investigate cases where B/(1+A2) is less than, greater that, and equal to one.
Please make "x vs t" and "y vs x" plots for three scenarios:
A=0.5, B=1, x_initial=0.1, y_initial=2
A=1, B=3, x_initial=1, y_initial=1
A=2, B=5, x_initial=1, y_initial=1
In all scenarios please use Nstep=10000 and tau=0.01.