function xdot = vdpODE(t,x,mu) ;
xdot(1)=x(2) ;
xdot(2)=-x(1)+mu*(1-x(1)^2)*x(2) ;
%xdot=xdot( :) ;
