function udot=ex4_ode(t,u,a)
% udot=ex4_ode(t,u,a)
% van der Pol ode with parameter a

% mms $Id: $

if nargin < 3
  a=1;
end

% udot must be a column vector
udot=[       u(2)
     -a*(u(1)^2-1)*u(2)-u(1)];