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

global a b 

 xdot=y;
 ydot=-x-a*(x^2-1)*y;
% xdot=x+exp(-y);
% ydot=-y;