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

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

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