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

global e f;

%van der poll with forcing

xdot=y;
ydot=-x-e*(x^2-1)*y+f*cos(t);

end