function dx = romiononlinear(t,y)

global alpha beta

dx = zeros(2,1);    % a column vector
%%%y(1) = Romio
%%%y(2) = Juliet

dx(1) = alpha*y(1);
dx(2) = -beta*(y(1)^2 - 1)*y(2) - y(1);
