function dy = romio(t,y)

global a b c d

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

dy(1) = a*y(1) + b*y(2);
dy(2) = c*y(1) + d*y(2);

