alg1
10
alg2
10
r1
1.3
r2
1.3
k1
2000
k2
2000
c1
0.5
c2
1












pt
alg1
r1*pt*(k1-c1*qt-pt)/k1
qt
alg2
r2*qt*(k2-c2*pt-qt)/k2