mathematica二重数值积分

被积函数太过复杂，而且积分区域是曲形区域，也就是变上下限，用传统的NIntegrate求不出来数值解，网上用复合辛普森方法算的都是矩形区域，这个到底怎么求。