# Blood cell dynamics
#

dx/dt = a*y/(1+y^2) -x
dy/dt = x-b*y

par a=2, b=1

done