>>588
だいたい合ってる
p=q のところに特異性があるから、それを考えると
div((q-p)/|q-p|^3) = -4πδ(q)
(これもガウスの定理から出る)

これを体積積分すれば最初の面積分が出るけど、
これ本当にガウスの定理を使えって問題？