function rhs = vdp_period_f(t, x, mu)
% odefile for vdp period finder code (to be called by ode45)
%
% 2003-11-07 cbm
% no rights reserved: public domain code
  
rhs = [x(2); -mu*(x(1)^2 - 1)*x(2) - x(1)];

