function y = fd( x )
%drift function 
    y=zeros(2,1);
%     y(1)=-x(2);
%     y(2)=x(1);
%circular
    
%     y(1)=x(2);
%     y(2)=(1-x(1)^2)*x(2)-x(1);
%van der Pol

      y(1)=1;
      y(2)=x(2)*(2-x(2));
      if (y(2)<0) y(2)=0;
 %fellel logistic
    
end



