function ydot=vdpol(t,y)
mu = 2; % set to a value greater than zero
y1dot = y(2);
y2dot = mu*(1-y(1)^2)*y(2)-y(1);
ydot = [y1dot; y2dot];