# Bliss-Painter-Marr oscillator
 
dx1/dt=a/(1+x3)-b1*x1
dx2/dt=b1*x1-b2*x2
dx3/dt=b2*x2-c*x3/(K+x3)


p a=150, b1=.2, b2=.2, c=16, K=1

done