Algoritmo de Choleski

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function [R] = choleski(A)
 n = length(A(:,1));
 Av = A(2:n,2:n);
 if n == 1 then
 R = sqrt(A);
 else
 R(1,1) = sqrt(A(1,1));
 r = 1/R(1,1)*A(2:n,1);
 Rv = choleski(Av - r*r');
 for j = 2:n
 R(1,j) = r(j-1);
 R(j,1) = 0;
 for i = 2:n
 R(j,i) = Rv(j-1,i-1);
 end
 end
 end
endfunction