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

zophim

用matlab的两个函数trapz、quad2d分别求二重定积分的结果相差近10倍，
trapzq求解的速度快想用它代替掉quad2d。但是发现两个函数积分结果不一样。可能是tarapz那出错了，但不知道错在哪一步。请教大神，非常感谢！
, `2 J3 ~0 W& T$ m/ Kx = 0.1:0.01:0.5;
y= 0.3:0.001:1;
[X,Y] = meshgrid(x,y);
F = cos((pi*Y).*(sin(X*pi)));
7 ~  i% }3 w$ f3 w$ i+ M1 |I = trapz(y,trapz(x,F,2));

9 f0 R/ O9 W/ ef =@(x,y)cos((y).*(sin(x)));
5 M& V" U( H' ?2 H7 o" t1 ra = quad2d(f,0.1*pi,0.5*pi,0.3*pi,pi);
% T3 b# z$ \7 O- e! [  F

Allevi

可以随便更改 积分区间？简单放缩可以？
假设 f(x)=x,[0,2],可以简单化为  f(x)=2*x,[0,1]？？？？
) S  I& Q* r0 P- i/ Z) ~% X; n
1 v2 x0 L/ S  Z. Ya=quad2d(@(X,Y)cos((pi*Y).*sin(X*pi)),0.1,0.5,0.3,1)

6 M, d% n' q& m- }a = 0.0102
与I接近。