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