需要求这个函数的临界点，求大神指导
h[q, r] = -3/(4*Pi)*
Log[2] + (2*Log[1 - r^2] + Log[1 - 2*(J - r^2)])/(4*
Pi) - (Log[((r^2)*(1 - Cos[2*q]))/(1 + r^4 -
2*r^2*Cos[2*q])] +
2*Log[(2*J - r^2 -
2*r*Sqrt[2*(J - r^2)]*
Cos[q])/(2*((1 + 2*r^2 (J - r^2) -
2*r*Sqrt[2*(J - r^2)]*Cos[q])))])/(4*Pi) /. J -> 0.456;
Solve[\!\(
\*SubscriptBox[\(\[PartialD]\), \(q\)]\(h[q, r]\)\) == 0 && \!\(
\*SubscriptBox[\(\[PartialD]\), \(r\)]\(h[q, r]\)\) == 0, {q, r}]

或者求近似解

首先，确保方程本身正确，然后，FindRoot慢慢试初值吧。

这个很简单，先这样然后那样，最后这样就解出来啦