%
% << fo >> 
%
function xdot=fo(t,x)
% fourche
global mu alpha seuil
if abs(x)> seuil
	xdot=0;
else
	xdot=mu*x-alpha*x^3; 	
end
return
