看起來很簡單的數值積分
一維的積分我會轉成矩陣
但是重積分時好像沒有這麼簡單
請問各位有什麼方法
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