%
% << nc >> 
%
function xdot=nc(t,x)
% noeud-col
global mu alpha seuil
if abs(x)> seuil
	xdot=0;
else
	xdot=mu-alpha*x.^2; 
end
return
