fc(n)=bc*n/(kc+n)
fb(c)=bb*c/(kb+c)
par ni=90,del=.6
par bc=3.3,kc=4.3,bb=2.25,kb=15,m=.055,lam=.4,eps=.25
@ bound=10000,total=200
#n'=del*(ni-n)-fc(n)*c
n=del*ni/(del+bc*c)
c'=fc(n)*c-fb(c)*r/eps-del*c
r'=fb(c)*r-(del+m)*r
init c=1,r=1
done