双重积分问题

huahekulang

想请教关于二重积分问题。我写出了程序。
【对于 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 的积分，该怎么编程呢？ 想了很久，还是没有答案。