s'=-b*i*s
i'=b*i*s-g*i
par b=1,g=1
@ xp=s,yp=i
@ xlo=-.5,ylo=-.5,xhi=2,yhi=2
@ nmesh=200
done