马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
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
[email=fH3=@(x)(Cnx(x,m,T0,lamb,Gamma,pjun).*dBnxdx(x,n,T0,lamb,Gamma,pjun]fH3=@(x)(Cnx(x,m,T0,lamb,Gamma,pjun).*dBnxdx(x,n,T0,lamb,Gamma,pjun[/email])+...
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 |