function dF=brusselator(t,x)
global A B
dF=[-B-1 x(1)^2;B -x(1)^2]*x+[A 0]';
end