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

% mms $Id: $

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