Th Cholesky code doesn't work for symbolic cases, I think. Try this (loading linearalgebra isn't needed):
(%i21) cholesky(matrix([14,32,596],[32,77,1217],[596,1217,36245]));
(%o21) matrix([sqrt(14),0,0],[32/sqrt(14),3^(3/2)/sqrt(7),0],[596/sqrt(14),-(113*sqrt(3))/sqrt(7),5*6^(3/2)])
(%i22) %.transpose(%);
(%o22) matrix([14,32,596],[32,77,1217],[596,1217,36245])
--Barton
________________________________
T
So, as you can see, the cholesky function is not loaded. This is confirmed by the simple example I introduced above:
(%i3) cholesky(matrix([1, p12, p13], [p12, 1, p23], [p13, p23, 1]));
Unable to find the Cholesky factorization
-- an error. To debug this try: debugmode(true);
Actually, what this shows is that the cholesky function IS loaded, since it gives an error message.
It was unable to find the factorization though. I have not looked at your matrix at all.
________________________________
This email is subject to a disclaimer, please click on the following link or cut and paste the link into the address bar of your browser.
https://www.wellsfargo.com/com/disclaimer/wfd