% コレスキー分解を用いて連立一次方程式を解く
% 3次の正定対称行列 M から b = M*(与えたベクトル) を計算
% 行列 M をコレスキー分解し，連立一次方程式 Mx = b を解く．
% ベクトル x が与えたベクトルに一致することを確認する．

M = [ 4, -2, -2;...
     -2,  2,  0;...
     -2,  0,  3 ];
b = M*[ 1; 2; -1 ];

U = chol(M);
y = U'\b;
x = U\y;
x
