I cant get the correct answer through this code, and I cannot find out where is my mistake.

Also, how to create matrices A = (aij ) with aij = 1/(i+j1) with n = 3,4,5.

= cholesky (A) function L = cholesky(A) [~,n]=size(A); L = A;] for k = 1:n for i = 1:k-1 S = 0; for j = 1:1-1 s + L(i,j)*L(k,j); end L(k, i) = (L(k,i)-s)/L(i,i); end V = 0; for j = 1:k-1 V = V + L(k, j)^2; end L(k,k) = sqrt(L(k, k)-v); end end function CHOLESKY-FACTORIZATION(C) Factors C = LLT, assuming C is symmetric and positive definite LAC > This algorithm destructively replaces C with L for k + 1, 2, ...,n D Back-substitute to place at the beginning of row k for it 1,...,k - 1 Current element i of ik Sf0 = D Iterate over L11; j