alg1
900
alg2
500
r1
1.1
r2
1.05
k1
8000
k2
16000
c1
1
c2
1












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