问一下，我们老师要求做一个点电荷的电势分布，不能直接用积分符号做，本来我以为挺简单，直接用梯形积分就行了，再把点连起来就好了，结果发现一个点就要从x积分到无穷远，根本不能算啊，更何况那么多点

看具体函数了
估计有极限

1. 单个点电荷电势
f[k_, Q_, r_] := k Q/r
Show[VectorPlot3D[{x, y, z}, {x, -10, 10}, {y, -10, 10}, {z, -10,
10}], SphericalPlot3D[
Evaluate@Table[f[8.9, 1, r], {r, 1, 3}], {\[Theta], 0, Pi}, {\[Phi],
0, 3 Pi/2}], PlotRange -> All]
2. 两个点电荷的电势
electroStaticPotential[q_, p_, r_] :=
Sum[ q[[i]]/Norm[r - p[[i]]], {i, Length[q]}]
electricField[{q1_,
q2_}, {{x1_, y1_, z1_}, {x2_, y2_, z2_}}] = -D[
electroStaticPotential[{q1,
q2}, {{x1, y1, z1}, {x2, y2, z2}}, {x, y, z}], {{x, y, z}}] /.
Derivative[1][Abs][x_] :> x/Sqrt[x^2];
c = ContourPlot3D[
Evaluate[electroStaticPotential[{1, -1}, {{-1, 0, 0}, {1, 0,
0}}, {x, y, z}]], {x, -6, 6}, {y, 0, 4}, {z, -4, 4},
Contours -> {-0.75, -0.25, -0.1, 0, 0.1, 0.25, 0.75},
ContourStyle -> Table[Hue[i/7, 1, 1, 0.5], {i, 0, 6}], Mesh -> None]
v = VectorPlot3D[
electricField[{1, -1}, {{-1, 0, 0}, {1, 0, 0}}], {x, -6, 6}, {y, 0,
4}, {z, -4, 4}, VectorStyle -> "Arrow3D", VectorPoints -> 5,
VectorScale -> {.3, Scaled[0.3]}]
Show[c, v]

插一句。你们老师提的这个问题其实是避重就轻的，因为点荷电势在空间上的衰减非常快所以只要适当截断就可以积了，但更一般的无穷域上的积分其实没那么简单，参看自带帮助中的教程《NIntegrate 积分策略》（NIntegrateIntegrationStrategies）