自编Python程序实现数值计算矩形区域二重积分，使用复化辛普森法。以函数f=xsiny在0<x<1和pi/2<y<pi区域上的积分为例。网格节点数m，n需为2的倍数。对于非矩形区域可以使用虚拟节点和区域，填补为矩形区域后计算，填补区域上的函数值为0.