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


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

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