[A,B]=matrix();
[r,L]=size(A);
j=1;
for n=1:L
i=1;
for m=1:r-n;
i=i+1;
A(n+m,:)=A(n+m,:)-(A(n+m,n)./A(n,n)).*A(j,:);
if n+m==r
j=j+1;
end
if i>=L
break
end
end
end
上三角矩阵生成法.zip （4个子文件） 
Run_main.m 237B matrix.asv 151B Run_main.asv 242B matrix.m 151B
