function [ xdot , ydot] = func2D(x,y,t)

global gamma om;


% talantotis me aposvesi
xdot =  y;
ydot = -2*gamma*y - om^2*x;

% van der poll
% xdot =  y;
% ydot = gamma*(1-x^2)*y - om^2*x;