q1'=beta1*q1*(M-c-q2-2*q1)
q2'=beta2*q2*(M-c-q1-2*q2)


par beta1=1 beta2=1 M=20 c=1

@ xp=q1,yp=q2,xlo=0,ylo=0,xhi=25,yhi=25
done