# epidemic with periodic forcing
s'=-r*s*i+del*(1-s)
i'=r*s*i-gam*i
par del=.1,gam=.2
r=r0*(1+a*u)
u'=u*(1-u*u-v*v)-w*v
v'=v*(1-u*u-v*v)+w*u
init u=1,s=.196,i=.405
par r0=1,a=.2,w=1
@ total=6.28,dt=.02
done