function yp = vanderpol(y)
% Returns the right hand sides of the Van der Pol first order system

% damping constant
c = 1;

% set the right hand sides
yp = [y(2); -c*(y(1)^2-1)*y(2)-y(1)];

endfunction
