function x=Chol_Solve(A,b)
%利用對稱正定矩陣之Cholesky分解求解線性方程組Ax=b
n=length(b);
l=Cholesky(A);
x=ones(1,n);
y=ones(1,n);
for i=1:n
    z=0;
    for k=1:i-1
        z=z+l(i,k)*y(k);
    end
    y(i)=(b(i)-z)/l(i,i);
end
for i=n:-1:1
    z=0;
    for k=i+1:n
        z=z+l(k,i)*x(k);
    end
    x(i)=(y(i)-z)/l(i,i);
end