q1'=q1+beta1*q1*(M-c-q2-2*q1)
q2'=q2+beta2*q2*(M-c-q1-2*q2)
init q1=1, q2=3
par   beta1=0.1, beta2=0.2, M=10, c=5
@ maxstor=100000, meth=discrete
done