function xdot = vdpfun(t,x)

global e;

x1dot = x(2);
x2dot = -x(1)+e*(1-x(1)^2)*x(2);
xdot = [x1dot x2dot]';

