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
c'=fc(n)*c-fb(c)*b/eps-del*c
r'=fb(c)*r-(del+m+lam)*r
b'=fb(c)*r-(del+m)*b
init n=90,c=1,r=1,b=1
done