function rhs = Fode12Wb(t,y)

%
%  save file as Fode12Wb.m
%     - right-side forcing for ode45 solver
%

global ODE r1 a1 r2 a2 s1 s2

rhs = zeros(size(y));

switch ODE
    case{1}
        rhs(1) =  r1*y(1)*(1 - a1*y(1)) - s1*y(1)*y(2);
        rhs(2) =  r2*y(2)*(1 - a2*y(2)) - s2*y(1)*y(2);
    case{2}
        rhs(1) =  y(1)*(r1 - a1*y(1) - s1*y(2));
        rhs(2) =  y(2)*(r2 - a2*y(1) - s2*y(2));        
end