function fy = vanderPol2(x,y,mu)
%
% van der Pol equation
%
fy = [y(2);mu*(y(1)-1)*(y(1)+1)*y(2)-y(1)];

