matlab对矩阵进行LDLT和cholesky分解

%**************************************************************************
%（第四题）求矩阵的LDLT分解及cholesky分解--3111312053龚小林 硕1080
%**************************************************************************
clear
clc
format short
%待分解计算矩阵A
for i=1:10
    for j=1:10
        if j==i
            A（ij）=i;
        else
            A（ij）=min（ij）-2;
        end
    end
end
disp（‘待分解矩阵A为‘）
A
AG=A;
%矩阵LDLT分解
for k=1:10
    for p=1:k-1
        r（p）=A（pp）*A（kp）;
    end
    sum=0;
    for p=1:k-1
        sum=sum+A（kp）*r（p）;  
    end
        A（kk）=A（kk）-sum;
    if A（kk）==0
        disp（‘错误‘）
        break;
    end
    D（k