多重积分避免循环

[cross]

我目前处理如下积分，
int（dx1 dy1 dz1 （int（dx2 dy2 dz2 （D））））
其中D为三维空间点阵中任意两点 （x1，y1，z1）和（x2，y2，z2）之间得距离。
对于数值积分我采用simpson法。流程如下：
1。任意取第一点（x1，y1，z1），求其对空间全部点的距离，并对dx2 dy2 dz2积分
2。跑遍全部的（x1，y1，z1），得到上一次积分的结果三维矩阵，最后对 dx1 dy1 dz1

积分

目前方法的问题在于，跑遍全部 （x1，y1，z1）的时候需要 对3个坐标作3重循环，大量

占用资源。 我没有想出办法能用矩阵运算代替大量循环，特此求教。谢谢。