% 改進平方根法解對稱系數陣線性方程組
function b=MRoot_squaring(A,b,n)
% 求解LR分解的各元素，并以重址方式存儲在矩陣A中
for i=2:n
    for j=i:n
        S=A(i,j);
        for k=1:i-1
            S=S-A(k,i)*A(k,j)/A(k,k);
        end
        A(i,j)=S;
    end
end
% 解Ly=b  y也存放在b里
for i=2:n
    S=b(i,1);
    for j=1:i-1
        S=S-A(j,i)*b(j)/A(j,j);
    end
    b(i)=S;
end
%解Rx=y
b(n)=b(n)/A(n,n);
for i=n-1:-1:1
    S=b(i)/A(i,i);
    for j=i+1:n
        S=S-A(i,j)*b(j)/A(i,i);
    end
    b(i)=S;
end
return