# Goodwin oscillator

params trans=1,tl=1
param degM=1,degPo=0.5,degPt=1
params conv=0.5,n=10,unbin=10,bin=100

M'=  (trans*(1/((unbin/bin)+Pt))^n)-degM*M
Po'= tl*M-degPo*Po-conv*Po
Pt'= conv*Po-degPt*Pt

@ TOTAL=30,DT=.05,xlo=0,xhi=5,ylo=0,yhi=5
@ xplot=M,yplot=Po
init M=1,Po=1,Pt=2

done