1.pv=(q1^2+q2^2-q3^2)/[2(q1+q2-q3)]        q2=5.3415是q 取值的上限
2.pv+(1-2a)q+2F=2(1-a)C(q)       at q=q1 and q=q3。

F,a是parameters,这里先取F=0.015641,a=0.25
C(q)是一个convex 的cost function，这里先用ln(q+1)

把第一个是公式带到第二个里面，给定F,a，C(q)就可以解q1,q3了。
我用的mathematica,也是初学，没有解出来，下面是我的input和output

不理解system给的红字。。。。
matlab也一知半解，求大神帮忙解出来这两个q，q的上限q2是5.3415，所以希望q1，q3<5.3415,可以调试F,0<a<1,C(q) 来找到合适的解。
论坛币好商量！