% factorization_cholesky.m
% Cholesky分解示例
% 產生n階pascal矩陣,pascal矩陣為pascal三角陣發展而來，為對稱正定陣
disp('正定對稱陣的分解：')
n=6;
X=pascal(n)
R=chol(X)
disp('非正定對稱陣的分解：')
X(n,n)=X(n,n)-1;
X
[R, p]=chol(X)
Xn_1=R'*R
disp('利用cholesky分解結果計算線性方程組')
A=pascal(6);
b=[1 3 5 7 9 4]';
x=A\b
R=chol(A);