function step = Xcubic(t,y)

%
%  like the nonlinear BVP problem
%

step = zeros(size(y));

step(1) = y(2);
step(2) = y(1) - y(1)^3;