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

A(n,n-1)=-F(2);A(n,n)=2*F(2);A(n,n 1)=-F(2);
end
%第二段与第三段接触点的位置
A(LEN2,LEN2-1)=-2*D(2)*lam(2)/dx(2);

A(LEN2,LEN2)=2*D(2)*(lam(2)/dx(2) lam(3)/dx(3));
A(LEN2,LEN2 1)=-2*D(2)*lam(3)/dx(3);
for n=LEN2 1:LEN3-1