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

zophim

用matlab的两个函数trapz、quad2d分别求二重定积分的结果相差近10倍，
/ a! V6 m* w8 y) S' c3 T3 o4 ctrapzq求解的速度快想用它代替掉quad2d。但是发现两个函数积分结果不一样。可能是tarapz那出错了，但不知道错在哪一步。请教大神，非常感谢！
x = 0.1:0.01:0.5;
y= 0.3:0.001:1;
, @" g# Z, {4 s6 H# m" z[X,Y] = meshgrid(x,y);
" W$ _/ n! T* G* o3 p; T, w) QF = cos((pi*Y).*(sin(X*pi)));
; K) ~6 l; l+ Q+ `: y) C3 |- G: RI = trapz(y,trapz(x,F,2));

f =@(x,y)cos((y).*(sin(x)));
& E- Z: y. C: }  [0 `" y* Ia = quad2d(f,0.1*pi,0.5*pi,0.3*pi,pi);

Allevi

可以随便更改 积分区间？简单放缩可以？
假设 f(x)=x,[0,2],可以简单化为  f(x)=2*x,[0,1]？？？？
5 s( a7 _: a/ G4 T2 c" n
. G7 K, }7 Z5 F) e3 r3 _a=quad2d(@(X,Y)cos((pi*Y).*sin(X*pi)),0.1,0.5,0.3,1)
/ r! Y# g" \0 Y3 S- Z* @$ e
4 C) P: t% J# i5 h$ P. f% p4 a; Na = 0.0102
与I接近。