#	Math865L_example7.ode (XPPAUT file) written by Chiu-Yen Kao

#	This is an XPP file to numerically integrate the equations for SIR model


#	SYSTEM OF EQUATIONS
	
	dS/dt =  PI-mu*S-beta*S*I
	dE/dt = beta*S*I-(mu+k)*E	
	dI/dt =  k*E-(mu+alpha)*I

#	PARAMETERS

	parameter beta=0.001, alpha=0.1, k=0.1, mu=0.1


#	INITIAL CONDITIONS

	init S=499,E=200, I=1.0


#	CHANGES FROM XPP'S DEFAULT VALUES

@ TOTAL=20.0,DT=0.01,XLO=0.0,XHI=20.0,YLO=0,YHI=1000s

	done