LU分解函数

 function [L,U]=mylu(A)
[m,n]=size(A);
if m~=n
    display('请输入方阵');
end
if det(A)==0
    display('输入错误 请输入可逆矩阵');
end
for i=1:n-1
   for j=1: n
        l(j,j,i)=1;  
    end
end
B(:,:,1)=A;
for i=1:n-1
   for j=i: n
     l(j,i,i)=B(j,i,i)/B(i,i,i);
 end
 B(:,:,i+1)=inv(l(:,:,i))*B(:,:,i);
end
L=eye(size(A));
for i=1:n-1
    L=L*l(:,:,i);
end
for i=1:n-1
    U=inv(l(:,:,i))*B(:,:,i);
end