k^2 = (4 \[Rho]*r*Sin [\[Theta]])/(\[Rho]^2 + r^2 + 2 \[Rho]*r*Sin[\[Theta]]); V[r_, \[Theta]_] := Integrate[-\[Rho]/Sqrt[r^2 + \[Rho]^2 + 2 \[Rho]*r*Sin[\[Theta]]]* Integrate[1/Sqrt[ 1 - k^2*Sin[\[Psi]]^2], {\[Psi], -(\[Pi]/2), \[Pi]/2}], {r, 0, 1}];
2020年02月18日 15点02分