alg1
900
alg2
400
r1
0.01
r2
0.005
k1
1500
k2
3000
c1
1
c2
2












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