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

zophim

用matlab的两个函数trapz、quad2d分别求二重定积分的结果相差近10倍，
4 g% |+ V; H3 |0 I- B$ |( ttrapzq求解的速度快想用它代替掉quad2d。但是发现两个函数积分结果不一样。可能是tarapz那出错了，但不知道错在哪一步。请教大神，非常感谢！
x = 0.1:0.01:0.5;
0 b, S7 Z: Z" `5 Y3 u* s" dy= 0.3:0.001:1;
% B* ^: f8 V0 {0 ?[X,Y] = meshgrid(x,y);
" Q- M0 s, `% t/ x+ SF = cos((pi*Y).*(sin(X*pi)));
I = trapz(y,trapz(x,F,2));
% H5 _5 M6 q5 e% _) A5 A
3 O2 Q+ }; ^" [4 Ef =@(x,y)cos((y).*(sin(x)));
a = quad2d(f,0.1*pi,0.5*pi,0.3*pi,pi);
5 Q# n  }+ T" A, Y5 @7 `

Allevi

可以随便更改 积分区间？简单放缩可以？
: v; C  C5 J% R4 P1 J假设 f(x)=x,[0,2],可以简单化为  f(x)=2*x,[0,1]？？？？
4 b( M  G5 Y7 L# K
. s" q8 r& J5 R. ^& Qa=quad2d(@(X,Y)cos((pi*Y).*sin(X*pi)),0.1,0.5,0.3,1)

a = 0.0102
与I接近。