### 5.50.1 Cholesky decomposition : cholesky

cholesky takes as argument a square symmetric positive definite matrix M of size n.
cholesky returns a symbolic or numeric matrix P. P is a lower triangular matrix such that :

tran(P)*P=M

Input :

cholesky([[1,1],[1,5]])

Output :

[[1,0],[1,2]]

Input :

cholesky([[3,1],[1,4]])

Output :

[[sqrt(3),0],[(sqrt(3))/3,(sqrt(33))/3]]

Input :

cholesky([[1,1],[1,4]])

Output :

[[1,0],[1,sqrt(3)]]

Warning If the matrix argument A is not a symmetric matrix, cholesky does not return an error, instead cholesky will use the symmetric matrix B of the the quadratic form q corresponding to the (non symmetric) bilinear form of the matrix A.
Input :

cholesky([[1,-1],[-1,4]])

or :

cholesky([[1,-3],[1,4]])

Output :

[[1,0],[-1,sqrt(3)]]