# Goodwin oscillator
 
dx1/dt=1/(1+x3^p)-b1*x1
dx2/dt=b2*(x1-x2)
dx3/dt=b3*(x2-x3)

p p=10, b1=.1, b2=.1, b3=.1

done