用于二维域的快速、完全矢量化的Simpson方法版本。该代码避免使用任何for循环等,可在给定精度级别下比dblquad快一个数量级或更快。代码示例提供了如何使用。功能: ans = simp2D('func',xs,xe,ys,ye,NX,NY)输入参数: func - 接受向量输入的二维函数(否则可能导致错误结果) xs, xe - x积分极限 ys, ye - y积分极限 NX - x方向的积分区间数(应为偶数) NY - y方向的积分区间数(应为偶数)