function F_prime=deRHS(t,y)
global gamma kr N

s = y(1);
i = y(2);

Sp = -gamma*i*s+kr*(N-i-s);
Ip = gamma*i*s -i;

F_prime = [Sp;Ip];