findpass 发表于 2009-12-10 16:00

matlab积分问题

function H3jmn=H3jmn(T0,lamb,Gamma,pjun,N)
H3jmn=zeros(N,N,N);
for jj=1:N
    for n=1:N
      for m=1:N
            fH3=@(x)(Cnx(x,m,T0,lamb,Gamma,pjun).*dBnxdx(x,n,T0,lamb,Gamma,pjun)+...
                Bnx(x,n,T0,lamb,Gamma,pjun).*dCnxdx(x,m,T0,lamb,Gamma,pjun)).*...
                WB(x,T0,lamb).*Bnx(x,jj,T0,lamb,Gamma,pjun);
            H3jmn(m,n,jj)=double(quadl(fH3,0,1));
      end
    end
end

就是T0,lamb,Gamma,pjun,N是参数,然后计算H3jmn是个三维数组,就是一定要(m,n,j)位置上分别对应值,怎么才能把这个for的循环去掉呢?:@L
页: [1]
查看完整版本: matlab积分问题