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

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

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