On the complex field

> B:=matrix([[3, 2+I, 3], [2-I, 0, 5+I], [3, 5-I, 0]]);

> C:=CHOLESKY(B,L,complejo);

Testing the output