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

dt=0.002;t=2000; %时间步长
Tout=75;Tin=37; %初始温度
LEN1=int8(x(1)/dx(1)) 1; %记录每段接触点的位置，外界与第一段接触点的位置

LEN2=LEN1 x(2)/dx(2);%第一段与第二段接触点的位置
LEN3=LEN2 x(3)/dx(3);%第二段与第三段接触点的位置
LEN4=LEN3 x(4)/dx(4);%第三段与第四段接触点的位置