2018A高温作业服设计 (显式 隐式差分代码)(13)

B(LEN4,1)=2*dt*h(2)*Tin/(rho(4)*c(4)*dx(4));
%以下是重复出现的数，先计算方便书写
F=lam./c./rho.*dt./dx./dx;

D=[dt/(dx(1)*rho(1)*c(1) dx(2)*rho(2)*c(2)),dt/(dx(3)*rho(3)*c(3) dx(2)*rho(2)*c(2)),dt/(dx(3)*rho(3)*c(3) dx(4)*rho(4)*c(4))];