%this routine solves the system of odes for an SIR epidemic

clear
global gamma kr N

%specify parameter values
gamma = 1.275;
gamma = 2.01;
kr = 2;
N = 2;


%specify the temporal output points
tspan = [0:.1:50];

%specify initial data

y0 = [0.999,.001];

[T,Y] = ode23('SIR_deRHS',tspan, y0);  
figure(1)
plot(Y(:,1),Y(:,2))
figure(2)
plot(T,Y)