clc, clear all, close all

%*** first matrix
A=[1,-1,0;-1,1,1;0,-1,2]
L=cholesky(A);
L*L'
pause
%*** second matrix
disp('**************************')
A=[1,2,1;2,1,1;1,1,2]
L=cholesky(A);
L*L'
pause
%*** third matrix
disp('**************************')
A=[1,2,1;2,8,2;1,2,2]
L=cholesky(A);
L*L'
