# Chemostat
#

dx/dt = -r*xo*x*y/(1+x) + xo - x
dy/dt = r*x*y/(1+x) - y

par xo=1, r=.5

done