function rhs = Fode12Wa(t,y)

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

global ODE r1 a1

rhs = zeros(size(y));

rhs(1) =  r1*y(1)*(1 - a1*y(1));
