#The differential equations

S'=mu*(1-S)-beta*S*Y
Y'=beta*S*Y-(mu+gamma)*Y
R'=gamma*Y-mu*R

#These are comments they mean nothing

#initial conditions
init S=.99 Y=.01 R=0

#parameters
par mu=0.00005 beta=0.5 gamma=0.14

#extra crap
@ dt=0.02 total=100 bounds=200000 maxstor=200000

done