在柱面坐标系下，点电荷的相互作用能是ｚ和ｐ的具有极点的二重积分，此类积分有简单的解析解，但不容易数值计算。Ｍａｔｈｅｍａｔｉｃａ３．０和ＭａｐｌｅⅤ不能求解此类积分；ＧａｕｓｓＬｅｇｅｎｄｒｅ算法......