С помощью факторизации Холецкого решить Ax = b, как решить для х

Question

A = [2 -0.3 -0.2;-0.3 2 -0.1;-0.2 -0.1 2];
b=[7; 5; 3;];
L=Cholesky(A);
disp(A)
disp(b)
disp(L)
x = A\b;
disp(x)

Используемый мной Cholesky:

function L=Cholesky(A)
L=A;
n=size(A,2);
for k=1:n
    L(k,k)=sqrt(L(k,k));
    for i=k+1:n
        L(i,k)=L(i,k)/L(k,k);
    end
    for j=k+1:n
        for i=k+1:n
            L(i,j)=L(i,j)-L(i,k)*L(j,k);
        end
    end
end
for i=1:n-1
    for j=i+1:n
        L(i,j)=0;
    end
end

это код, который я пишу, но L не соответствует x

Это проблема, которую я пытаюсь решить