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

zophim

用matlab的两个函数trapz、quad2d分别求二重定积分的结果相差近10倍，
6 P* z) |% T& P0 W1 |- otrapzq求解的速度快想用它代替掉quad2d。但是发现两个函数积分结果不一样。可能是tarapz那出错了，但不知道错在哪一步。请教大神，非常感谢！
x = 0.1:0.01:0.5;
; l+ ^4 v; i7 k& P- cy= 0.3:0.001:1;
4 Q; e! Y- Y( H+ \+ B# q% ^% Z/ l[X,Y] = meshgrid(x,y);
F = cos((pi*Y).*(sin(X*pi)));
I = trapz(y,trapz(x,F,2));
" I( ^, x9 s) v/ z% [0 s
4 c) @1 ^3 X" z  gf =@(x,y)cos((y).*(sin(x)));
6 D- ]! C0 k+ B% ]6 s! Z; D' na = quad2d(f,0.1*pi,0.5*pi,0.3*pi,pi);
( Z, o+ ~' v. S: V2 h8 c

Allevi

可以随便更改 积分区间？简单放缩可以？
0 S' W( h1 p) C! i3 ]. x8 Q假设 f(x)=x,[0,2],可以简单化为  f(x)=2*x,[0,1]？？？？
* W; T  Y" e! d, W' M
2 ~6 z* b8 w8 Ka=quad2d(@(X,Y)cos((pi*Y).*sin(X*pi)),0.1,0.5,0.3,1)
4 C4 q) Z) l# [+ J4 o
a = 0.0102
3 ^+ m  P. L, M+ X7 o7 m与I接近。