不知道错在哪一步怎么办？

zophim

用matlab的两个函数trapz、quad2d分别求二重定积分的结果相差近10倍，
trapzq求解的速度快想用它代替掉quad2d。但是发现两个函数积分结果不一样。可能是tarapz那出错了，但不知道错在哪一步。请教大神，非常感谢！
x = 0.1:0.01:0.5;
" `2 k5 s: w2 r0 d9 Uy= 0.3:0.001:1;
/ b* R3 G4 k4 L  X' e[X,Y] = meshgrid(x,y);
F = cos((pi*Y).*(sin(X*pi)));
I = trapz(y,trapz(x,F,2));
2 T9 y% D3 G4 {" {. ]4 M
" C: K& _" O9 O- qf =@(x,y)cos((y).*(sin(x)));
5 [2 \- _! `& x' L% M3 I" H/ _) _$ V9 Aa = quad2d(f,0.1*pi,0.5*pi,0.3*pi,pi);
& N" m* d3 t. P/ O" e* [1 r% {

Allevi

可以随便更改 积分区间？简单放缩可以？
7 e+ e! W$ `& ~" k+ H% `假设 f(x)=x,[0,2],可以简单化为  f(x)=2*x,[0,1]？？？？

a=quad2d(@(X,Y)cos((pi*Y).*sin(X*pi)),0.1,0.5,0.3,1)
* U+ i; X) K! f1 b! g0 C) H
3 |- c$ a/ u" Q. Ga = 0.0102
: {$ \( r9 h7 A/ G/ a2 v# E: k与I接近。