function L=chol6245(R)
% function L=chol6245(R)
%
% fault-less version of chol.m 
[L,r]=chol(R);
if r>0,
    [V,D]=eig(R);
    L=V*diag(sqrt(abs(diag(D))))*V';
end