Added Cholesky factorisation and equation-solving to linalg. Modified linalg.diag to extract entries of main-diagonal if input array is 2D