function yprime = vdpol(t,y,FLAG,mu,omegasq);
%  Vdpol(t,y) returns the state derivatives of the Van der Pol equation
%
yprime = [y(2); mu*y(2)*(1.0 - y(1)^2)-omegasq*y(1)];