#Put the differential equations
S'=mu*(1-S)-beta*S*Y
E'=-mu*E-alpha*E+beta*S*Y
Y'=alpha*E-(mu+gamma)*Y
R'=gamma*y-mu*R
#Set the initial conditions
init S=.99 Y=.01 E=0 R=00.000
#Set parameters
par mu=0.00005 alpha=.25 beta=.5 gamma=.14
#Put all that extra crap
@ dt=0.02 total=100 bounds=200000 maxstor=200000
done