function dx = elemi(t,x,mu, omega)
dx = zeros(2,1);
dx(1) = mu-x(1)^2;
dx(2) = -1*x(2);


