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

global e f;

%van der poll with forcing

xdot=y;
ydot=-x-e*x^3;

end