双重积分问题
想请教关于二重积分问题。我写出了程序。【对于 r 和 z1 进行双重积分】
function X = integ(ri,ro,ni,y,z,L1,L2,n)
dz1=(L2-L1)/(n-1);
dr1=(ro-ri)/(ni-1);
X=0;
for k=1:n-1
z1=L1+dz1*(k-1);
z2=L1+dz1*k;
f1=I1(r,y,z,z1); % I1 is function of r y z z1
f2=I1(r,y,z,z2);
F=0.5*(f1+f2)*dz1;
X=X+F;
end
上面的程序只是对 z1 进行了积分。对于 r 的积分，该怎么编程呢？ 想了很久，还是没有答案。
页:
[1]