alg1
10
alg2
10
r1
1.4
r2
1.2
k1
1000
k2
2000
c1
1
c2
1












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