Cholesky分解法的matlab程序

function[L,flag]=chol_factor(A)n=length(A);L=zeros(n);flag='OK';fork=1:ndelta=A(k,k);forj=1:k-1delta=delta-L(k,j)^2;endifdelta1e-10flag='failure';return;endL(k,k)=sqrt(delta);fori=k+1:nL(i,k)=A(i,k);forj=1:k-1L(i,k)=L(i,k)-L(i,j)*L(k,j);endL(i,k)=L(i,k)/L(k,k);endend