function f=function_2d(x,t)

global a e om

x1=x(1);x2=x(2);

f1=x2;
f2=-om^2*x1-2*e*x2;%e*(x1^2-1)*x2;

f=[f1;f2];