function [L,R,X]=Doolittle(A,B) %%%输入n*n的方阵A, %%%通过Doolittle分解; %%%输出单位下三角矩阵L和上三角矩阵R。 [n,~]=size(A); L=eye(n); R=zeros(n); for k=1:n for j=k:n R(k,j)=A(k,j)-L(k,1:k-1)*R(1:k-1,j); end for i=k+1:n L(i,k)=(A(i,k)-L...