用于 3D 域的 Simpson 方法的快速、完全矢量化版本。 此代码避免使用任何 for 循环等。对于给定的精度水平，它可以比三重四边形快一个数量级或更快。 代码中提供了示例。 用途： ans = simp3D('func3D',xs,xe,ys,ye,zs,ze,NX,NY,NZ) 输入参数： func3D - 接受向量输入的 3D 函数（否则可能会出现误导性结果）
xs, xe - x 积分极限ys, ye - 积分的 y 极限zs, ze - z 积分极限
NX - x 方向的积分区间数（应该是偶数） NY - y 方向的积分区间数（应该是偶数） NZ - z 方向的积分区间数（应该是偶数）
