alg1
900
alg2
500
r1
1.2
r2
1.1
k1
8000
k2
16000
c1
2
c2
2












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