看起来很简单的数值积分
一维的积分我会转成矩阵
但是重积分时好像没有这么简单
请问各位有什么方法
for j1=1:101
for j2=1:101
q=0;
for j3=1:101
for j4=1:101
w=1/101*exp(j*k*z)*exp(j*k/2/z*(x0(j3)^2+y0(j4)^2))...
*exp(-j*2*pi/l/z*(x0(j3)*x1(j1)+y0(j4)*y1(j2)))*a(j3,j4);
q=q+w;
end
end
b(j1,j2)=q;
end
end