你这就是解2元非线性方程组。phi[x_, y_] := (m*x + n*x^2)*(p*y + q*y^2)*(x - y)/a*(m*x + n*x^2) + b*(p*y + q*y^2)In[3]:= Solve[{D[phi[x, y], x] == 0, D[phi[x, y], y] == 0}, {x, y}]Out[3]= {{x -> 0, y -> -(p/(2 q))}, <<1>>, <<17>>,