function dy=sircovf(t,y)
global N beta gamma

dy(1,1)=-beta/N*y(1).*y(2);
dy(2,1)=beta/N*y(1).*y(2)-gamma*y(2);
dy(3,1)=gamma*y(2);