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=1/del*(del*ni-del*kc-bc*c+sqrt(del*del*ni*ni+2.0*del*del*ni*kc-2.0*del*ni*bc*c+del*del*kc*kc+2.0*del*kc*bc*c+bc*bc*c*c))/2.0
c'=fc(n)*c-fb(c)*r/eps-del*c
r'=fb(c)*r-(del+m)*r
init c=1,r=1
done