A=[6,7,5;7,13,8;5,8,6];
b=[9 10 9]';
n=length(b);%方程个数n
G=zeros(n,n);
G(1,1)=sqrt(A(1,1));
G(2:n,1)=A(2:n,1)/G(1,1);
for j=2:n-1
    G(j,j)=sqrt(A(j,j)-sum(G(j,1:j-1).^2));
    for i=j+1:n
        G(i,j)=(A(i,j)-sum(G(i,1:j-1).*G(j,1:j-1)))/G(j,j);
    end
end
G(n,n)=sqrt(A(n,n)-sum((G(n,1:n-1)).^2));
G
%--------用A=GG’分解求解方程组Ax=b（平方根法）----
%Gy=b，G'x=y;
x=zeros(n,1);%未知向量
y=zeros(n,1);%中间向量
y(1)=b(1)/G(1,1);
for i=2:n
    y(i)=(b(i)-sum(G(i,1:i-1)'.*y(1:i-1)))/G(i,i);
end
y
%---------由G'x=y求出x----------
%方法类似于U由x=y解出x
G=G';
x(n)=y(n)/G(n,n);
for i=n-1:-1:1
    x(i)=(y(i)-sum(G(i,i+1:n)'.*x(1:i+1)))/G(i,i);
end
x


用平方根法求解方程组，可是求解结果x的值不对，求帮忙，先谢谢了！

我已经找出来了，不过还是谢谢你了！

x(i)=(y(i)-sum(G(i,i+1:n)'.*x(1:i+1)))/G(i,i);改成
x(i)=(y(i)-sum(G(i,i+1:n)'.*x(i+1:n)))/G(i,i);