点电荷相互作用能与一类二重积分的算法

在柱面坐标系下，点电荷的相互作用能是ｚ和ρ的具有极点的二重积分。此类积分有简单的解析解，但不容易数值计算。Ｍａｔｈｅｍａｔｉｃａ３．０和ＭａｐｌｅＶ不能求解此类积分；Ｇａｕｓｓ－Ｌｅｇｅｎｄｒｅ算法只能获得粗略结果，Ｒｕｎｇｅ－Ｋｕｔａ算法虽能成功地计算，但不能交换积分顺序。