mu = 5;
omega = 1;
[t, phi] = ode45(@(t,x)elemi(t,x,mu,omega),[0 10],[1 1]);
plot(phi(:,1),phi(:,2));