function [L,d,x]=cholesky(A,b)
[n,m]=size(A);
k=1;
if n==m
    for i=1:n;%判斷是否為對稱矩陣
        for j=1:n;
            if A(i,j)~=A(j,i)
                k=0;
                disp('此矩陣不是對稱矩陣');
                break;
            end
        end
        if k==0
            break;
        end
    end
    if k~=0%判斷是否為正定矩陣,正定矩陣的充要條件為A為對稱矩陣,且所有的順序主子式大于零
        for i=1:n;
            B=A(1:i,1:i);
            if det(B)<=0
                disp('此矩陣的順序主子式不全大于零,不為正定據陣');
                k=0;
            end 
        end
    end
    if k~=0%如果A為正定矩陣,則進行下一步操作
        L=diag(ones(1,n));
        for i=1:n
            for j=1:i-1
                s=0;
                for l=1:j-1
                   s=s+d(l)*L(i,l)*L(j,l);
                end
                L(i,j)=(A(i,j)-s)/d(j);
            end
            s=0;
            for l=1:i-1
                s=s+d(l)*L(i,l)^2;
            end
            d(i)=A(i,i)-s;
        end
        d=diag(d);
        y(1)=b(1);
        for i=2:n
            y(i)=b(i)-L(i,1:i-1)*y(1:i-1)';
        end
        B=inv(d)*y';
        C=L';
        x(n)=B(n);
        for i=n-1:-1:1
            x(i)=B(i)-C(i,i+1:n)*x(i+1:n)';
        end
    end
else
    disp('此矩陣不為方陣');
end