function dxdt=sch_ode(t, x, mu)
% Schnakenberg egyenlet
dxdt=zeros(2,1);
dxdt(1)=1/8-x(1)+x(1)^2*x(2);
dxdt(2)=mu-x(1)^2*x(2);
end