function dx = bifurcation(t,x)

dx = zeros(2,1);
%state equation
mu = 1;
dx(1) = mu-x(1)^2;
dx(2) = -x(2);