% Esta funcion realiza la factorizacion de Cholesky de
% la matriz L haciendo un preordenamiento de la misma.
%
% Sintaxis:
%           [R,ord] = lchol(L);
% Donde:
%           R - Matriz de la factorizacion de Cholesky.
%         ord - Vector de permutaciones para la factorizacion
%               de Cholesky.
%
% Para mas informacion ver Readme.doc. 
%
% Autor: Nelson J. Trujillo Barreto.

function [R,ord] = lchol(L);
ord=symrcm(L);
R=chol(L(ord,ord));