function step = Xvdp(t,y)

%
%  like the nonlinear matlab example problem
%

step = zeros(size(y));

step(1) = y(2);
step(2) = 100*(1 - y(1)^2)*y(2) - y(1);