M = [5 1.2 0.3 -0.6; 1.2 6 -0.4 0.9; 0.3 -0.4 8 1.7; -0.6 0.9 1.7 10];
b=rand(4,1);
[x]=chelosky_dejen_2018(M,b);