【正文】
for k=L1:1:1 %求解 x(k,1)=(c(k,L+1)c(k,k+1:L)*x(k+1:L))/c(k,k)。 c(i,:)=c(i,:)c(j,:)*m。 end c=[a b]。 b(q)=b(zhuh)。 a(zhuh,:)=p。 end end p=a(q,:)。for q=1:n big=max(max(abs(a(q:n,q:n))))。 L=length(b)。b=[1,0,0]39。 。 end end end x=y。 end y=[1:n(1)]39。 end end x(L,1)=c(L,L+