%
% << ho >> 
%
function zdot=ho(t,z)
% hopf
global mu alpha omega seuil
x=z(1); y=z(2);
r2=(x*x+y*y);
if r2 > seuil*seuil
	zdot(1)=0; zdot(2)=0;
else
	ar2=alpha*r2;
	zdot1=mu*x-omega*y-ar2*x;
	zdot2=mu*y+omega*x-ar2*y; 
	zdot=[zdot1;zdot2];	
end
return
