>>588,595
p が固定されてて、q が球面上を動くのか
じゃあ、r=(x,y,z), p=(p_x,p_y,p_z) として、同じ計算で
div((r-p)/|r-p|^3) = 4πδ(r-p)
(>>592 はδ関数の中身も間違えてる)

V を S の内部として
∫[S] (q-p)/|q-p|^3・nds
= ∫[V] div((r-p)/|r-p|^3) d^3r
= ∫[V] 4πδ(r-p) d^3r
= 4π