function dV=SIR(t,V,alpha,beta)

S=V(1);
I=V(2);

dS=-alpha*S*I;
dI=alpha*S*I-beta*I;

dV=[dS;dI];
end